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ABSTRACT 


Full  details  are  given  of  the  seismic  measurements  undertaken 
on  the  series  of  explosive  charges  detonated  at  Suffield  Experimental 
Station  in  the  period  195&-1963.  Seismograms  are  included  for  charges 
ranging  in  size  from  8  lb  to  200,000  lb  TNT,  surface  burst,  and  for  a 
variety  of  buried  charges. 

The  P  wave  is  analysed  only  to  the  extent  of  identification 
and  the  derivation  of  an  approximate  layering  system  for  the  test  site. 

The  surface  waves  are  analysed  in  some  detail  in  terms  of  their 
particle  trajectories.  In  the  radial- vertical  phase  the  seisnograms  from 
surface-burst  charges  are  divided  into  two  distinct  regions  by  the  arrival 
of  the  airblast.  Preceding  the  airblast  a  dispersed  surface  wave  is  evi¬ 
dent,  which,  in  general,  is  retrograde  in  habit  at  short  distances  and 
prograde  at  longer  distances  from  the  explosion.  After  the  arrival  of 
the  airblast  the  seisnograms  consist  of  a  series  of  large  amplitude  pulses, 
apparently  coupled  to  the  airblast,  upon  which  is  superimposed  an  air- 
coupled  wave  of  Press-Ewing  type. 

The  results  from  the  surface-burst  charges  are  compared  with 
those  from  buried  charges.  It  is  shown  that  the  buried  charges  generate 
a  succession  of  phases  which  alternate  in  habit  between  the  retrograde 
and  prograde  forms.  A  form  of  dispersion  analysis  similar  to  that  of 
Sato  is  described,  and  applied  to  certain  of  the  phases.  It  is  shown 
that  the  experimental  dispersion  curves  from  the  buried  shots  cross  the 
sonic  velocity,  thus  indicating  that  Press-Ewing  coupling  may  occur  to 
at  least  one  phase. 
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Theoretical  Rayleigh  Wave  dispersion  curves  are  computed 
by  the  Anderson  programme  for  models  which  approximate  to  the  test 
site.  It  is  shown  that  the  major  features  of  the  seismograms  can 
be  interpreted  in  terms  of  dispersion  curves  of  the  type  computed, 
but  that  a  full  analysis  will  require  detailed  analysis  by  similar 
methods  of  all  the  phases  detectable  in  the  records. 

It  is  shown  that  in  the  transverse  component  there  is 
repeatability  in  the  wave  form  in  a  given  azimuth  when  shots  are 
repeated  on  the  same  ground  zero,  and  that  there  is  some  form  of 
coupling  between  the  radial-vertical  modes  and  the  transverse  modes. 

The  data  given  provide  a  clear  picture  of  the  general 
pattern  to  be  expected  in  the  seismograms  from  very  large  detonations 
in  the  region  of  the  Western  Plains.  The  research  described  has  led 
to  the  publication  of  a  wide  range  of  reports  and  papers  by  the 
present  writer  and  others,  dealing  with  limited  aspects  of  the  problem, 
and  these  are  listed  or  quoted  at  the  relevant  points.  The  research 
has  developed  to  the  stage  where  further  interpretation  requires  a  long 
period  of  routine  analysis  of  the  spectral  content  and  the  dispersion 
characteristics  of  the  individual  phases. 

Extensions  of  the  experimental  programme  to  investigate 
special  aspects  of  the  problem  are  suggested. 

In  addition  to  the  strong  motion  seismograms,  some  experimental 
data  are  given  which  allow  the  correlation  of  the  ground  shock  and  airblast 
right  in  to  the  crater  area  to  be  demonstrated. 
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STRONG  MOTION  SEISMIC  EFFECTS  OF  THE  SUFFIELD  TNT  EXPLOSIONS 

by 

G.  He  S.  Jones 

1.00  INTRODUCTION 

1 .10  General  Background  to  the  Study 

1.11  Many  of  the  complexities  of  the  half-space  problem 
are  due  to  the  generation  of  surface  waves  at  the  boundary,  directly 
and  by  the  conversion  of  P  and  S  waves.  Lord  Rayleigh  (1885  and  1900) 
predicted  the  existence  of  a  free  surface  wave,  but  Lamb  (1904  and 
1916a,  b)  was  one  of  the  first  to  study  in  detail  the  propagation  of 
pulses  in  the  half-space. 

The  general  literature  has  been  well  surveyed  by 
Davies  (1953) s  Miklowitz  (I960),  McLean  and  Caless  (1961)  and  by 
Ewing,  Jardetzky  and  Press  (1957).  Miklowitz  lists  135  references, 
and  Ewing  et  al.  give  a  voluminous  bibliography  and  have  reduced  much 
of  the  work  to  textbook  status.  As  is  becoming  common,  the  latter  is 
abbreviated  ?EJPf  in  this  dissertation. 

Useful  surveys  in  abstract  form  have  been  given  by 
Bartunek  and  Igel  (1955)  and  by  Nestler  (n.d.).  These  surveys  list 
seme  of  the  limited  circulation  reports  which  are  fairly  generally 
available.  Attention  is  also  drawn  to  the  Vela  Uniform  Periodic 
Information  Digest  (Vesiac)  and  the  Preliminary  Bibliography  of 
Soviet  Seismology  and  Seismometry  (A.I.D.  1961).  Any  fresh  survey, 
so  soon  after  the  above,  would  merely  be  a  transcription,  and 
attention  is  therefore  limited  to  noting  some  papers  which  are  of 
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particular  interest  in  the  present  research. 

1.12  Following  upon  the  pioneer  work  of  Rayleigh  and  Lamb  (l.c.) 
and  of  Stoneley(l924)  on  interface  waves,  Sneddon  (1951)  made  a  theoretical 
study  by  Transform  methods  of  the  seismic  effect  of  a  travelling  source  of 
pressure  moving  over  the  boundary.  This  was  for  the  case  of  a  pulse  travel¬ 
ling  at  a  velocity  less  than  that  of  shear  waves  in  the  medium  (sub-sonic 
case)  and  the  treatment  was  extended  by  Cole  and  Huth  (195$)  to  the  tran¬ 
sonic  case  (velocity  between  that  of  shear  and  compressional  waves)  and 

the  supersonic  case  (velocity  greater  than  that  of  compressional  waves). 

The  nomenclature  is  not  free  from  ambiguity. 

In  these  cases  the  load  was  assumed  to  have  been  moving  for 
some  time  so  that  the  steady  state  solution  was  the  one  of  interest. 

In  such  models  the  displacement  tends  to  infinity  at  the  singular  point 
of  the  Rayleigh  velocity,  and  also  contains  arbitrary  constants  as  time 
tends  to  infinity.  These  inherent  ambiguities  may,  in  theory,  be  removed 
by  using  an  axially  symmetric  model,  with  peak  overpressure  and  the  velo¬ 
city  of  the  pulse  both  decaying  with  distance  from  the  origin.  This 
approach  has  been  undertaken  by  Miles  (i960),  but  it  is  necessary  to 
introduce  many  simplifying  assumptions  before  the  problem  becomes  amenable 
to  solution. 

1.13  Bateman  (193$)  discussed  the  interaction  between  ground 
and  air  resulting  in  free  vibrations  having  the  characteristics  of 
Rayleigh  waves  except  that  the  velocity  is  slightly  less  than  the  velo¬ 
city  of  sound  in  air,  rather  than  less  than  that  of  the  shear  wave  in 

the  ground.  Press  et  al.  (1951)  showed  experimentally  that,  under  certain 
circumstances,  resonant  coupling  may  exist,  producing  a  wave  superficially 
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similar  to  that  of  Bateman,  and  gave  an  explanation  based  on  an  extension 
of  the  Lamb  approach  for  a  travelling  line  source.  This  treatment  allowed 
the  main  features  to  be  predicted,  and,  in  the  few  years  subsequent  to 
this,  theoretical  and  experimental  approaches  were  made  by  Press  and 
Ewing  (1951a,  b,  c),  Dobrin,  Simon  and  Lawrence  (1951),  Haskell  (1951)  and 
Jardetzky  and  Press  (1952) .  These  showed  that  the  explanation  was  adequate, 
at  least  for  the  case  where  the  pressure  pulse  in  the  air  travelled  at 
essentially  constant  velocity. 

It  appears  to  be  generally  agreed  (see  e.g.  EJP,  page  230 
et  seq.)  that  the  necessary  and  sufficient  condition  for  the  production 
of  TPress-EwingT  type  coupled  waves  is  a  layered  medium,  in  which  dis¬ 
persive  Rayleigh  waves  can  exist  and  can  propagate  at  some  frequency 
whose  phase  velocity  is  equal  to  that  of  the  velocity  of  the  travelling 
pulse  of  pressure  in  the  air. 

Under  these  circumstances  a  train  of  constant  frequency  waves 
is  generated  at  the  surface,  with  Rayleigh-type  orbits,  starting  roughly 
in  coincidence  with  the  arrival  of  the  airblast  and  lasting  until  a  time 
consistent  with  an  arrival  at  the  group  velocity  corresponding  with  the 
selected  frequency. 

1.14  The  cited  papers  do  not  make  explicit  reference  to  the 
source  strength  or  dimensions,  and  it  would  appear  that  the  solutions 
offered  Imply  a  constant  frequency  train  independent  of  charge  size. 

The  period  is  fixed  by  the  particular  member  of  the  dispersed  train 
which  has  a  phase  velocity  equal  to  the  constant  pressure-pulse  velo¬ 
city.  Dr  Dobrin  appears  to  agree  with  this  interpretation  of  the  air- 
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coupled  wave  discussed  in  his  papers  .  Experimentally,  the  theory  has 
been  verified  using  small,  constant  size  charges,  monitored  at  points 
sufficiently  remote  to  allow  the  airblast  velocity  to  decay  to  a  con¬ 
stant  value,  as  envisaged  by  the  theory. 

Miles,  however,  in  his  theory  introduces  a  pulse  of  decaying 
velocity  and  also,  implicitly,  the  size  of  the  charge.  However,  he 
makes  no  reference  to  the  air-coupled  wave,  presumably  because  he  is 
dealing  with  a  semi-infinite  medium  and  assumes  a  non-dispersive 
Rayleigh  wave.  In  this  case  a  singularity  arises  only  at  the  single 
instant  in  time  when  the  blast  wave  velocity  is  equal  to  the  Rayleigh 
velocity. 

Dix  (1955)  suggests  that  long-period  waves  such  as  those 
seen  in  ground  roll,  can  be  generated  by  a  form  of  energy  retention 
at  the  source,  through  backward  propagation  of  waves  deriving  from 
the  interaction  of  the  spreading  wave  front  and  the  various  boun¬ 
daries.  Eventually  the  energy  is  dissipated  by  transmission  in  all 
the  available  modes,  including  the  Rayleigh  wave.  Dix,  however,  made 
no  mention  of  air-coupling  and  his  approach  is  not  really  applicable. 

Kogan,  Pasechnik  and  Sultanov  (1959)  discuss  the  differences 
in  periods  associated  with  explosions  and  earthquakes,  and  suggest  that 

*  In  private  communication* 

TIt  would  be  quite  essential  to  show  that  the  frequency  of  the 
air-coupled  wave  is  the  same  as  that  predicted  from  the  dispersion 
curve  from  a  hole  shot  at  the  velocity  of  sound  in  airT...’I  feel 
that  such  air-coupling  requires  a  dispersion  curve  which  crosses 
the  velocity  of  sound  in  air,  giving  a  predictable  frequency  for 
the  air-coupled  wave  corresponding  to  the  point  of  crossing.’ 
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the  short  periods  are  due  to  the  small  size  of  the  explosive  source. 

This  suggests  that  the  dominant  periods  will  change  with  charge  size, 
but  they  do  not  make  the  point  explicitly. 

1.15  The  Suffield  test  sites  were  chosen  primarily  as  being 
suitable  for  studying  the  propagation  of  blast  waves  from  large  charges 
to  low  overpressures  over  flat  uniform  ground.  However,  the  low  seismic 
surface  velocity  and  the  horizontal  layering  make  the  site  very  suitable 
for  study  of  surface  waves,  with  particular  attention  to  the  occur¬ 
rence  of  classical  or  ’Press-Ewing’  coupling  between  the  surface  wave  and 
the  airblast.  Early  work  showed  that  air-coupling  certainly  existed,  but 
it  appeared  that  there  were  departures  from  the  classical  form  and  an 
apparent  variation  in  the  period,  depending  upon  the  charge  size.  The 
seismic  participation  was  therefore  planned  essentially  in  terms  of 
these  air-coupled  waves,  though  without  prejudice  to  the  possibility  of 
fruitful  excursions  from  the  main  theme.  To  the  best  of  the  writer*  s 
knowLedge,  this  is  the  only  series  to  date  in  which  routine  measurements 
have  been  made  on  surface  waves  close  to  the  source,  with  charges  of 
increasing  size  all  on  virtually  identical  flat  ground.  It  would  appear 
that  such  a  series  is  essential  for  the  complete  verification  of  many 
theories  which  have  been  propounded  on  the  basis  of  limited  evidence  and 
limited  verification.  This  requirement  justifies  the  publication  of  the 
complete  range  of  seismograms  in  a  way  suitable  for  further  analysis  by 
other  workers  with  specialized  interests.  Some  seianograms  have  been 
published  separately,  but  for  the  sake  of  completeness,  reduced  copies  of 
all  the  main  series  and  many  subsidiary  seismograms  are  included  in  this 
dissertation  as  a  separate  section. 
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1.16  It  is  convenient  at  this  point  to  introduce  several  con¬ 
cepts  which  are  relevant  but  which  are  not  treated  extensively  in  the 
following  sections o  The  term  air-coupled  wave  should  not  be  confused 
with  the  tern  coupled  wave,  introduced  by  Leet  (e.g.,  1939,  1950,  I960). 
The  trajectory  of  the  latter  wave  is  a  diagonal  line  in  the  sense 
push-left-up  then  pull-right-down,  or  the  inverse.  It  is  claimed  to 
have  been  observed  in  the  records  from  dynamite  blasts  and  the  first 
nuclear  test.  The  advent  of  nuclear  energy  stimulated  interest  in  the 
strong-motion  region,  and  historically  interesting  papers  are  Taylor 
(1940),  Penney  and  Reines  (1945),  Coon,  Houghton  and  Nobles  (1945)  and 
Leet  (1945a,  b) .  The  coupled  wave,  christened  ’C’  and  described  by 
Leet,  may  be  treated  as  a  combination  of  a  P  and  S  wave.  There  is  some 
argument  about  the  reality  of  the  wave  as  a  distinct  entity,  and  the 
present  writer  considers  it  likely  that  there  have  been  at  least  some 
cases  of  wrong  identification.  However,  an  apparent  coupling  between  a 
transverse  mode  and  a  radial-vertical  mode  surface  wave  is  invariably 
present  in  the  records  from  the  Suf field  site,  and  was  also  obtained  by 
the  writer  at  the  Nevada  Test  Site.  A  preliminary  discussion  of  this 
effect  has  been  given  elsewhere  (Jones,  1962c). 

Leet  (1944  and  1962) ,  and  others,  also  discuss  a  surface  wave 
with  prograde*  orbital  motion,  but  otherwise  similar  to  Rayleigh  waves. 
This  Leet  terms  a  Hydrodynamic  wave,  H,  and  observed  experimentally  at 

jJj 

This  spelling  appears  to  be  common  usage  among  seianologists.  Neither 
Webster  nor  the  New  Oxford  list  this  word,  though  both  list  TprogredeT 
(to  move  forward).  Webster  indicates  that  ’progrede’  is  ’rare1.  Pro¬ 
grade  appears  to  be  an  adoptive  English  antonym  to  retrograde,  rather 
than  a  derivative  of  progredi.  The  listed  antonym  to  retrograde  is 
’ progressive’ . 
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the  first  nuclear  shot.  Such  prograde  -waves  are  distinct  features  of  the 
Suffield  records,  and  are  discussed  in  some  detail  herein,  and  have  also 
been  obtained  elsewhere.  They  are  usually  interpreted,  after  Sezawa  and 
Kanai  (1939),  as  equivalent  to  the  antisymmetric  vibration  of  a  plate. 

An  elegant  treatment  by  ray  theory  has  been  given  by  Tolstoy  and  Usdin 

(1953). 

Leet  (i960)  is,  incidentally,  guilty  of  a  common  error  in 
stating  (on  page  75)  that  the  air  pressure  wave  from  an  explosion  is 
first  supersonic  "but  within  a  distance  of  about  20  to  50  times  the 
diameter  of  the  solid  explosive,  has  levelled  off  at  the  speed  of  sound 
in  air".  This  is  quite  incorrect.  The  distances  correspond  to  about 
6.7  ft  and  16.6  ft  from  a  1  lb  charge,  ground  burst.  At  these  distances 
the  respective  Mach  numbers  of  the  shock  are  1.503  and  1.10$.  In  the 
case  of  the  100  ton  charge  (roughly  20  ft  diameter),  the  corresponding 
distances  are  401  ft  and  101$  ft.  Even  at  a  distance  of  59,900  ft  from 
this  charge  the  Mach  number  was  still  1.0009"T'. 

It  is  possible  to  use  records  such  as  those  included  herein 
to  determine  the  motion  of  the  originating  source,  according  to  the 
methods  of  Sato  (1955a,  b,  1956  and  1959) .  This  method  depends  upon 
the  dispersion  in  the  wave  train,  and  has  been  applied  by  Kisslinger  (i960) 
to  small  explosions  at  Saint  Louis.  However,  due  to  the  intervening  effect 
of  the  airblast,  the  method  is  only  applicable  to  the  regions  of  the  traces 
ahead  of  the  airblast.  A  modification  of  this  approach  is  discussed  herein. 

The  discussions  of  Pekeris  (1955a,  b,  1959  and  I960)  are  relevant 
to  the  prediction  of  wave  forms  from  explosions.  Again,  it  is  probable  that 
only  portions  of  the  seismograms  may  be  used  to  verify  his  theory,  due  to 

* 

Data  supplied  by  T.  K.  Groves,  SES. 


the  masking  effect  of  the  airblast  which  is  not  included  in  the  theory. 

1.20  The  Mature  of  the  Force  Functions 

1.21  The  seismic  effects  are  created  by  forces  deriving  from 
the  detonation  of  an  explosive  charge  in  close  contact  with  the  surface. 

The  distribution  of  these  forces  cannot  be  given  in  a  theoretically- 
derived  form,  though  many  aspects  may  be  described  by  empirical  laws,  and 
some  theoretical  symmetry  laws  are  applicable  within  limits.  Initially, 
the  forces  derive  directly  from  the  expansion  of  two  hemispheres,  one  the 
shock  hanisphere  expanding  in  the  air  and  the  other  that  expanding  into, 
and  disrupting,  the  material  upon  which  the  charge  lies.  It  is  clear  that 
this  lower  hanisphere  very  rapidly  attenuates  into  elastic  stress  waves  of 
the  type  commonly  discussed  in  the  theory  of  seismology.  At  points  on  the 
surface  remote  from  Ground  Zero  (GZ)  there  will  be  weak  disturbances  due 
to  the  reflected  and  refracted  P  and  S  waves,  together  with  surface  waves 
which  will  depend  upon  the  detailed  structure  of  the  upper  layers  in  the 
ground.  Normally,  these  effects  are  recorded  by  instruments  of  high  to 
very  high  gain,  remote  from  GZ. 

YJhen  we  pass  beyond  the  region  of  disruption  of  the  ground 
(the  crater  and  plastic  zone),  we  must  consider  the  force  function  deriv¬ 
ing  from  the  expansion  of  the  shock  wave  in  the  air.  In  the  case  of  a 
charge  actually  on  the  surface,  this  is  an  expanding  hemisphere  at  all 
relevant  times.  If  the  charge  is  above  the  surface,  the  primary  shock 
from  the  charge  combines  with  the  reflected  shock  from  the  ’ image  point* 
to  form  an  almost  vertical  Mach  Stem.  At  some  distance  from  GZ  we  may 
treat  this  expansion  as  that  of  a  vertical  cylinder  whose  axis  runs  through 
GZ.  At  the  surface  the  two  cases  are  distinguishable  only  in  terms  of  the 
strength  of  the  shock  front. 


•  *0 1 
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1.22  Immediately  beyond  the  true  crater  the  forces  deriving 
from  the  expanding  shock  front  will  be  such  that  the  material  is 
strained  beyond  the  (effective)  elastic  limit  and  shows  permanent  set. 

The  annular  area  at  the  surface  which  shows  permanent  set  is  loosely 
defined  as  the  Plastic  Zone,  the  definition  being  a  loose  one,  since  the 
zone  of  permanent  set  may  not  coincide  precisely  with  the  onset  of  elastic 
transmission  of  energy  by  seisnic  waves.  In  addition,  though  the  air 
blast  tends  to  create  a  downward  permanent  set,  this  is  balanced  in  part 
by  the  upward  forces  in  the  ground,  and  the  set  at  a  given  point  may  be 
up,  down,  or  zero,  even  close  to  the  crater.  The  precise  balance  could 
possibly  be  defined  for  an  ideal  medium,  but  such  speculation  is  fruitless 
in  the  real  case,  as  has  been  shown  by  Penney  and  Raines  (1945) ,  Jones  and 
Krohn  (1960a,  b)  and  Jones,  Krohn  and  Dewey  (1962). 

1.23  No  theoretical  law  exists  for  the  rate  of  decay  of  pressure 
with  distance  from  GZ,  for  the  complete  range.  It  is  common  practice  to  fit 
the  data  empirically  to  a  law  of  the  type 

Z  =  A/r  +  B/r2  +  C/r3  where  Z  is  the  peak  pressure  and  r  is  the 
distance  from  GZ,  A,  B  and  C  being  fitting  constants  to  which  we  cannot 
attach  physical  meaning.  However,  this  law  is  better  for  nuclear  shots  than 
HE  shots  and  a  form  which  fits  HE  data  better  is 

Z  =  A  +  B  +  C 

(r  +  A0)  (r  +  Bq)2  (r  +  G0)'5 

In  practice,  the  pressure  values  are  derived  from  the  Rankine-Hugoniot 

relationship  Z  =  M*  (Mf  +  2)  (_ )  where  MT  =  (Mach  Number  -  1) 

Y  +  1 

This  we  get  from  the  radius -"time  relationship 


..... 
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and  thus  M’  = 


t  +  b 


(t  +  b)  ln(t  +  b) 1 /2 


+  .+ 


where  aQ,  a2,  a3  . b  ... 

of  sound  ahead  of  the  shock* 


b 


are  fitting  constants  and  c  is  the  speed 


This  relationship  is  far  more  general  than 


the  first,  and  may  be  used  to  fit  data  from  all  types  of  explosions.  For 
different  size  charges,  the  law  of  geometrical  similarity,  originally  for- 


jljjt 

mulated  for  explosive  charges  by  Hopkinson  is  obeyed  very  closely  by  all 


explosives,  provided  the  air  mass  entrained  is  more  than  ten  times  the  ex¬ 
plosive  mass.  A  discussion  of  this  scaling  law,  following  the  well-known 
methods  of  dimensional  analysis  (cf  Hubbert,  1937)  has  been  given  by 
Jones  (1953) .  Taking  this  scale  law  into  effect  we  may  re-write  the  law 
of  pressure  decay,  for  the  simplest  case  as 


A  +  B  +  C 
r/w1'''3  (r/w1/3)3  (r/w'/3)3 


and  similarly  modify  r  and  the  fitting  constants  in  the  other  equations. 
In  most  experiments  several  sets  of  data  are  acquired,  and  the  constants 
vary  slightly  for  the  sets  of  differently  acquired  data.  The  data  them¬ 
selves  are  only  distinguishable  by  this  ’best  fit’  constants  variation. 
Detailed  discussion  of  the  curves  is  apparently  only  available  in  classi¬ 
fied  literature,  but  results  from  the  Suffield  explosions  are  available 
(Groves  1963a,  b,  c) . 

Fig.  1  is  typical  of  the  type  of  prediction  plot  used,  being 
that  for  a  200,000  lb  charge  based  on  data  from  a  AO, 000  lb  charge 


* 

** 


T.K.  Groves.  Private  Communication, 
circa  1914*  Precise  reference  not  located 
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similarly  detonated  on  the  Suf field  site.  Data  giving  a  wider  range 
are  shown  in  Table  1.  Table  2  gives  the  Hopkinson’ s  Law  scaling 
factors  for  a  variety  of  charges.  Fig.  2  is  a  photograph  of  a  large- 
scale  trial,  typical  of  the  cine  photos  from  which  travel  time  data 
are  obtained.  The  expanding  hemispherical  shock  is  clearly  visible. 
Fig.  3  shows  a  typical  crater  produced  by  a  large  surface-burst 
charge. 

1.24  The  decay  of  pressure  with  time  at  a  given  point 
after  the  passage  of  the  shock  front  must  also  be  treated  more  or 
less  empirically.  It  is  usually  expressed  in  the  form  of  a  (modified) 
Friedlander  Equation  (Friedlander  1946)  as 

P  =  Ps(l  -  t/r)e~at/V 

where  P  is  the  instantaneous  pressure,  Ps  is  the  peak  pressure  and  T3 
is  the  positive  duration,  t  being  measured  from  the  time  of  arrival 
of  the  shock.  This  form  fits  the  region  of  positive  pressure 
reasonably  well,  but  ignores  the  low  pressure  region,  the  second 
shock,  and  the  negative  phase  which  is  usually  present  in  practice. 
The  constant  a  is  derived  empirically  from  the  experimental  pressure¬ 
time  curves.  The  applicability  of  this  equation  to  wide  ranges  of 
shock  strengths- is  in  doubt,  but  the  simple  form  given  above  is 
usually  considered  in  the  open  literature,  and  is  probably  valid  in 
the  low  pressure  regions.  Note  that  Hopkinson’ s  Law  is  incorporated 
through  the  parameter  T,  which  itself  scales  according  to  this  law. 

Some  of  the  techniques  of  measurement  used  in  these  trials 
have  been  given  by  Harvey  (1957),  Mitalas  and  Harvey  (1958),  Groves 
(i960)  and  Dewey  (1961  and  1962) .  A  detailed  discussion  of  the  pro¬ 
pagation  of  blast  waves  may  be  found  in  Bethe  et  al.  (1951).  A 


4,000 


200,000  LB.  T.N.T.  SURFACE  BURST 


TYPICAL  SHOCK -BUBBLE  PHOTOGRAPH 


FIG.  3 

40,000  LB.  CRATER  24  HRS.  AFTER  DETONATION 
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simplified  description  of  cylindrical  and  spherical  blast  waves  has  been 
given  by  Friedman  (i960),  while  an  approach  to  the  problem  of  plane 
waves  in  granular  materials  has  been  made  by  Katz  and  Woeber  (1960a,  b) . 
Parkin  (195$)  has  discussed  the  dimensional  scaling  of  some  aspects 
of  ground  shock. 

1.30  Statement  of  the  Problem  and  the  Scope  of  the  Dissertation 

Many  of  the  mathematical  approaches  cited  above  are  relevant 
to  the  interpretation  of  the  strong  motion  records.  In  particular 
the  formulation  by  Miles  of  the  Blast  Wave  problem,  Pekeris’s  formu¬ 
lation  of  the  surface  seismic  pulse  and  the  EJP  discussion  of  Rayleigh 
and  air-coupled  waves  in  layered  media  and  the  papers  by  Sezawa  and 
Kanai  and  by  Tolstoy  and  Usdin  have  direct  bearing.  Even  when  we  admit 
the  gross  simplifications  inherent  in  all  these  theories,  the  problem 
is  complex  and  not  always  amenable  to  explicit  solution.  The  complete 
problem,  allowing  for  the  non-analytic  form  of  the  blast  propagation, 
plastic  deformation  in  an  area  comparable  in  diameter  to  the  shot-detector 
distance,  and  the  non- ideal  medium,  is  mathematically  over- complex . 
Nevertheless,  the  partial  theories  are  approximately  correct  and  allow 
satisfactory  prediction  of  limited  aspects  of  real  seismograms.  Generally 
speaking,  in  this  dissertation  the  interpretation  is  based  upon  these 
theories. 

A  comprehensive  body  of  data  from  the  Suffield  sites  is 
included  in  such  a  way  that  additional  interpretation  may  be  undertaken 
by  other  workers.  In  addition,  some  data  obtained  in  subsidiary  experi¬ 
ments  both  on  and  off  the  main  test  site  are  included  and  discussed. 


■ 
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Theoretical  dispersion  curves  for  a  variety  of  models , computed 
for  the  writer  at  the  Computation  Centre,  University  of  Toronto,  using 
the  programme  provided  to  the  writer  by  Don  L.  Anderson  ,  are  included 
and  compared  with  the  experimental  data.  The  models  chosen  were 
selected  by  the  writer  to  cover  both  the  best  estimate  of  the  actual 
ground  conditions  and  a  number  of  comparatively  standard  systems  likely 
to  be  of  more  general  interest. 

It  is  shown  that  the  data  produced  at  SES  are  internally 
consistent,  and  demonstrate  a  number  of  effects  which  have  not  been 
detected  in  a  definitive  way  previously. 

The  modified  air-coupled  wave  has  not  been  reported 
previously,  so  far  as  the  writer  is  aware,  and  no  explicit  theoretical 
treatment  exists.  However,  the  Press-Ewing  theory  implies  that,  on 
changing  the  assumption  of  constant  velocity  in  the  air  pulse  to  a 
decaying  velocity,  the  equivalent  theory  would  indicate  that  a 
succession  of  frequencies  would  be  excited  as  a  modified  fom  of 
air-coupled  train.  Two  variants  of  the  classical  theory  of  air¬ 
coupling  are  presented  in  Appendix  A. 

For  the  data  to  be  of  value  in  further  interpretation,  it 
is  essential  to  define  the  experimental  conditions  in  some  detail, 
and  for  this  reason  both  the  test  site  and  the  field  operations  are 
described.  An  appendix  gives  the  main  conclusions  obtained  from 
refraction  and  velocity  log  surveys  carried  out  on  the  test  site  by 
the  Century  Geophysical  Corporation  of  Canada  in  association  with  the 
writ  er . 

* 


Seismology  Lab.,  Pasadena.  See  Anderson  (1962)  and  Harkrider  and 
Anderson  (1962) . 
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Some  interpretation  of  the  transverse  component  of  motion  has 
been  attempted,  and  is  given  in  a  separate  section.  It  is  shown  that 
the  surface  wave  transverse  motion  is  repeatable  in  both  amplitude  and 
detailed  profile  in  a  given  azimuth  from  repeated  shots  on  the  same 
ground  zero,  even  though  the  crater  has  been  excavated  and  refilled 
between  shots.  A  note  on  this  effect  has  been  published  (Jones,  1962b) 
and  it  is  understood  that  a  similar  effect  has  been  noted  in  charges 
detonated  in  water-filled  craters  on  the  Florissant  Site  by  the  Saint 
Louis  Group. 

The  correlation  between  the  Love  and  Rayleigh  type  events  is 
discussed,  and  has  also  been  discussed  elsewhere  by  the  writer  (Jones, 
1962  c). 

It  is  shown  that  a  ’synthetic’  seismogram,  showing  the 
characteristics  of  a  multi-ton  surface-burst  charge  seismogram  may 
be  produced  by  combining  two  seismograms  from  a  pair  of  similar 
1000  lb  charges,  one  detonated  on  the  surface  and  one  detonated  fully 
contained.  A  brief  discussion  of  this  problem  has  been  given  by 
Jones,  Maureau  and  Cyganik  (1963) . 

Several  distinct  phases  of  Rayleigh  type  may  be  detected  in 
the  seismograms  from  sub-surface  charges.  The  phases  alternate 
between  the  pro grade  and  retrograde  habit,  the  two  dominant  phases 
being  a  prograde  phase  preceding  a  retrograde  phase.  No  reference 
has  been  located  in  the  literature  to  this  type  of  alternation,  though 
both  prograde  and  retrograde  modes  are  known  experimentally  and  pre¬ 
dicted  by  theory. 

Many  aspects  of  the  ground  shock  problem,  dealing  with 
cratering  and  plastic  deformation  as  distinct  from  purely  seismological 
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data,  have  been  published  separately  by  the  -writer  and  are  referred 
to  in  the  body  of  the  dissertation. 

Analysis  of  some  of  the  seismograms  by  a  method  similar  to 
that  of  Sato  is  included,  and  this  analysis  reveals  several  features  of 
interest  which  have  a  bearing  on  the  interpretation  of  the  multi-ton 
seismograms  and  on  the  identification  of  explosion  sources. 
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2.00  EXPERIMENTAL  WORK 

2.10  Equipment  and  Method 

2.11  The  experimental  methods  used  in  the  present  work  have 
been  described  in  previously  published  reports  (Jones,  I960;  Jones  and 
Winfield,  I960;  Jones,  Kisslinger  and  Cyganik,  1961;  and  Vesso,  Sanders 
and  Winfield,  1963).  The  description  herein  is  a  general  one  of  the 
developed  system,  and  it  should  be  realized  that  some  of  the  early  experi¬ 
ments  used  simpler  systems. 

2.12  The  main  instrumentation  consisted  of  from  one  to  eleven 
Sprengnether  three-component  portable  seismographs.  These  instruments 
have  a  natural  period  of  0.75  sec.  and  are  damped  to  0.55  critical.  There 
are  two  magnetically  damped  inverted  pendula  for  the  radial  and  transverse 
components  and  a  magnetically  damped  spring  supported  beam  for  the  vertical 
component.  The  instruments  had  no  provision  for  recording  zero-time  or 
the  time  of  arrival  of  the  airblast,  but  modifications  were  introduced  to 
provide  these  additional  data.  Details  have  been  given  by  Vesso,  Sanders 
and  Winfield  (1963).  Basically,  they  consist  of  installing  two  neon  lamps 
in  each  magazine,  the  light  passing  through  snail  holes  to  record  on  the 
photographic  paper  at  the  correct  position.  There  are  two  independent 
circuits «  the  zero  time  circuit  which  links  all  instruments  and  applies 

a  triggering  voltage  to  discharge  a  capacitor  through  the  neon;  secondly, 
individual  circuits  associated  with  each  instrument,  in  which  a  microphone- 
operated  relay  triggers  a  similar  capacitor  system.  These  time  marks  are 
impressed  with  an  accuracy  of  a  few  milliseconds,  certainly  better  than  the 
precision  of  reading  events  from  the  paper  record. 
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2.13  The  paper  speed  is  not  governed  and  varies  slightly 
with  the  ambient  conditions,  the  amount  of  paper  in  the  magazines  and 
the  state  of  the  batteries.  Timing  marks  are  impressed  every  0.02  sec., 
with  every  0.2  mark  being  thickened.  These  time  marks  appear  to  have 

a  cumulative  error  of  less  than  0.01  sec.  in  ten  seconds.  A  large 
number  of  consecutive  tests  on  all  instruments  gave  a  true  time  of 
0.2  +  .0005  sec.  for  the  longer  time  interval  .  This  appeared  to 
be  independent  of  the  state  of  the  batteries  provided  they  were 
reasonably  well-charged,  as  was  the  case  in  all  the  trials,  with  a 
single  known  exception  (seismogram  54 W) .  The  effect  of  large  changes 
in  ambient  temperature  has  not  been  investigated  since  the  majority  of 
the  trials  were  held  in  summer. 

2.14  The  majority  of  the  instruments  had  static  magnifi¬ 
cations  of  about  fifty,  though  some  with  magnifications  of  7.5  and 
of  about  150  were  used  in  a  few  of  the  later  trials.  Details  are 
given  in  Table  U,  together  with  the  other  acquisition  data.  The 
instruments  appear  to  satisfy  the  design  criteria  for  portable 
seismographs  given  by  Duval  (i960) ,  and  have  the  standard  response 
curve  shown  in  Fig.  4.  Mechanical  seismographs  of  this  type  act  as 
displacement  meters  for  periods  less  than  about  half  the  natural 
period,  and  as  accelerometers  for  periods  greater  than  about  twice 
the  natural  period.  The  dominant  periods  recorded  in  the  present 
work  were  0.16  to  0.34  for  500  lb  charges,  0.35  to  0.60  for  10,000  lb 
charges,  0.40  to  O.65  for  40,000  lb  charges  and  0.40  to  0.90  for  the 
200,000  lb  charge.  This  bracketing  is  unfortunate,  but  the  effect  is 
not  serious  as  may  be  seen  from  Fig.  5>  which  is  a  nomogram  after 

*  Detemined  by  B.  Sanders  for  the  writer. 
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Richter  (195$).  It  will  be  seen  that  the  dynamic  magnification  is 
virtually  constant,  and  the  phase  shift  is  about  ±  10  deg.  centred 
around  60  degrees.  For  any  given  charge  size,  the  periods  are 
virtually  constant,  and  no  great  error  is  introduced  by  using  the 
direct  seismograph  readings  in  deriving  particle  trajectories. 

Care  should  be  taken  when  we  attempt  to  compare  the  phase  relation¬ 
ship  between  waves  of  different  periods.  If  we  disregard  the  records 
from  the  highest  pressure  levels,  the  records  are  not  of  the  impulsive 
start  type,  and  so  the  dynamic  magnification  curve  for  non-impulsive 
conditions  is  acceptable. 

2.2 0  Description  of  Test  Sites 

Two  test  sites  were  used  for  the  main  series  of  experiments, 
the  sites  lying  alongside  one  another,  and  known  as  Drowning  Ford  and 
Watching  Hill  respectively.  Little  detailed  geological  information  is 
available  about  the  area,  though  there  are  two  Geological  Survey  Memoirs 
covering  the  Southern  Plains  area  (Williams  and  Dyer,  1930J  and  Russell 
and  Landes,  19AO) .  The  Pleistocene  geology  is  also  described  briefly 
by  Common  (1962) . 

Towards  the  end  of  the  present  series  of  experiments,  a 
velocity  and  refraction  survey  was  prepared  by  the  Century  Geophysical 
Corporation  of  Canada,  in  association  with  the  writer.  The  relevant 
details  of  this  operation  are  discussed  in  Appendix  B.  The  following 
description  is  based  upon  the  writer’s  field  observations. 

The  sites  are  towards  the  southern  end  of  the  Ross  Depres¬ 
sion  (see  Fig.  6)  which  appears  to  be  an  old  drainage  channel  draining 
the  large  morainal  area  in  the  vicinity.  It  appears  probable  that  the 


' 


TOPOGRAPHICAL  MAP  OF  TEST  SITES 


-  19  - 


deposits  to  a  depth  of  some  200  ft  (on  the  Watching  Hill  site)  are 
lacustrine,  being  formed  at  a  time  when  the  Ross  Depression  and  an 
area  to  the  South  and  West  contained  a  rather  large  lake.  What 
appears  to  be  a  late  boundary  of  this  lake  is  visible  running  from 
the  river  circling  North  to  North-West,  around  the  Ross  Depression* 

This  boundary  has  the  form  of  a  distinct  ’ shoreline®  of  heavy  boulder 
drift,  the  boundary  showing  a  sharp  transition  between  the  almost 
stone-free  lacustrine  deposit  and  the  heavy  drift  beyond.  Juniper 
Flats  and  Dishpan  Lake  appear  to  be  non-perennial  vestiges  of  this 
lake. 

The  Drowning  Ford  site  has  been  cut  down  to  its  present  level 
by  retreating  meanders  of  the  present  South  Saskatchewan  river.  Stalker 
(1961)  indicates  that  the  area  is  close  to  the  confluence  of  the  pre- 
Glacial  Milk  River  and  the  pre-Glacial  Oldman  River.  However,  on  p.  11 
of  the  cited  reference,  it  is  stated  that  this  particular  branch  of 
the  buried  valley  has  not  been  investigated. 

The  deposits,  being  typically  lacustrine,  consist  mainly  of 
uniform  beds  of  clay  and  silt,  with  occasional  sand  lenses.  These 
lenses  appear  to  be  of  small  lateral  extent  except  for  one  at  a  depth 
of  about  50  ft  on  the  Watching  Hill  site.  This  particular  bed  is  a 
thin  stratum  of  sand-gravel  and  is  waterfilled,  though  apparently  it 
does  not  draw  from  any  appreciable  body  of  water  as  it  may  be  drained 
easily.  Comments  on  this  point  by  Century  Geophysical  Corporation  are 
given  in  Appendix  B,  together  with  some  drillers®  logs.  There  is  a 
variation  in  detail,  but  a  typical  soil  classification  down  to  the 
waterbearing  gravel  is  given  in  Table  3.  This  refers  to  an  exploratory 
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log  taken  when  exploring  the  ground  to  select  the  GZ  of  the  200,000  lb 
trial .  The  classification  is  based  upon  samples  provided  by  the 
writer  to  the  Waterways  Experiment  Station,  Vicksburg.  Details  of  the 
velocity  layering  are  not  given  here,  since  they  are  included  and 
discussed  in  para  3*20  and  Appendix  B. 

Since  some  periods  of  nearly  one  second  have  been  noted,  at 
velocities  around  1000  ft/sec.,  it  would  appear  that  depths  somewhat 
greater  than  the  Pleistocene  deposits  are  of  significance,  though  for 
the  majority  of  the  records,  the  Pleistocene  geology  will  be  dominant. 

The  Pleistocene  deposits  are  underlain  by  Upper  Cretaceous  beds  of  the 
Foremost  formation,  which,  due  to  their  horizontal  attitude  (dip  in  a  few 
feet  per  mile)  are  effectively  ’conformal’,  as  far  as  the  layering  is 
concerned,  with  the  Pleistocene  deposits.  These  beds  are  mainly 
arenaceous,  with  soft  ’Pale  Beds’  alternating  with  non- continuous 
indurated  beds.  Some  shale  and  coal  seams  are  typical  of  the  system, 
but  these  are  of  small  thickness  and  the  dominant  seismic  members  are 
probably  the  ’Pale  Beds’  since  these  are  comparatively  massive  and 
continuous.  At  a  depth  which  is  greater  than  that  significant  for 
surface  waves  of  the  type  under  study,  but  effective  in  refracting 
P  and  S  waves,  the  area  is  underlain  by  Mississippian  limestone  with 
a  velocity  of  about  20,000  ft/sec. 

2.30  The  Primary  Experiments 

All  the  seismograms  from  the  primary  series  of  experiments 
are  included  in  Section  5»00.  Information  on  the  coding,  experimental 
conditions  and  so  on  are  given  in  Table  4. 
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2.31  $  lb  Charge  Series.  This  series  was  carried  out  on  the 
Drowning  Ford  site,  using  a  single  instrument  at  various  distances  from 
the  charges.  Each  charge  was  on  undisturbed  prairie,  and  cratering 
effects  have  been  reported  by  Jones,  Spackman  and  Winfield  (1959).  The 
zero  timer  was  under  development  so  few  zeros  are  recorded.  There  is  no 
indication  of  the  arrival  of  the  airblast.  The  instrument  was  not  really 
suitable  for  these  charges,  due  to  the  small  gain  and  low  time  resolution. 
Nevertheless,  the  seismograms  have  some  value  in  association  with  other 
charges  on  this  site  and  elsewhere.  The  three  components  of  Seismogram  4D, 
at  75  feet  from  the  charge,  are  included.  The  remaining  seismograms  at 
this  scale  are  only  shown  as  single,  radial  components.  (At  this  charge 
size  only,  the  vertical  and  transverse  components  are  less  consistent 

than  the  radial  in  profile) . 

2.32  30  lb  Charge  Series.  These  charges  were  better  suited 
to  the  gain  and  paper  speed.  Zero  marks  are  available  on  all  records, 
but  the  time  resolution  is  only  about  0.01  sec.  No  time  of  arrival  of 
the  airblast  is  recorded.  The  charges  were  cast  hemispheres  initiated 
at  the  centre  of  the  flat  base,  which  was  in  contact  with  the  prairie. 
Cratering  data  have  been  given  by  Jones  et  al.  (1959).  The  seismograms 
are  shown  in  sets  representing  one  component  of  motion  at  a  variety  of 
distances. 

2.33  60  lb  Charge  Series.  These  charges  were  centrally 
initiated  cast  spheres  resting  on  the  prairie.  Zero  time  is  available 
on  all  records,  but  not  the  air  blast  arrival  time.  Cratering  data 
have  been  given  by  Jones  et  al.  (1959).  The  seismograms  are  similar 
in  form  and  distance  to  the  30  lb  records,  and  are  given  in  sets 
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showing  one  component  at  a  time.  The  vertical  component  at  all  dis¬ 
tances  shows  what  is  postulated  to  be  a  refracted  P  wave  as  a  first 
arrival,  a  single  crest  showing  at  the  shortest  distance  and  a  train 
of  crests  and  troughs  at  the  longest  distance.  At  this  distance, 

$00  feet,  the  train  of  refracted  P  waves  has  not  separated  completely 
from  the  strong  motion.  The  phase  also  shows  at  much  reduced  ampli¬ 
tude  in  the  radial  component.  Similar  records  have  been  obtained 
from  identical  charges  detonated  with  their  diametral  planes  at 
ground  level.  The  records  are  basically  similar  to  those  from  the 
surface  bursts  and  are  not  discussed  further  in  this  dissertation. 

2.34  500  lb  Charge  Series.  These  charges  were  tran¬ 

sitional  in  type  between  the  cast  charges  used  for  the  smaller  trials 
and  the  block-built  charges  used  for  the  larger  trials.  This 
transition  was,  at  the  time,  the  main  purpose  of  the  series.  Charges 
30D  and  31D  were  built  from  small  blocks  of  TNT,  charge  32D  was  a 
charge  cast  in  ’slices’  in  a  ring  mould  and  charge  330  was  a  solid 
cast  hemisphere.  Details  of  the  cratering  have  been  given  by  Jones 
et  al.  (1959),  and  details  of  the  methods  of  fabrication  have  been 
given  by  Pennie,  Philips  and  Holdsworth  (i960)  and  by  Philips,  Ditto 
and  Holdsworth  (i960) .  The  experiments  showed  that  the  block  method 
of  construction  was  entirely  adequate,  and  produced  results  indis¬ 
tinguishable  from  those  of  solid  cast  charges.  Two  distances  were 
used  for  the  seismic  work,  500  ft  and  about  1000  ft.  No  zero  marker 
was  available.  The  records  at  1000  ft  are  so  similar  that  they  are 
indistinguishable,  thus  verifying  not  only  the  uniformity  of  the 
charges  and  ground,  but  also  that  of  the  seismic  recording  system. 
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The  seismogram  for  500  ft  (30D)  and  one  of  the  1000  ft  seismograms 
(32D)  are  included  in  Section  3.00.  The  P  wave  shows  up  very  clearly 
in  the  vertical  and  less  distinctly  in  the  radial  trace.  On  the  ori¬ 
ginal  record  there  is  a  faint  displacement  in  the  transverse  trace 
also.  At  1000  ft  the  P  train  has  not  quite  separated  from  the  strong 
motion. 

In  1961  an  additional  500  lb  charge  was  detonated  on  the 
Watching  Hill  site,  as  an  extension  of  the  programme  of  additional 
experiments  suggested  by  C.  Kisslinger.  The  records  are  discussed 
later,  in  association  with  the  buried  charges  of  the  same  programme. 

2.35  10,000  Charge  Series.  Many  records  are  available  at 

this  charge  size,  since  in  addition  to  the  direct  series  detonated 
for  the  SES  programme,  the  Station  acted  as  agent  for  Bell  Telephone 
Laboratories  in  a  further  series.  Advantage  was  taken  of  all  these 
charges  for  seismic  recording.  The  charges  were  divided  fairly  evenly 
between  the  Drowning  Ford  site  and  the  Watching  Hill  site,  and  thus 
afford  a  means  of  comparing  the  two  ranges.  This  is  essential  if  we 
are  to  consider  results  from  the  two  sites  as  a  single  series. 

All  the  charges  were  block-built  charges,  initiated  at 
the  centre  of  the  base.  Cratering  data  for  many  of  these  charges 
have  been  given  by  Jones  et  al.  (1959)  and  Jones  and  Krohn  (1960a) . 
Some  of  the  seismograms,  and  preliminary  discussion,  are  contained 
in  a  preliminary  report  (Jones,  I960). 

The  ten  seismograms  34D-43D  were  the  earliest  ones  obtained 
on  the  Drowning  Ford  site.  All  the  records  are  closely  similar  in 
form  with  the  exception  of  the  closest,  that  at  2000  ft.  This  record 
shows  a  double  pulse  of  large  amplitude  in  the  vertical  and  radial 
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traces  at  the  start  of  the  strong  motion.  This  may  be  compared  with 
a  pulse  of  similar  profile  obtained  at  the  Nevada  Test  Site  (see  para 
2.40).  Many  of  the  records  have  both  zero  time  and  the  airblast 
arrival  recorded.  Although  a  fresh  GZ  was  chosen  for  each  shot,  the 
seismograms  may  be  treated  as  a  single  series  since  all  the  measure¬ 
ments  were  taken  along  the  same  seismic  line.  In  one  shot,  a  rotation 
of  two  or  three  degrees  was  required  in  the  instrument  alignment. 
Seismograms  44W-46W  were  obtained  on  the  Watching  Hill  site,  in  a 
similar  way  to  the  Drowning  Ford  records.  Seismograms  47W  to  5&W 
were  also  obtained  on  the  Watching  Hill  site  but  for  this  series  re¬ 
consolidated  GZfs  were  used,  and  the  instruments  were  disposed  in 
azimuth  around  GZ  as  a  first  attempt  to  study  the  symmetry  of  the 
wave  propagation. 

2.36  405000  lb  Charge  Series.  The  first  40,000  lb  charge 

was  detonated  on  the  Watching  Hill  site  in  August,  I960,  as  the  focus 
of  a  large  international  trial  involving  agencies  from  Canada,  the 
United  States  and  Britain.  Details  of  the  experiments  have  been 
published  (SES  I960,  1962) .  Data  relating  to  cratering  effects  and 
permanent  and  transient  displacements  are  contained  in  Jones  and 
Krohn  (1960b) .  Data  on  the  power  spectrum  in  the  plastic  zone  have 
been  given  by  Barton  (i960) . 

Five  seismographs  were  used  on  a  single  seismic  line, 
from  3000  ft  to  10,000  ft.  Copies  are  included  as  59W-63W  herein. 

The  seismic  line  coincided  approximately  with  that  of  44W,  45W  and 
46W  at  the  10,000  lb  level.  In  addition  to  these  main  series 
seismograms,  records  were  obtained  from  a  25  cps  geophone  spread 
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at  55*000  ft  on  an  extension  of  the  same  line.  The  latter  records 
are  discussed  in  para  2. 43. 

The  second  40,000  lb  charge  was  detonated  on  the  Drowning 
Ford  site  on  17  August,  1962.  The  charge  was  used  primarily  for 
test  purposes  by  the  Bell  Telephone  Laboratories,  but  some  partici¬ 
pation  of  a  pure  research  nature  was  included  by  SES  teams.  The 
seismic  participation  consisted  of  a  single  line  of  instruments 
every  1000  ft  from  3000  ft  to  7000  ft  from  GZ.  Excellent  records 
were  obtained  which  included  zero  time  marks  and  time-of-arrival 
marks  on  all  records.  These  records  are  included  herein  as  75D  to 
79D  and  provide  a  direct  comparison  between  the  Watching  Hill  site 
(records  59W-63W)  and  the  Drowning  Ford  site  at  identical  distances 
from  two  similar  charges. 

In  addition  to  the  on-site  participation,  seismic  records 
were  obtained  by  teams  from  the  University  of  Alberta  (U.  of  A. — SES 
MOHO  project*)  in  the  vicinity  of  High  River.  Some  aspects  of  the 
records  obtained  are  discussed  elsewhere  in  this  dissertation.** 

2.37  200,000  lb  Charge  Experiment 

This  charge  was  detonated  at  17s 30*1,  GMT  on  3  August,  196l, 
as  the  focus  of  the  largest  international  trial  yet  held  at  Suf field. 
It  is  believed  that  up  to  the  present  this  is  the  largest  non-nuclear 

A  programme  undertaken  by  the  University  of  Alberta  to  investigate 
the  true  attitude  of  the  MOHO  and  Conrad  discontinuities  under  the 
western  plains.  The  University  was  supported  by  grants  from  DRB 
and  from  A.R.P.A.,  and  by  Munitions  and  Field  teams  from  SES. 


A  third  40,000  lb  charge  was  fired  July,  1963.  See  127D-131D 
Table  4. 
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detonation  of  this  controlled  type  in  history*.  Field  data,  and 
a  listing  of  projects  have  been  published  (SES  1962) .  Obviously 
it  will  be  several  years  before  all  the  data  have  been  published. 
Papers**  available  to  date  include  Jones,  Krohn  and  Dewey  (1962), 
Jones,  Kisslinger  and  Cyganik  (1961),  Jones,  Reiniger  and  Cyganik 
(1962),  Jones  (1961),  Kisslinger  (1961,1562)  >  Kempster  (1962), 

Winfield  (1962),  Johnson  (1962),  Beare  (1962),  Dewey  (1962)  and 
Groves  (1962) .  The  high  pressure  level  ground  shock  spectrum  was 
recorded  by  Space  Technology  Laboratories  and  a  preliminary  report 
is  available  (Halsey  and  Barton,  1961) .  It  is  doubtful  that  this 
spectrum  is  directly  relevant  to  the  present  analysis. 

Excellent  records  were  obtained  by  G.  Frantti  of  the 
body  waves  and  surface  waves  near  the  village  of  Hays,  100  Km  from 
the  test  site.  Instruments  used  were  of  the  Benioff  three-component 
type. 

The  strong  motion  seismic  recording  was  undertaken  jointly 
by  the  writer  and  Dr  Kisslinger,  Saint  Louis  University.  By  combining 
the  instruments  from  the  two  establishments  in  a  single  network, 
exceptionally  good  coverage  was  obtained.  Eleven  instruments  were 
used  on  three  lines  radiating  from  GZ  as  shown  in  Fig.  7.  The  main 
line  was  similar  to  that  used  on  smaller  trials,  and  included  all  the 


*  A  similar  charge  of  1,000,000  lb  TNT  is  scheduled  for  July,  1964  at  SES. 
Classified  papers  not  listed. 
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Suf field  instruments.  The  other  two  lines  were  intended  to  detect 
azimuthal  asymmetry,  the  main  interest  of  the  Saint  Louis  team.  The 
2000  ft  instrument  was  of  low  magnification,  7*5  static,  and  was 
installed  in  a  shallow  pit,  being  firmly  bonded  to  the  ground.  The 
pit  was  then  covered  with  a  heavy  steel  plate  and  the  edges  sealed 
to  exclude  at  least  part  of  the  airblast.  The  remaining  instruments 
were  mounted  in  the  standard  way  directly  on  lightly  scraped  prairie. 

No  zero  time  was  recorded  on  the  trial  (i.e.  directly  on 
the  seismic  records),  due  to  the  failure  of  a  surge  suppressing  diode 
after  the  dress  rehearsal.  However,  zero  can  be  obtained  with  ade- 
quate  precision  by  a  variety  of  indirect  methods  (see  Jones,  Kisslinger 
and  Cyganik,  1961) . 

In  spite  of  the  very  long  records  (in  excess  of  45  seconds) 
it  was  found  that  the  trace  positions  at  the  end  of  the  run  were 
identical  with  the  pre-shot  positions.  There  was  no  detectable 
instrumental  drift  and  obviously  no  instrument  was  subjected  to 
accelerations  large  enough  to  cause  bodily  displacement  relative  to 
the  ground.  It  is  believed  that  the  low  magnification  instrument  at 
2000  ft  was  in  as  high  a  pressure  region  as  any  previously  used  for 
recording  of  this  type. 

The  seismograms  obtained  are  given  in  reduced  form  in 
Section  5.00.  It  will  be  appreciated  that,  due  to  the  length  of 
these  records,  considerable  reduction  in  size  is  necessary  for  printing 
purposes.  Nevertheless,  these  prints  are  adequate  for  measurement 
purposes  provided  very  minor  events  are  not  under  study. 
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2.40  Subsidiary  Experiments 

2.41  Contained  and  Venting  Charges  at  SES 
In  order  to  elucidate  certain  aspects  of  the  multi-ton 
records,  and  with  the  auxiliary  purpose  of  providing  sources  for  the 
U.  of  A. -SES  MOHO  project,  several  charges  of  1000  lb  TNT  were 
detonated  in  geometries  different  from  that  used  for  the  main  series 
of  experiments.  These  are  enumerated  and  the  records  presented  without 
comment  in  this  section.  Discussion  of  the  records  is  contained  in 
Section  3*50. 

Charge  FE  549-4  was  a  1000  lb  cast  sphere  of  TNT  buried  in  a 
shaft  dug  approximately  on  the  main  seismic  line  of  the  200,000  lb 
trial,  in  the  region  between  Watching  Hill  and  Juniper  flats.  The 
charge  centre  was  at  a  depth  of  23.5  ft,  chosen  to  correspond  with  one 
of  the  scaled  depths  of  the  Stagecoach  shots  (see  Section  2.42),  and 
the  shaft  was  stemmed  with  the  excavated  material  to  the  original  surface 
level.  Seismographs  were  placed  along  a  single  line  every  1000  ft  from 
1000  to  5000  ft.  The  charge  produced  a  large  cavity  collapse  crater 
with  some  venting,  but  the  airblast  pressure  was  very  low,  equivalent 
to  that  from  less  than  half  a  pound  of  TNT  at  the  surface  zero*.  The 
seismograms  obtained  are  included  herein  as  94J-98J. 

Charges  FE  549-1  and  -3  were  similar  charges  consisting  of 
cylinders  of  cast  TNT,  lightly  cased  in  a  case  6  ft  long  by  1$  inches 
in  diameter,  weighing  1000  lb  (TNT)  and  detonated  at  the  bottom  of  fully 
tamped  wells  53  ft  deep  on  the  flank  of  Watching  Hill.  On  the  first 
charge  seismographs  were  used  at  every  2000  ft  from  2000  ft  to  10,000  ft. 
The  recording  at  2000  ft  was  not  successful,  but  the  remaining  four 

T.K.  Groves  -  unpublished  data. 
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records  are  included  herein  as  100H-103H.  For  the  second  trial  the 
seismographs  were  disposed  in  azimuth  at  a  distance  of  2000  ft  from 
the  surface  zero.  Seismograms  are  included  as  104H-10$H. 

Charges  FE  552-7  and  -£  were  similar  in  geometry  and 
location  to  the  last  two  trials,  with  the  exception  that  water  tamping 
was  employed  instead  of  soil  tamping.  On  the  first  trial  of  this  pair 
the  seismographs  were  disposed  in  azimuth  at  2000  ft,  as  for  FE  549-3* 
while  for  the  second  trial  only  two  seismographs,  similarly  disposed 
were  employed  for  confirmation  purposes.  The  records  are  included  as 
109H  to  116H. 

LT  2$5-l  consisted  of  a  charge  of  296  lb  detonated  at  a 
depth  of  3  lambda  (20  ft)  on  the  Drowning  Ford  test  site.  This 
charge  vented,  though  previous  experience  on  other  test  sites  had 
indicated  that  the  charge  was  at  containment  depth.  This  may  have  been 
partly  due  to  the  fact  that  the  wooden  cribbing  was  not  removed  from 
the  shaft  during  the  stemming  process.  Records  were  obtained  at  500  ft 
and  at  every  1000  ft  from  1000  to  5000.  Records  are  included  as  $9D  to 
93D. 

2.42  40 >000  lb  Charges  at  the  Nevada  Test  Site 

The  writer  was  fortunate  in  being  invited  to  lead  a  team 
from  Canada  to  participate  in  Operation  Stagecoach  at  the  Nevada  Test  Site 
during  March  I960.  The  operation  was  carried  out  in  Area  10  of  the  test 
site,  in  close  proximity  to  the  craters  of  the  Teapot  S  and  Buster-Jangle 
Nuclear  shots,  and  close  to  the  later  trial  of  1,000,000  lb  TNT  called 
Operation  Scooter.  The  shallow  geology  of  the  site  is  one  of  a  thick 
deposit  of  desert  alluvium,  consisting  mainly  of  lightly  cemented  sands 
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and  gravels.  The  upper  levels,  to  a  depth  of  between  50  and  100  ft 
derive  mainly  from  a  limestone  source,  while  the  deeper  levels  to 
at  least  500  ft  derive  mainly  from  Rhyolite  and  Tuff. 

The  three  shots  of  Stagecoach  were  identical,  with  the 
exception  of  the  depth  of  burial  of  the  charges.  In  each  case,  a 
cylindrical  shaft  was  sunk  to  the  required  depth,  a  spherical  cavity 
excavated,  and  the  charge  built  in  this  cavity  as  a  nominally  spherical 
charge  of  40,000  lb  TWT.  Initiation  was  by  a  single  electric  detonator 
in  a  spherical  booster  charge.  The  excavated  material  was  returned  to 
the  shaft  as  stemming,  and  compacted  approximately  to  the  original 
density.  The  three  charges  were  placed  at  depths  $0  ft,  17.1  ft  and 
34.2  ft,  corresponding  to  2.3  lambda,  0.5  lambda  and  1.0  lambda.  The 
three  charges  were  functioned  at  intervals  of  about  one  week  in  this 
order. 

For  each  charge  a  single  seismic  line,  using  three  instruments, 
was  established,  instruments  being  at  3000,  4000  and  5000  ft  in  each 
case,  and  the  azimuth  being  kept  constant.  A  trial  record  covering 
the  operation  has  been  issued  (Jones  and  Winfield,  I960) .  The  records 
obtained  are  included  herein  as  SOW  to  SSN,  and  are  discussed  elsewhere 
in  this  dissertation. 

2.43  Records  at  Longer  Distances  from  SES  Charges 

Many  records  have  been  obtained  at  distances  from  the  SES 
charges  greater  than  those  discussed  in  the  previous  sections. 

Generally,  these  are  not  relevant  to  the  present  dissertation,  but 
certain  records  are  of  interest  in  showing  the  existence  of  air¬ 
coupling  at  quite  long  ranges. 
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During  the  I960  trial  of  40,000  lb,  the  writer  had  the  use 
of  a  Texaco  seismic  truck  belonging  to  the  University  of  British 
Columbia.  This  equipment  was  used  at  a  nominal  distance  of  56,000  ft 
from  GZ  on  an  extension  of  the  main  seismic  line.  The  ground  between 
GZ  and  the  recording  site  was  mainly  flat,  but  some  rolling  country 
intervened.  The  site  consisted  of  a  shallow  saucer- shaped  depression 
below  a  feature  known  as  the  Hogback.  The  site  lies  in  a  generally 
sandy  area,  but  is  itself  a  silt  deposit  not  unlike  that  at  Watching 
Hill. 

Twenty-two  geophones  of  25  cps  type  were  available,  and  these 
were  used  in  a  spread  some  2000  ft  long,  with  the  recording  truck  in 
the  centre  of  the  spread.  No  pre-suppression  was  used,  and  the  filter 
circuits  were  inoperative,  so  that  the  low  frequency  cutoff  was 
dependent  upon  the  geophones.  The  gain  was  set  approximately  equal 
on  each  channel,  and  was  approximately  the  maximum  available,  since 
the  area  was  found  to  be  generally  quiet.  In  view  of  the  anticipated 
record  length,  the  paper  was  allowed  to  run  freely  from  the  camera  and 
was  processed  immediately  before  opening  the  (stiflingly  hot)  cab. 

A  few  prelirainary  shots,  with  1/2  lb  charges  air-burst  a  few 
hundred  yards  from  the  spread,  produced  typical  records  showing  re¬ 
fraction  arrivals  well  separated  from  distinct  air-coupled  waves  of 
the  type  recorded  by  Dobrin  et  al.  in  their  investigation  of  air-coupled 
waves  (loc  cit) . 

The  radio  count- down  from  the  control  bunker  was  monitored 
at  the  site,  and  a  zero-time  accurate  within  the  operator’s  reaction 
time  was  recorded  by  releasing  two  held  down  ’line  test’  switches. 
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one  in  each  bank  of  recorders.  Visual  observation  of  the  monitor 
panel  detected  the  passage  of  a  refracted  wave  which  was  later  found 
to  be  barely  above  the  noise  level.  An  approximate  velocity  of 
15,000  fps  was  calculated  for  this  arrival.  The  precision  is  low, 
and  beyond  assuming  that  this  probably  represents  a  refraction  from 
the  Mississippian  Limestone,  little  can  be  said  of  it. 

Approximately  45  seconds  after  zero,  the  first  geophone 
detected  an  arrival,  and  the  others  picked  it  up  in  succession.  As 
the  signal  transferred  from  the  inside  bank  to  the  outside  bank,  the 
airblast  arrived  at  the  recording  vehicle.  An  assistant  outside 
heard  the  detonation  distinctly,  and  the  writer  inside  was  somewhat 
startled  by  the  violence  of  the  arrival.  The  vehicle  acted  as  an 
efficient  seismoscope.  The  relevant  portion  of  the  record  is  in¬ 
cluded  herein  as  Fig.  $.  The  velocity  of  the  arrival  over  the 
spread  length  turns  out  to  be  2.S0S/2.52,  or  about  1120  fps,  a 
reasonable  value  for  the  sound  velocity.  The  typical  period  of 
the  early  part  of  the  train  was  0.03  sec.,  but  in  view  of  the 
high  frequency  geophones,  this  is  not  significant.  The  record 
length  for  this  phase  was  around  ten  seconds  before  the  signal 
dropped  to  noise  level.  This  indicated  an  air coupled  train  of  at 
least  11,000  ft  -  a  fifth  of  the  epicentral  distance. 

The  writer,  both  from  the  record  and  from  his  physiological 
sensations  at  the  time,  is  in  no  doubt  that  this  was  classical  Press- 
Ewing  coupling.  The  pressure  pulse  at  this  distance  would  be  no 
more  than  a  few  hundredths  of  a  psi,  and  the  front  would  be  of 
essentially  constant  speed.  At  the  time  this  was  the  most  distant 


i 

twt 


-  33  - 


positive  identification  of  air-coupling  from  a  surface  burst 
explosive  charge.  However,  records  which  may  show  air-coupling 
were  obtained  by  G.  Frannti  at  Hays,  100  kms  from  the  200,000  lb 
GZ.  Benioff,  Ewing  and  Press  (1951)  also  give  evidence  of  the 
arrival  of  an  air  wave  at  265  km  from  a  natural  earthquake,  which 
they  interpret  as  being  an  air-coupled  wave  produced  in  a  track  length 
of  55  km  over  low  velocity  sediments,  and  then  refracted  down  to  the 
recording  site  from  a  vertical  height  of  40  km.  An  interesting  illus¬ 
tration  of  the  duration  of  the  seismic  effects,  and  of  the  large  sur¬ 
face  area  undergoing  disturbance  simultaneously,  is  provided  by  the 
fact  that  the  Sprengnethers  at  10,000  ft  were  still  recording  quite 
strong  motion  while  these  effects  were  being  recorded  at  56,000  ft* 

The  records  shown  in  Fig.  $  indicate  two  distinct  bursts 
of  energy  of  greater  intensity  than  the  mean  value.  The  data  are 
not  sufficient  to  allow  positive  interpretation  of  this  effect,  but 
it  is  probable  that  the  second  burst  is  an  indication  of  atmospheric 
focussing  at  some  point  in  the  track. 

During  the  same  trial,  records  were  made  by  Imperial  Oil 
Ltd.,  at  Calgary*. 

As  mentioned  in  para.  2.37j  a  team  from  the  University  of 
Michigan  made  recordings  of  the  seismic  effects  of  the  200,000  lb 
charge,  using  Benioff  equipment  near  the  village  of  Hays,  100  km 

A.  Mair,  Imperial  Oil  Ltd.,  personal  communication*  "We  heard 
two  distinct  arrivals  and  recorded  three".  (See  also,  Mair  (i960). 
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from  the  test  site.  Some  interpretation  of  the  records  has  been 
given  by  Frantti*  (1962) .  An  air-coupled  ground  roll,  presumably 
of  classical  Press-Ewing  type  was  detected,  but  it  is  interesting 
to  note  that  Frantti  was  able  to  show  close  similarities  between 
these  waves  and  the  general  microseismic  background  in  the  plains 
areas.  Whether  or  not  this  is  an  indication  that  much  of  the 
microseismic  noise  in  these  regions  is  produced  by  coupling  to 
wind  is  outside  the  scope  of  this  dissertation. 

Partly  at  the  writer’ s  request,  and  partly  to  provide 
data  for  another  study  ,  equipment  used  by  the  U.  of  A.  MOHO 
project  was  allowed  to  run  in  various  locations  until  times  con¬ 
sonant  with  the  arrival  of  the  airblast  from  certain  surface 
bursts. 


Paper  presented  at  S.E.G.  meeting,  Calgary,  1962.  (See  also  foot¬ 
note  p.  90) . 

**  E.  Reinelt,  Doctoral  Dissertation,  U.  of  A.,  Edmonton,  Alberta. 
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3.00  GENERAL  ANALYSIS 

3.10  Zero  Time  and  the  Arrival  of  the  Air blast 

The  correlation  of  zero  time,  and  the  time  of  arrival 
of  the  airblast,  is  of  importance  in  the  interpretation  of  the 
seismograms.  In  some  cases  these  data  are  available  directly 
from  the  impressed  marks  on  the  records.  In  all  cases  other  than 
the  very  small  charge  records,  the  arrival  of  the  airblast  is  also 
indicated  by  a  distinct  TblipT  on  the  records  caused  by  the  action 
of  the  airblast  on  the  case.  If  the  travel  time  of  the  airblast 
is  known,  zero  may  thus  be  obtained  even  in  the  absence  of  a  zero 
marker.  Where  the  refracted  P  wave  is  detectable  as  a  first  arrival, 
zero  time  may  also  be  estimated  from  the  known  travel-time  curve. 

Table  5  gives  the  airblast  data  abstracted  from  the 
seianograms,  the  same  data  corrected  for  ambient  conditions,  and 
the  corrected  data  scaled  according  to  Hopkinson’ s  Law  to  a 
10,000  lb  charge.  Corrections  for  ambient  conditions  are  not 
wholly  satisfactory  with  the  type  of  data  available.  No  correction 
was  made  for  ambient  pressure,  since  the  changes  about  the  standard 
of  13.69  for  the  site  were  small.  If  data  from  other  sites  are  to 
be  correlated  with  these,  however,  conversion  to  sea  level  would 
be  appropriate.  Corrections  were  made  for  the  mean  wind  by  using 
local  measurements  at  two  metre  height,  the  measurements  being 
taken  close  to  GZ  and  as  near  zero  time  as  practicable.  Similarly, 
corrections  for  the  mean  temperature  at  two  metres  were  made  to 
convert  the  readings  to  a  standard  of  32°F. 

^  • 
psi. 
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After  preparing  Table  5>  the  writer  attempted  to  obtain, 
for  comparison,  a  ’standard’  travel  time  curve.  It  appeared  that, 
while  ample  data  exist  out  to  about  600  ft  from  a  10,000  lb  charge, 
few  data  are  available  beyond  this  distance.  The  data  given  herein 
are,  therefore,  the  best  available  experimental  data,  and  an  expansion 
of  previous  data  by  a  factor  of  nearly  fifteen.  It  will  be  seen  that, 
for  these  moderately  short  distances,  the  velocity  of  the  front  is  well 
in  excess  of  the  sonic  velocity,  and  errors  would  be  involved  in 
assuming  a  constant  velocity  front  even  at  the  out-stations. 

3.20  The  P  Wave  and  Postulated  S  Wave 

At  points  sufficiently  remote  from  GZ,  the  first  arrival 
will  be  a  refracted  longitudinal  body  wave,  or  P  wave.  This  state¬ 
ment  is  true  provided  the  seismic  velocity  increases  in  some  way 
with  depth,  and  provided  the  gain  of  the  instrument  is  sufficient  to 
detect  the  low  energy  arrival.  In  the  present  dissertation,  analysis 
of  the  P  wave  is  restricted  to  identification  and  the  abstraction  of 
incidental  information,  since  the  experimental  system  was  quite 
unsuitable  for  a  study  of  this  event. 

The  P  wave  may  be  identified  as  the  short  train  visible  in 
the  vertical  component  of  Seismogram  35D,  starting  about  0.4  sec. 
after  zero.  From  this  point  it  is  possible  to  trace  the  event 
forward  and  back  In  the  various  records.  Table  6  gives  the  travel 
time  data  abstracted  for  the  event.  No  P  wave  data  are  quoted 
for  the  small  seismograms  (before  12D)  since,  in  the  interpretation 
of  the  writer,  the  airblast  arrival  is  the  first  recorded  event. 
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This  could  be  disputed,  but  confirmation  may  be  obtained  by  com¬ 
paring  the  amplitude  of  the  P  wave  in  11-13D  with  the  amplitude 
of  the  first  arrival  in  4-7D.  The  arrival  times  quoted  for  the 
P  wave  in  the  200,000  lb  trial  are  based  upon  the  calculated  zero 
(see  Jones,  Kisslinger  and  Cyganik,  1961). 

Since  there  are  two  test  sites  involved,  the  P  data  are 
not  necessarily  reconcilable  between  the  sites,  and  are  therefore  plotted 
separately  in  Fig.  9.  Actually,  the  Watching  Hill  data  agree  well  with 
the  Drowning  Ford  data  for  the  high  velocity  layer.  Although  not  very 
suitable  for  a  refraction  analysis  they  may  be  interpreted  in  the 
standard  way  (see,  for  example,  Dobrin,  I960).  This  is  worth  doing 
in  view  of  the  normal  seismic  survey  conducted  by  Century  (see 
Appendix  B) . 

For  the  Drowning  Ford  Site  the  data  indicate  a  three  layer 
system  consisting  of  unconsolidated  sediments  with  a  velocity  of  less 
than  2000  fps  down  to  70  ft.  A  velocity  of  about  6900  fps  between 
70  ft  and  300  ft  represents  the  horizontal  Foremost  beds  of  soft 
shales  and  sandstones.  Below  this  level  the  velocity  rises  to  about 
$400  fps,  possibly  indicating  the  somewhat  more  competent  Basal  Fore¬ 
most.  This  interpretation  may  be  compared  with  the  Century  data  given 
in  Appendix  B,  which  give  2200  fps  to  75  ft,  6000  fps  down  to  225  ft 
and  7700  below  this  depth.  The  agreement  is  satisfactory  considering  that 
the  present  analysis  is  based  upon  a  large  variety  of  seismograms,  taken 
at  different  times  with  instruments  of  low  gain.  The  P  wave  amplitude  is 
barely  above  the  recording  threshold  of  the  instruments,  and  the  strong 
motion  region  has  signal  strengths  at  least  an  order  of  magnitude  greater. 
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An  event,  marked  A  in  36D,  is  occasionally  to  be  seen  as 
an  abrupt  arrival  in  the  radial  trace,  associated  with  a  synchronous 
arrival  in  the  transverse  trace.  The  amplitude  is  low,  about  that 
of  the  P  wave,  and  precise  measurements  are  not  possible.  The 
travel  time  plot  is  shown  in  Fig.  9,  but  the  scatter  tends  to  show 
that  some  erroneous  identifications  are  included.  The  mean  step  out 
velocity  between  2000  ft  and  6000  ft  is  about  2160. 

This  event  could  be  interpreted  as  merely  the  start  of 
the  surface  wave  discussed  in  the  next  section,  but  it  is  more 
likely  to  be  a  distinct  phase. 

The  event  may  be  a  refracted  S  wave,  as  in  unconsolidated 
sediments  there  is  no  reason  to  suppose  that  the  a/3  ratio  is  that 
of  a  "Poisson1'  solid.  Most  of  the  seismologists  who  have  looked  at 
these  records  favour  this  interpretation,  despite  the  fact  that  this 
requires  the  sub-surface  layers  to  be  unconsolidated. 

3.30  On  Dispersion  in  Strong  Motion  Seismograms 

3.31  On  Group  and  Phase  Velocity 

Despite  the  apparent  simplicity  of  the  relationship  between 
Group  and  Phase  velocities,  much  confusion  is  introduced  by  the  fact 
that  the  phase  velocity  is  frequently  identified  with  the  velocity  of 
the  physical  wavelet  or  ’crests1  in  the  recorded  train.  It  seems 
expedient,  therefore,  before  discussing  the  strong  motion  regions  of 
the  seismograms,  to  clarify  precisely  what  is  to  be  understood  in  the 
present  context.  The  essential  point  is  that  the  common  practice 
among  seismologists  of  equating  wavelet  velocity  and  phase  velocity 
is  acceptable  only  in  regions  remote  from  the  source. 
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Let  us  define  Phase  velocity  in  the  usual  way  as  the  velocity 
of  propagation  of  an  infinitely  long  constant  frequency  constant  ampli¬ 
tude  sine  or  cosine  wave  -  essentially  a  Fourier  component.  We  bear 
in  mind  that  only  in  rare  circumstances  is  this  wave  directly  observable. 
In  a  non-dispersive  medium  the  Fourier  components  will  have  the  same  con¬ 
stant  velocity  independent  of  frequency.  Hence  if  we  compound  several 
components,  with  different  amplitudes  and  in  a  given  initial  phase 
relationship,  the  compound  wave  will  be  an  infinitely  repeating  train  of 
amplitude  groups.  The  train  of  groups  will  propagate  without  change  of 
form,  in  the  common  velocity.  If,  by  some  means,  we  can  produce  a  single 
group  from  the  train,  that  is  a  packet  starting  and  ending  with  zero 
amplitude  and  in  isolation,  we  may  represent  this  group,  within  its 
space  time  interval  only,  by  the  same  set  of  Fourier  components  by 
introducing  the  fiction  that  the  group  repeats  indefinitely  outside 
its  interval.  Two  situations  may  arise  in  such  a  case,  depending 
upon  the  actual  form  of  the  isolated  packet.  In  some  cases  it  will 
be  possible  to  replace  the  group  by  a  finite  number  of  components, 
to  the  necessary  accuracy  in  any  practical  situation.  On  the  other 
hand,  in  the  case  common  to  explosion  waves,  we  may  have  an  abrupt 
change  in  profile  at  some  point  within  the  group.  In  such  a  case  a 
strictly  infinite  number  of  components  is  required  to  describe  the 
wave  packet,  and  even  in  practice  the  packet  cannot  be  described  by 
a  small  number  of  components.  The  difficulty  can  be  overcome  by 
replacing  the  Fourier  component  by  the  Heaviside  function  if  this  tech¬ 
nique  is  acceptable  for  other  reasons.  (For  a  reasonably  full 
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discussion  of  this  aspect,  see  Pekeris,  1959). 

Now  let  us  consider  the  case  in  which  the  medium  is 
dispersive  -  that  is  to  say  each  frequency  of  Fourier  component 
propagates  at  its  own,  distinct  velocity.  We  assume  in  all  that 
follows  that  the  attenuation  of  the  Fourier  components  with  dis¬ 
tance  is  either  zero  or  at  a  rate  common  and  constant  for  all 
frequencies.  We  combine  several  such  components  in  a  given  ampli¬ 
tude  and  phase  relationship  at  time  zero  (as  each  component  is 
assumed  of  infinite  length  we  need  not  specify  the  location  of  the 
space  origin) .  This  superposition  of  Fourier  components  is  our 
initial  condition,  and  is  the  situation  considered  in  a  similar  way 
by  Lamb  (1932,  p.  3^0).  The  subsequent  motion  is  in  general  made  up 
of  a  system  of  component  waves,  each  travelling  at  its  own  velocity, 
moving  in  either  direction,  independently  of  each  other.  The  result 
will  be  a  compound  wave  which  will  exhibit  gross  maxima  and  minima 
of  amplitude,  the  envelope  of  which  may  enclose  many  individual 
physical  wavelets.  The  shape  of  the  envelope  and  of  the  individual 
wavelets  will  continuously  alter,  so  that  if,  for  example,  we  fix 
our  attention  upon  a  given  amplitude  maximum  in  the  envelope,  and 
move  with  it,  we  shall  see  the  envelope  curve  lengthen  with  time, 
and  the  shape  of  the  individual  wavelets  within  the  envelope  change 
in  keeping  with  the  envelope.  In  the  case  of  surface  waves  on  water, 
for  example,  we  see  amplitude  groups  move  across  the  surface,  each 
amplitude  group  being  made  up  of  wavelets  which  form  at  the  rear  and 
die  away  at  the  front,  the  wavelets  travelling  through  the  group  with 
a  velocity  greater  than  the  gross  velocity  of  the  group.  The  wave  let 


velocity  is  the  velocity  with  which  these  wavelets  travel,  and  this 
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is  not  the  same  as  the  phase  velocity  of  any  of  the  individual  Fourier 
components  into  which  the  group  may  be  analysed.  In  many  circumstances, 
however,  common  in  long  range  seismology,  where  we  are  observing  the 
wave  packet  well  removed  from  the  source  region,  dispersion  has  already 
acted  to  such  an  extent  that  a  given  region  of  a  wave  packet  contains 
only  a  very  narrow  range  of  frequencies  in  the  Fourier  components. 

In  this  case,  little  error  is  introduced  by  equating  phase  and  wavelet 
velocity.  This  is  emphatically  not  the  case  in  the  strong  motion  region 
where  active  dispersion  is  changing  the  wave  profile. 

We  may,  after  Lamb,  consider  the  group  obtained  by  super¬ 
posing  two  equal  amplitude  Fourier  components.  Each  component  will, 
of  course,  have  a  definite  wavelength  independent  of  whether  the 
wavelength  is  measured  as  the  wave  passes  a  fixed  point  at  the  phase 
velocity,  or  whether  we  look  at  the  actual  wave  profile  in  space. 

Summing  the  two  components  we  get* 

?  =  a  sin(kx  -  wt)  +  a  sin(kTx  -  w’t) 

=  2a  cos  ^l/2(k  -  kT)x  -  l/2(w  -  wT)tj  sin  [l/2(k  +  kT)x  - 

l/2(w  +  wT)t-j 

where  k  =  —  =  ~  ~  for  the  single  Mk”  component, 

c  cT  \ 

and  k’  =  ™  =  i~n  for  the  single  kf  component. 

c?  cTTT  \T 

Now  if  kT  is  very  nearly  the  same  as  k,  the  cosine  term  varies  but 
slowly  with  x.  The  surface  therefore  has  a  sine  wave  profile  in 
which  the  amplitude  changes  periodically  between  the  values  0  and  2a. 
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The  distance  between  the  centres  of  two  successive  groups  is 


2n 


k  -  k’  ' 


and  the  system  will  shift  through  this  space  in  a  time 


2U 


w  -  w’ 


We  may  then  define  as  the  Rroup  velocity  the  value  w.  .~..w  .  • 

k  -  k’ 

In  the  limit  as  we  bring  k  and  k’  together,  this  reduces  to  b  =  ^  . 

dk 

If  we  wish,  in  this  limit,  to  put  the  group  velocity  in  terms  of  the 
phase  wavelength  for  small  perturbations  about  a  fixed  component,  we 
have 

b  =  d(kc)  =c+k^!=c-Xdc 

dk  dk  dX 


(1) 


In  this  limiting  case,  we  have  the  usual  situation  for  a  dispersed 
train  well  removed  from  the  origin.  The  wavelet  profile  at  any 
point  in  the  group  is  almost  exactly  a  sine  wave,  while  the  whole 
group  merely  looks  like  a  sine  wave  train  which  is  being  progressively 
’  stretched1  as  we  go  along  the  train.  In  this  circumstance,  the 
velocity  of  propagation  of  the  wavelets  is, to  all  practical  intents 
and  purposes,  the  phase  velocity  for  a  sine  wave  of  the  local  period 
Consider,  however,  the  case  where  active  dispersion  is 
taking  place  in  a  complex  waveform.  Individual  wavelets  will  not  be 
of  sine  wave  form,  but  we  may  nevertheless  assign  a  ’period’  and  a 
’wavelength’  to  any  given  wavelet.  In  this  case  it  is  more  instructive 
to  consider  two  expressions  as  below* 

First  we  express  the  rate  at  which  two  successive  wave  crests 
are  separating  due  to  the  dispersion. 


<3X 
dt  + 


4  =  ^ 

dx  dX 


ex 

dx 
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Secondly,  by  definition  we  have  that  the  group  velocity  is  the 
velocity  at  which  a  given  wave  profile  moves.  Thus  in  the  region 
of  a  geometrical  point  travelling  with  the  group  velocity  b,  the 
’apparent  period’  or  ’group  period’  or  ’group  wavelength’  is  con¬ 
stant 


§\  .  Jh  _  n 
at  0 

In  general,  there  is  a  difference  in  interpretation  between  the  two 
’wavelengths’  used  in  (2)  and  (3).  In  the  limiting  case,  we  are 
again  led  to  (l)  by  combining  (2)  and  (3)  and  submerging  the  dif¬ 
ference  in  the  meaning  of  X. 

Brillouin  (i960)  has  shown  that,  if  we  wish  to  be  precise 
in  our  discussions,  we  should  distinguish  between  signal  velocity, 
front  velocity,  group  velocity,  wave  velocity  and  phase  velocity  - 
all  are  distinct. 

Evidently,  as  stated  by  Stoker  (1957)  it  is  important  for 
the  usual  discussion  of  group  velocity  that  the  motion  discussed 
should  consist  of  a  superposition  of  waves  differing  only  slightly 
in  frequency  and  amplitude.  In  practice,  particularly  for  the 
strong  motion  region,  the  motions  obtained  are  the  results  of  the 
superposition  of  waves  whose  frequencies  range  from  zero  to  infinity 
and  whose  amplitudes  also  vary  widely. 

As  Stoker  points  out,  although  it  is  usual  to  describe  the 


(3) 


transmission  of  energy  in  terms  of  the  group  velocity,  this  is  logical 
only  in  the  limiting  case,  and  the  kinematic  concept  of  group  velocity 
is  of  primary  significance. 
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Inspection  of  a  seismogram  may  easily  lead  to  wrong  impressions 
about  the  actual  ’ground  profile’,  in  the  case  of  a  seismic  surface  wave, 
even  in  the  case  of  a  perfect  displacement  seismograph.  For  example,  a 
constant  period  in  the  seismogram  does  not  of  necessity  imply  a  constant 
wavelength  in  the  ground  profile.  The  effective  wavelength  of  the 
profile,  considered  as  an  entity  passing  the  seismograph,  is  given 
by  the  product  of  the  period  and  the  group  velocities  at  individual 
parts  of  the  train.  If  we  assume  the  duration  of  the  train  is  short 
compared  with  the  travel  time  to  the  seianograph,  the  group  velocity 
decreases  monotonically  as  we  move  along  the  train.  In  the  strong 
motion  region  the  effect  is  complicated  by  the  fact  that  the  travel 
time  is  usually  not  long  compared  with  the  duration  of  the  source 
motion.  For  example,  at  2000  ft  we  may  have  a  train  with  two  to 
three  seconds  duration,  with  a  travel  time  to  the  start  of  the  train 
of  about  two  seconds.  If  we  treat  the  whole  train  as  a  ’wave  packet’ 
we  may  talk  of  a  group  velocity  of  about  $00  fps.  We  would  not  be 
justified  in  quoting  group  velocities  for  the  start  and  end  of  the 
train,  calculated  as  the  epicentral  distance  divided  by  the  travel 
time  from  zero  to  the  start  and  end  of  the  train.  This  would  be 
acceptable  for  the  front  of  the  train,  but  the  travel  time  to  the 
end  of  the  train  is  unknown,  since  the  duration  of  the  originating 
motion  is  unknown,  and  may  well  be  of  the  same  order  as  the  train 
duration . 

Obviously,  in  such  a  discussion  we  must  ignore  the  source 
motion,  at  least  in  the  first  instance.  We  must  consider  the  seismogram 
at  the  seismograph  closest  to  the  source  as  the  signature  of  the  ground 
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motion  at  that  point,  and  then  investigate  the  transmission  of  this 
motion  from  this  point  outwards.  Of  course,  having  at  a  later  stage 
succeeded  in  describing  the  change  in  profile  with  distance,  we  can 
work  backwards  to  the  source.  This  is  essentially  Sato’s  method  of 
analysing  dispersed  surface  waves.  Care  should  also  be  taken  to 
distinguish  the  ’profile’,  worked  out  from  a  seismogram  at  a  fixed 
range  r,  from  the  ground  profile  existing  at  time  t,  equal  to  the 
arrival  time  of  the  first  motion.  Where  dispersion  is  slow,  the  two 
profiles  are  virtually  identical,  but  they  differ  quite  markedly  in 
the  case  of  a  rapidly  dispersing  train. 

In  many  practical  cases  there  is  an  additional  complication, 
due  to  the  fact  that  the  recorded  motion  may  be  compounded  of  more  than 
one  mode,  travelling  according  to  different  dispersion  curves.  Close 
to  the  source  the  ’wave  packet’  will  appear  to  be  a  single  entity. 

As  we  proceed  outwards,  however,  we  note  that  the  wave  packet 
separates  into  two  or  more  groups,  possibly  similar  in  character, 
which  diverge  more  and  more.  The  quiescent  region  between  the  groups 
is  actually  composed  of  a  superposition  of  short-period  waves  being 
shed  by  the  leading  pulse  and  long-period  waves  being  shed  by  the 
following  pulse.  Seismograms  89D-92D  furnish  the  first  example  of 
this  type  of  effect  recorded  in  the  present  research. 

In  the  special  case  where  air-coupling  is  occurring, 
additional  complications  are  introduced  by  the  fact  that  two  essen¬ 
tially  different  wave  groups  are  being  superimposed  one  upon  the  other. 
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We  may  use  this  to  illustrate  our  discussion  graphically. 

Let  us  assume  a  hypothetical  medium  in  which  the  phase 
velocity  varies  with  phase  period  in  the  way  shown  in  Fig.  10. 
Provided  we  are  considering  groups  in  which  only  a  small  range 
of  frequencies  is  involved,  we  may  derive  a  group  velocity-period 
curve  from  this  phase  velocity  curve  by  differentiation.  This, 
however,  presupposes  that  the  actual  group  is  approximately 
sinusoidal  in  wave  form,  so  that  we  can  assign  a  ’group  period’ 
to  the  wavelets.  The  physical  reason  for  the  difference  between 
the  phase  velocity  of  a  given  period  and  the  group  velocity  of 
the  particular  region  in  the  wave  packet  in  the  region  of  this 
period  is  that  the  amplitude  only  of  the  wavelets  change  sharply 
through  the  group.  In  a  region  where  the  train  is  of  reasonably 
constant  amplitude  and  only  slowly  varying  period,  the  phase  velocitjr, 
wave  velocity  and  group  velocity  are  very  closely  similar.  In  the 
strong  motion  region  this  is  not  the  case.  Rapid  dispersion  is 
taking  place,  and  the  wave  group  is  neither  approximately  sinu¬ 
soidal  in  form  nor  constant  in  amplitude,  even  from  one  crest  to 
the  next  adjacent.  Thus,  close  in  there  is  little  meaning  in 
assigning  a  ’period’  to  a  specific  group  velocity.  Let  us,  however, 
move  rather  farther  out  so  that  some  appreciable  velocity  sorting 
has  occurred,  so  that  the  individual  wavelengths  look  roughly 
sinusoidal,  and  we  can  assign  a  period  to  each  wavelet.  Generally 
speaking,  we  shall  see  that  the  early  part  of  the  train  shows  longer 
periods  than  later  parts.  If  we  consider  two  seismograms  from 
stations  sufficiently  close  together  that  there  is  a  recognizable 
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relationship  between  the  two  profiles,  we  see  that  the  periods  associated 
with  comparable  regions  of  the  seismograms  have  changed.  Thus  it  is  not 
feasible  to  derive  ’phase  velocity’  directly  from  the  record.  We  can, 
however,  assign  group  velocities  to  the  comparable  regions,  but  the 
group  velocities  for  the  comparable  regions  will  be  slightly  different, 
since  the  assigned  periods  are  slightly  different.  From  two  such 
seismograms  we  could  plot  a  curve  such  as  that  in  Fig.  10,  where  the 
crosses  refer  to  values  calculated  at  one  station  and  "dots”  values 
calculated  at  the  other.  Each  pair  of  crosses  and  "dots"  is  calculated 
for  the  ’same’  wavelet,  but,  due  to  dispersion,  the  values  of  group 
velocity  and  period  will  both  differ.  A  practical  illustration  is 
afforded  by  seismograms  59W  and  60W,  in  the  region  ahead  of  the  air- 
blast.  The  first  surface  wave  crest  in  59W  has  a  ’group  velocity’  of 
1600  fps,  while  the  "same"  crest  in  60W  has  a  group  velocity  of  1$00  fps. 
The  wavelet  has  a  wave  velocity  of  2560  fps  -  but  this  cannot  be  assigned 
as  the  phase  velocity  of  any  specific  period.  What  can  we  mean  by 
’wavelength’  in  this  case?  Obviously,  we  can  contrive  to  compute  at 
least  three  different  wavelengths,  at  any  one  station.  First,  we  can 
consider  the  wavelength  defined  by  the  wavelet  speed  at  that  point 
multiplied  by  the  period.  This  has  a  distinct  physical  interpretation. 
Similarly,  we  can  compute  a  wavelength  for  the  actual  ground  profile 
when  the  wavelet  is  just  at  its  peak  at  the  seismogram.  This  will  be 
the  group  velocity  multiplied  by  the  period.  Finally,  if  we  knew  the 
phase  velocity  of  a  Fourier  component  with  this  period,  we  could  compute 
a  wavelength  given  by  multiplying  the  measured  period  by  the  phase  velocity. 
Does  this  latter  have  any  physical  meaning?  The  writer  appears  to  differ 
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sharply  in  opinion  here  with  Dobrin  et  al.  (1954) ,  also  quoted  in 
EJP  page  208,  who  show  a  group  \ravelength  computed  in  this  way. 

It  is,  of  course,  merely  a  matter  of  presentation,  but  it  appears 
devoid  of  physical  meaning  to  the  present  writer.  We  may  avoid  this 
difficulty  -  -which  is  really  only  semantic  -  by  confining  our  atten¬ 
tion  to  the  periods  as  actually  measured  on  the  seismogram.  If  we 
do  this,  we  may  interpret  physically  precisely  what  happens  in  the 
case  of  air-coupling  to  a  rapidly  dispersing  wave  packet.  It 
should  be  stressed  that  in  the  following  discussion  we  are  dealing 
with  the  coupling  of  an  ’  air- induced’  train  of  waves  to  a  dispersed 
train  of  waves  which  is  physically  present  in  the  ground.  This 
differs  quite  markedly  from  the  case  where  the  only  train  existing 
is  the  air- coupled,  or  air- induced,  train  itself. 

Referring  again  to  Fig.  10,  let  the  horizontal  line 
marked  ’sonic’  be  the  constant  sonic  velocity  in  the  air  above 
the  medium,  and  let  us  assume  that  a  pressure  pulse  moves  in  the 
air  above  the  surface  at  this  velocity,  starting  simultaneously 
with  the  generation  of  a  strong  motion  pulse  in  the  ground  at  the 
same  origin.  We  attempt  to  predict  what  will  happen  at  two  points, 
one  comparatively  close  and  one  comparatively  far  from  the  source 
(illustrative  examples  will  indicate  the  interpretation  of  close 
and  far) . 

Consider  first  the  point  close  to  the  origin.  The  seismo¬ 
graph  will  first  detect  long  period  waves  arriving  at  high  velocity, 
corresponding  to  point  A  in  the  group  velocity  curve.  Successively 
shorter  period  waves  will  arrive  with  successively  lower  velocities. 
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the  periods  and  velocities  following  the  group  velocity  curve.  When 
we  reach  point  B  of  the  group  velocity  curve,  however,  the  air blast 
wave  arrives  at  the  seismograph.  The  period  in  the  dispersed  wave 
is  then  T^.  However,  there  is  a  phase  wave  of  period  T2  which  has 
the  same  phase  velocity  (point  G  in  Fig.  10).  The  air  blast  will  have 
been  generating  this  phase  wave  continuously  as  it  passed  over  the 
surface,  and  therefore  an  almost  pure  sine  wave  of  this  period  will 
begin  passing  through  the  seismometer  position  simultaneously  with 
the  arrival  of  the  airblast.  Thus  in  the  region  immediately  after 
the  arrival  of  the  airblast,  the  seismogram  wall  be  a  composite 
waveform  of  two  apparent  periods,  Tq  and  T£j  with  the  added  period  T2 
considerably  shorter  than  the  period  of  the  original  ground  wave. 
Seismogram  66W  furnishes  a  good  example  of  the  abrupt  arrival  of  a 
relatively  short  period  wave  in  this  way.  The  nature  of  the  wave 
form  in  the  succeeding  time  interval  will  depend  upon  the  relative 
amplitude,  periods,  and  phase  relationship  of  the  two  superimposed 
wave  groups.  If  the  air-induced  wave  is  dominant,  a  virtually  con¬ 
stant  period  wave,  of  period  close  to  T2  and  virtually  constant 
amplitude  will  persist  until  a  time  corresponding  to  the  point  D 
in  Fig.  10,  which  represents  a  point  in  the  normal  dispersed  wave 
which  has  a  group  period  equal  to  the  period  of  the  phase  wave  of 
the  sonic  velocity.  The  actual  duration  of  the  train  will  thus 
depend  upon  the  distance  from  the  source  and  be  proportional  also 
to  the  intercept  CD.  Thus,  at  a  range  R  the  duration  of  the  train 
will  be  given  by 
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At  this  time  the  air- induced  wave  will  stop,  and  the  dispersed  train 
will  continue  with  successively  shorter  periods  along  the  group 
velocity  curve.  The  situation  is  clearly  evident  in  seianogram  66W. 

Now  consider  a  point  considerably  farther  from  the  source. 
The  duration  of  the  dispersed  train  ahead  of  the  airblast  will  be  much 
longer,  since  the  long  period  waves  will  be  continually  gaining  ground 
on  the  airblast.  The  situation  at  the  arrival  of  the  airblast  will  be 
similar  to  the  previous  case,  but  the  duration  of  the  composite  train 
will  be  longer.  The  situation  changes  if  we  care  to  consider  an  air 
pulse  of  decreasing  velocity.  In  this  case  the  airblast  will  be 
associated  with  gradually  decreasing  periods  in  the  dispersed  train, 
that  is  to  say,  the  airblast  will  suffer  phase  retardation  compared 
with  the  dispersed  train.  Also,  of  necessity,  the  intensity  of  the 
airblast,  and  thus  of  the  air- induced  train,  will  decrease  with 
distance.  Eventually  we  shall  approach  the  true  sonic  velocity,  and 
the  previous  situation  will  be  restored,  with  the  exception  that  the 
effect  of  the  airblast  may  be  either  the  sole  remaining  effect,  or 
alternatively,  may  be  negligible  compared  to  the  amplitude  of  the 
dispersed  wave.  As  we  shall  see  in  the  more  detailed  discussions 
of  the  individual  seismograms,  these  situations  are,  in  fact,  met 
in  practice. 

Finally,  provided  we  have  seismograms  from  at  least  two 
positions  on  a  radial  line  from  the  epicentre,  we  may  apply  SatoTs 
method  and  calculate  the  actual  phase  velocities  for  various  com¬ 
ponents.  In  Appendix  E  a  modified  version  of  Sato’s  method,  suitable 
for  our  present  purposes  and  adapted  to  the  Stantec  Zebra  computer 
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is  described  in  detail.  The  method  is  applied  to  certain  of  the 
seismograms  from  the  present  series  in  Section  3*50  and  elsewhere. 
3.40  The  Surface  Wave  Preceding  the  Airblast 

We  may  now  consider  the  strong  motion  region  of  the  seismo¬ 
grams,  bearing  the  above  discussion  in  mind.  We  deal  first  with  that 
part  of  the  surface  wave  which  is  recorded  before  the  arrival  of  the 
airblast  at  a  given  seismograph.  In  36D  and  37D  the  main  trough  of 
this  phase  is  marked  TBT,  and  we  confine  our  attention  to  the  region 
between  the  dubious  event  mentioned  earlier  and  the  arrival  of  the 
airblast.  This  region  has  been  discussed  in  some  detail  in  a  pre¬ 
liminary  report  (Jones,  Reiniger  and  Cyganik,  1962)  and  the  descrip¬ 
tive  part  of  this  section  is  based  to  some  extent  upon  that  report. 
The  report  includes  a  digitization  of  this  region  for  the  10, OCX)  lb, 
40,000  lb  and  200,000  lb  seismograms,  which  is  too  voluminous  to 
include  here. 

We  start  with  a  phenomological  description  of  the  region, 
based  upon  the  data  quoted  in  the  cited  report. 

Obviously,  from  the  geometry  of  the  situation,  the  records 
represent  surface  waves  travelling  outward  from  the  source,  with  a 
gradual  decrease  in  group  velocity  towards  the  local  air-pressure- 
pulse  velocity.  Insufficient  cycles  are  available  for  a  detailed 
dispersion  analysis,  but  the  trend  may  be  seen  in  Figs.  11  and  12 
for  the  two  larger  charge  sizes.  These  are  plotted  for  the  vertical 
component  only,  since  this  component  is  somewhat  more  sinusoidal 
than  the  radial  and  gives  a  more  consistent  picture.  Group  velocity 
and  periods  for  this  purpose  were  calculated  for  each  crossing  of 
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the  base  line  and  for  each  crest  and  trough. 

3.42  Turning  to  the  trajectories  for  the  region,  it  will 
be  seen  that,  in  the  vertical  trace,  there  is  first  a  slow  rise  in  the 
ground  level,  followed  by  an  abrupt  downward  movement.  In  the  radial 
trace,  a  slow  initial  movement  towards  the  charge  is  followed  by  a 
strong  movement  away  from  the  charge.  It  is  possible  that  the  true 
start  of  the  surface  wave  is  the  strong  down-away  motion,  and  that 
the  up- towards  motion  is  associated  with  the  previous  event. 

Figs.  13  and  14  show  the  nine  trajectories  obtained  from 
10,000  lb  trials  on  the  Drowning  Ford  site.  The  arrival  of  the 
airblast  is  the  last  point  on  the  trajectory.  In  every  case  the 
trajectory  is  a  horizontally  elongated  prograde  ellipse,  about  3/4 
of  a  cycle  at  the  shortest  distance  and  about  11/2  cycles  at  the 
longest  distance. 

Fig.  15  shows  the  trajectories  from  two  10,000  lb 
trials  on  Watching  Hill.  The  orbits  are  not  so  clearly  defined 
as  the  previous  ones,  but  they  do  indicate  a  prograde  elliptical 
orbit . 

Fig.  16  shows  five  trajectories  obtained  with  a  40,000  lb 
trial  on  Watching  Hill.  At  3000  ft  and  4000  ft  the  orbits  are 
single  retrograde  cycles,  approximating  horizontal  ellipses  with  a 
2:1  ratio  horizontal-vertical  axis.  At  5000  ft  and  7000  ft,  the 
orbit  is  still  generally  retrograde,  but  the  ellipse  has  narrowed 
and  tilted  so  that  the  motion  is  approximating  to  an  oscillation 
at  an  angle  of  30°,  the  high  point  being  away  from  the  charge. 

At  10,000  ft  the  orbit  has  changed  radically  into  a  single  prograde 
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cycle  similar  to  those  obtained  with  the  smaller  charges. 

Figs.  17,  1$  and  19  show  ten  trajectories  obtained  with 
the  200,000  lb  charge  on  Watching  Hill.  The  2000  ft  trajectory  is 
that  obtained  with  the  7.5  magnification  instrument,  and  is  possibly 
from  as  high  a  pressure  level  as  any  hitherto  obtained.  The  orbit 
is  by  no  means  clear,  and  it  is  probably  best  to  interpret  it  as 
two  independent  components,  a  general  trend  towards  elevation  of 
the  ground  together  with  an  oscillatory  radial  motion. 

At  4000  ft  a  distinct  retrograde  cycle  occurs,  rather 
more  circular  than  those  at  3000  ft  and  4000  ft  from  the  40,000  lb 
charge. 

At  6000  ft  the  orbit  is  basically  retrograde,  but  has 
narrowed  to  an  almost  linear  45  degrees  motion  similar  to  that  at 
5000  ft  from  the  40,000  lb  charge.  At  7000  ft  three  trajectories 
are  available  from  instruments  disposed  in  azimuth.  All  three  give 
a  consistent  picture  of  a  narrow,  elliptical  orbit  sloped  away  at 
45°  from  the  charge,  and  with  a  faint  tendency  to  convert  to 
prograde  motion. 

At  10,000  ft,  the  orbit  is  definitely  circular,  prograde 
motion,  about  2  1/4  cycles  being  present.  The  early  part  still  shows 
a  tendency  to  a  45°  slope,  but  by  the  end  of  the  motion  the  orbit 
is  almost  circular. 

At  12,000  ft,  two  trajectories  from  instruments  disposed 
in  azimuth  are  available.  Both  show  distinctly  prograde  motion, 
but  the  one  instrument  indicates  a  rather  flat  ellipse,  while  the 
other  shows  an  initially  circular  motion  becoming  elliptical  with 
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a  tilt  of  45° ,  the  high  part  being  towards  the  charge. 

At  15,000  ft  the  trajectory  consists  of  three  flat  lying 
prograde  cycle s,  a  natural  extension  of  the  orbit  at  12,000  ft  on 
the  same  line. 

Fig.  20  shows  the  particle  trajectories  plotted  from  the 
seismograms  obtained  at  distances  of  3000,  4000,  5000,  6000  and 
7000  ft  from  a  charge  of  40,000  lb  on  the  Drowning  Ford  site.  These 
may  be  compared  with  the  trajectories  obtained  from  the  similar  charge 
on  the  Watching  Hill  site  (Fig.  16) .  Although  the  two  test  sites  are 
not  identical,  it  is  obvious  that  a  very  similar  sequence  of  events  is 
indicated  by  the  trajectories.  The  trajectories  at  short  ranges  again 
indicate  a  retrograde  mode,  and  the  trajectories  from  the  region 
5000  -  7000  ft  indicate  a  transition  zone,  with  a  tendency  to  prograde 
motion.  No  seismogram  is  available  from  a  range  of  10,000  ft  on  the 
Drowning  Ford  site,  xdiere  we  would  expect  pure  prograde  motion. 

In  summary,  it  may  be  said  that  these  trajectories  indicate 
that  close  to  the  charge  the  motion  is  retrograde- elliptical.  At 
intermediate  distances  there  is  a  transition  zone  where  the  orbit 
narrows  to  a  linear  trajectory  sloped  away  from  the  charge.  Beyond 
this  point,  the  orbit  T crosses  over’  and  becomes  prograde,  initially 
a  narrow  ellipse  which  opens  out  into  a  circle  and  then  reduces  to  a 
flat  lying  ellipse  which  continues  undisturbed  and  becomes  more  clean 
in  character  up  to  the  arrival  of  the  airblast  at  the  site. 

3.43  Sezawa  and  Kanai  (1935a,  b)  showed  that  there  are 
generally  two  dispersion  curves  for  Rayleigh  Waves  in  a  stratified 
body,  and  in  1939  they  reconsidered  the  subject  and  showed  that  the 
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orbital  motion  was  reversed  between  the  two  cases.  The  fastest 
travelling  (M2)  wave,  according  to  these  authors,  should  show  a 
prograde  orbit,  while  the  normal  (Ml)  Rayleigh  waves  should  show 
the  characteristic  retrograde  orbit.  Nagamune  (1956),  however,  while 
agreeing  with  Sezawa  and  Kanai  in  general,  does  not  agree  that  a 
prograde  orbit  is  typical  of  the  142  branch.  His  calculations  for 
several  specific  cases  show  that  the  M2  branch  may  itself  be 
characterized  by  either  prograde  or  retrograde  orbits,  though  the 
ML  branch  will  always  show  a  retrograde  motion.  A  reversal  in  the 
M2  branch  may  occur,  due  to  changes  in  the  detail  of  the  layering 
along  the  path. 

Howell  (1949),  in  an  investigation  of  the  surface  waves 
from  very  small  charges,  divided  his  record  into  three  regions 
C,  Hand  R.  His  region  G  appears  to  correspond  with  the  region  now 
under  discussion.  Howell  measures  what  he  terms  the  ’group’ 
velocity  of  the  end  of  the  C  region  (and  the  start  of  the  H  region) 
and  quotes  345  metres/sec.,  the  sonic  velocity.  Haskell  (1951) 
discusses  this  paper  in  connection  with  air-coupling  and  points 
out  that,  in  fact,  Howell  was  measuring  the  phase  velocity  .  The 
R  region  is  stated  to  be  the  Airy  phase. 

Leet  (1939  and  1946)  discusses  a  surface  wave  with  prograde 
orbital  motion,  which  he  termed  the  Hydrodynamic  wave,  H.  Care  should 
be  taken  to  distinguish  between  the  Leet  and  Howell  H  waves,  which  are 

Neither  statement  appears  to  be  quite  correct  -  see  previous 

discussion. 
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quite  distinct.  Leet’ s  recording  of  prograde  orbital  motion  is  believed 
to  be  the  earliest  in  the  literature,  but  apparently  Leet  has  at  no  time 
suggested  that  this  might  be  the  Sezawa  and  Kanai  M2  wave. 

In  1953?  Howell,  Neuenschwander  and  Pierson,  in  an  investi¬ 
gation  of  surface  waves  in  the  Gulf  Coast,  showed  that  their  records 
could  be  divided  into  two  regions,  in  the  first  of  which  the  motion 
was  pro grade  and  in  the  second,  retrograde,  as  called  for  in  the  Sezawa 
and  Kanai  theory. 

The  Ml  and  M2  waves  are  frequently  termed  the  symmetric  and 
antisymmetric  branches,  by  analogy  with  the  vibrations  of  a  plate. 
Tolstoy  and  Usdin  (1953)?  noting  PekerisTs  (1948)  interpretation  of 
normal  mode  propagation  in  a  layer  as  an  interface  or  ’wave  guide’ 
phenomenon,  developed  a  very  elegant  treatment  by  ray  theory.  The  Ml 
and  M2  branches  appear  quite  simply  as  a  consequence  of  this  theory, 
and  the  theory  confirms  the  generality  of  Sezawa  and  Kanai  (and 
Nagamune’s)  numerical  treatment.  They  show  that  the  frequency  equation 
can  be  derived  from  ray  theory  in  a  far  simpler  way  than  from  the 
usual  method  of  expanding  the  determinant  (illustrated  in  Appendix  A). 

In  the  Tolstoy  and  Usdin  method,  the  form  of  the  equation  is  entirely 
physical,  and  each  coefficient  may  be  interpreted  in  terms  of  the 
reflection  coefficients  at  the  boundaries. 

A  superficially  similar  set  of  trajectories,  with  directional 
reversal  following  a  linear  transition  phase,  may  be  obtained  in  the 
case  of  the  normal  (ill)  wave  simply  by  tracing  the  trajectory  changes 
with  depth  at  a  given  point.  This  has  been  done,  for  example,  by 
Dobrin,  Simon  and  Lawrence  (1951)*  This  similarity  is,  of  course. 
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quite  superficial  and  the  mechanism  is  entirely  different  in  the  two 
cases. 

None  of  the  orbits  obtained  in  this  section  has  the  canonical 
form  of  a  Rayleigh  orbit  with  the  radial  axis  about  two-thirds  of  the 
vertical.  This  departure  from  the  canonical  form  is  usually  explained 
by  departures  from  the  Poisson  condition  (in  which  the  two  Lame  elastic 
constants  are  numerically  equal)  and  possibly  by  visco-elastic  effects. 
Horton  (1953)  has  investigated  Rayleigh  waves  on  a  visco-elastic  solid 
and  has  shown  that,  while  there  is  a  tendency  to  convert  from  a  verti¬ 
cal  ellipse  to  a  Lissajous  figure  tilted  towards  the  source,  the  magni¬ 
tude  of  this  effect  is,  in  theory,  very  much  less  than  the  experimentally 
observed  departures  from  the  canonical  state. 

Attention  is  also  drawn  to  the  theoretical  treatment  of 
Kuo  (195$)  of  the  propagation  of  Rayleigh- type  pulses  on  a  visco-elastic 
solid.  He  showed  that  competent  rocks  may  be  treated  as  purely  elastic, 
but  that  unconsolidated  sediments  should  be  treated  as  visco-elastic. 
However,  his  treatment  does  not  allow  for  a  departure  of  the  orbit  from 
the  canonical  retrograde  elliptical  form.  Kuo  also  discusses  the  experi¬ 
mental  verification  of  the  existence  of  Ml  and  M2  waves,  and  confirms 
that  the  major  axes  of  the  two  modes  are  tilted  relatively  to  each  other. 
In  the  present  region  of  the  seismograms  (that  is,  ahead  of  the  airblast), 
it  is  difficult  to  compare  the  results  with  the  limited  results  of  Kuo, 
and  it  seems  to  the  writer  that  the  gradual  rotation  of  the  major  axes 
in  these  cases  is  merely  due  to  the  proximity  of  the  transition  zone 
between  the  pulses.  Well  removed  from  the  transition  zone,  the  trajec¬ 
tories  appear  to  be  consistent  in  form,  as  may  be  seen  by  comparing  the 
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200,000  lb  trajectories  at  12,000  and  15,000  ft.  As  will  be  seen  in  a  later 
section,  however,  the  results  of  Kuo  are  confirmed  and  extended  for  buried 
charges. 

To  summarize  the  discussion  of  this  section,  we  may  therefore 
say  that  the  region  encloses  an  M2  vibration,  with  a  prograde  orbital  motion. 
However,  contrary  to  the  theory  of  Sezawa  and  Kanai,  for  a  single  model,  the 
prograde  M2  branch  arrives  after  a  retrograde  branch.  It  would  appear  that 
this  retrograde  branch  cannot  be  identified  as  the  equivalent  mode  of  the 
Ml  branch  in  the  same  system,  since  this  branch  would  arrive  after  the  M2. 

It  is  possible  that  we  are  dealing  here  with  a  transition  between  modes  of 
different  systems,  one  of  which  has  retrograde  and  the  other  prograde 
motion.  Should  this  be  indeed  the  case,  we  would  expect  that  the  prograde 
orbit  would,  in  due  course,  suffer  a  transition  into  a  retrograde  orbit. 
However,  the  prograde  M2  mode  is  interrupted  by  the  arrival  of  the  air- 
blast  and  we  must  therefore  defer  investigation  of  this  point  until  after 
a  discussion  of  the  event  associated  with  the  airblast. 

The  discussion  of  this  section  has  concentrated  on  the  three 
larger  charge  sizes.  It  will  be  seen  that  this  is  a  necessary  restriction, 
since  in  the  next  size  down  (30D  and  32D)  the  epicentral  distance  is  too 
small  to  allow  this  train  of  events  to  be  separated  from  the  preceding  P 
wave  and  the  following  airblast. 

Attention  is  drawn  to  the  work  of  Kisslinger  (1959),  (and 
Appendix  D),  using  similar  instrumentation,  and  small  charges,  at  Saint 
Louis.  There  is  a  close  similarity  between  his  buried  charge  records 
and  the  region  presently  under  consideration,  and  he  obtained  a  similar 
separation  into  prograde  and  retrograde  orbits. 


'■ 

l>  ■_ 

_ 


. 


-  59  - 


3.50  The  Surface  Waves  From  Contained  and  Venting  Charges 

3.51  It  is  convenient  at  this  point  to  discuss  the  Rayleigh 
and  Sezawa-type  surface  waves  produced  by  charges  detonated  in  a  geometry 
different  from  that  of  the  main  series.  These  charges  enable  us  to  look 
at  the  dispersed  surface  waves  and  the  air-coupled  waves  separately,  and 
thus  study  the  sequence  of  prograde  and  retrograde  phases  and  the  method 
of  coupling  to  the  airblast. 

3.52  Records  were  obtained  from  four  charges  of  1000  lb  TNT 
fully  contained  at  a  depth  of  50  ft  on  the  flank  of  Watching  Hill.  All 
these  charges  were  in  the  form  of  cased  cylinders,  6  ft  long  by  1$  ins 
in  diameter,  inserted  in  20  inch  diameter  wells  53  ft  deep.  Two  of  the 
charges  were  water-tamped,  and  two  fully  tamped  with  silt.  No  appreci¬ 
able  difference  was  detected  which  could  be  attributed  to  the  method  of 
tamping,  either  in  the  source  region  or  at  remote  points  where  the 
waves  were  monitored  by  U.  of  A.-SES  MDHO  project  teams. 

We  consider  first  a  pair  of  charges  (F.E.  549-1  and  F.E.  549-2) 
which  were  buried  at  50  ft  and  on  the  surface  respectively.  Records  100H- 
103H  were  obtained  at  distances  of  4000,  6000,  $000  and  10,000  ft  from 
the  GZ  of  the  buried  charge.  Attention  is  drawn  to  the  close-in  record, 
100H.  Cursory  examination  shows  that,  in  the  vertical  trace,  several 
phases  are  present,  and  the  train  ends  with  a  dispersed  train  of  decreas¬ 
ing  amplitude.  With  the  exception  of  the  final  region,  the  phases  show 
up  with  almost  equal  clarity  in  the  radial  and  transverse  traces. 

The  second  charge,  a  1000  lb  sphere  of  TNT,  was  detonated  on 
the  same  seismic  line,  on  the  surface,  but  some  150  ft  from  the  GZ  of 
the  buried  charge,  and  records  122H-126H  obtained  at  the  same  recording 
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sites  as  100H-103H.  Comparison  of  the  two  sets  of  seismograms  indi¬ 
cates  that  this  charge  did  not  feed  sufficient  energy  into  the  ground 
to  excite  the  dispersed  surface  wave  phases  with  measurable  amplitude. 
Immediately  following  the  airblast,  however,  we  see  a  wave  train  of 
approximately  constant  period,  though  not  of  particularly  sinusoidal 
form.  The  vertical  components  of  the  two  seismograms  at  a  nominal 
distance  of  4000  ft  from  these  two  charges  are  re-plotted  in  Fig.  21. 

If  we  assume  that,  in  the  case  of  the  surface-burst  charge,  the  whole 
of  the  record  is  a  Press-Ewing  coupled  wave,  we  may  estimate  the  phase 

velocity  from  the  duration  by  using  the  relationship  Vc  =  R  +  Vg  given 

T 

earlier  (section  3.31*  p.  49).  The  values  obtained  in  this  way  are, 
from  the  2000  ft  seisnogram,  122H:  T  =  2.60,  Vg  =  1280  and  hence 
Vc  =  2050  fps.  From  the  4000  ft  seismogram,  123H,  we  obtain  T  =  6.14* 
Vg  =  1200  and  hence  Vc  =  1855  fps.  It  will  be  seen  later  in  this 
section  that  these  velocity  pairs  are  reasonable  values  for  Press- 
Ewing  type  coupling  to  an  M2  mode  on  this  test  site.  The  distance  of 
4000  ft  corresponds  to  a  distance  of  24*000  ft  from  the  200,000  lb 
charge,  as  far  as  the  airblast  is  concerned,  and  so  is  well  outside  the 
strong  motion  region  investigated  in  the  case  of  the  larger  charges. 

It  is  of  interest  to  note  that  the  form  of  the  air-induced  train 
obtained  in  this  case  is  very  similar  to  the  records  obtained  at  the 
greatest  distances  from  the  larger  charges,  where  the  pre-airblast 
wave  has  fallen  below  the  recording  capability  (cf  42D,  43D  and  122H) . 
Anticipating  the  dispersion  curves  discussed  later  in  this  section,  we 
may  say  that  Press-Ewing  coupling  occurs  at  points  well  removed  from 
the  charge. 
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We  may  use  the  two  traces  plotted  in  Fig.  21  to  illustrate 
graphically  the  appearance  of  a  wave  which  is  a  combination  of  a 
Press-Ewing  train  and  a  pre-existing  dispersed  train,  though  of 
course  we  cannot  in  this  way  demonstrate  a  ’resonating’  effect.  The 
air-coupled  train  is  superimposed  directly  onto  the  dispersed  train 
at  a  point  corresponding  to  the  arrival  of  the  airblast  at  the  given 
distance.  The  process  is  then  repeated  with  the  time  of  arrival 
shifted  in  phase  relative  to  the  crest  of  the  dispersed  train.  Three 
examples  of  the  resulting  ’synthetic  seismograms’  are  shown  in  Fig.  22, 
which  is  a  photograph  of  the  computer  printout  of  the  trace  amplitudes 
added  at  intervals  of  0.02  sec.  Due  to  the  comparatively  low  amplitude 
of  the  air-induced  train,  the  effect  is  not  so  striking  as  in  the  multi¬ 
ton  charge  seianograms,  but  there  is  a  marked  similarity  which  indicates 
that  we  may  attempt  to  interpret  the  actual  seismograms  in  this  way. 

3.53  We  may  now  with  advantage  consider  the  seismograms 
obtained  from  charge  F.E.  549-4  (seismograms  94J-9&J).  This  charge 
was  a  1000  lb  sphere  buried  at  a  depth  of  23.5  ft  in  a  low  sandy  hill 
on  the  main  seismic  line  of  the  200,000  lb  charge  (see  para  2. 41). 
Several  distinct  phases  are  visible  in  the  seismograms.  If  only  the 
vertical  component  were  available  (as  is  usually  the  case  in  explora¬ 
tion  work)  one  might  judge  that  only  two  phases  existed.  However, 
when  one  correlates  the  radial  and  vertical  components  as  a  particle 
trajectory,  one  finds  that  there  are,  in  fact,  at  least  five  distinct 
phases,  which  alternate  between  the  retrograde  and  prograde  habit. 

Each  of  these  phases  has  associated  with  it  a  phase  in  the  transverse 
trace.  Fig.  23  shows  the  particle  trajectories  for  the  radial- vertical 
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components  of  seismogram  97<J>  in  -which  the  five  phases  are  clearly 
evident.  At  points  closer  to  the  charge,  this  separation  is  incomplete, 
and  the  interpretation  of  such  seismograms  must  always  be  complicated 
and  of  doubtful  validity.  Where  a  phase  is  clearly  separated,  it  is 
possible  to  identify  it  by  its  particle  trajectory  at  various  distances 
and  then  obtain  the  spectrum  and  analyse  the  dispersion  by  Sat8Ts  method. 
Where  the  phases  are  not  separated,  however,  one  can  be  led  into  serious 
difficulties  by  blind  attempts  to  apply  such  methods.  In  extreme  cases 
we  may  see  the  attempt  made  to  apply  a  Sat6-type  analysis  to  a  complete 
seismogram,  despite  the  fact  that  none  of  the  fundamental  assumptions  is 
valid. 


At  the  degree  of  separation  shown  in  Fig.  23  we  may  select  the 
following  phases? 

Phase  One.  Three-quarters  of  a  cycle,  retrograde  habit, 
flat  lying  ellipse. 

Phase  Two.  (a)  Four  cycles,  prograde  habit,  fairly  narrow 
ellipse  tilted  forward  60  deg.  from  the  vertical. 
This  breaks  into 

(b)  which  may  or  may  not  be  the  same  physical 
phase.  This  is  one  complete  cycle,  which  re¬ 
tains  the  prograde  habit  but  whose  tilt  is 
nearer  45  deg.  forward  of  the  vertical. 

Phase  Three.  Two  cycles,  retrograde  habit,  narrow 

ellipse  tilted  back  25  deg.  from  the  vertical. 
Phase  Four.  Two  cycles,  prograde  habit  which  retains  the 
same  general  shape  and  the  same  backward  tilt 


as  the  retrograde  phase  three 


. 
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Phase  Five.  Four  cycles,  retrograde  habit,  approximately 
vertical  ellipses  of  classical  form.  This 
phase  is  very  small  in  amplitude  compared 
with  the  other  phases. 

The  dominant  phases  are  the  second  and  third,  which  are 
prograde  and  retrograde  and  tilted  forwards  and  backwards  from  the 
vertical  respectively,  so  that  the  included  angle  is  approaching  a 
right  angle.  There  can  be  little  doubt  that  these  two  phases  con¬ 
stitute  an  M2-ML  pair  for  the  same  layering  system.  The  agreement  with 
the  systsn  discussed  by  Kuo  (1953)  is  virtually  exact.  The  other  phases 
are  less  easy  to  identify,  since  changes  of  habit  can  occur  for  a  variety 
of  reasons.  The  M2  motion  itself  may  become  retrograde  for  short  wave¬ 
lengths  in  a  fixed  layering  system,  but  in  such  a  case,  following  Kuo, 
we  would  expect  the  tilt  of  the  major  axis  to  remain  consistently  for¬ 
ward  of  the  vertical.  In  the  present  seismogram  no  retrograde  orbit 
has  the  same  forward  tilt  as  a  prograde  orbit.  It  is  noteworthy, 
however,  that  in  the  trajectory  illustrated,  a  retrograde  orbit  tilted 
backwards  from  the  vertical  transforms  into  a  prograde  orbit  of  similar 
inclination.  At  first  sight  one  might  think  that  this  represents  a 
prograde  Ml  wave,  but  there  is  no  evidence  whatsoever  in  the  literature 
or  in  the  present  work  to  justify  such  a  conclusion.  The  Rayleigh  (Ml) 
mode  is  essentially  retrograde  in  habit,  for  every  system  which  has  been 
calculated.  It  appears  to  be  generally  agreed  that  the  M2  and  Ml  orbits 
tend  to  be  orthogonal,  independent  of  the  fact  that  the  M2  mode  may 
itself  be  retrograde.  Kuo  states  that  the  M2  wave  is  tilted  forward 
of  the  vertical,  and  the  ML  backward,  and  this  conclusion  is  generally 
confirmed  by  the  present  work.  However,  Howell  (1943)  observed  that  in 
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his  experiments  the  two  phases,  while  orthogonal  to  each  other,  were 
tilted  in  the  reverse  direction  -  M2  backward  and  Ml  forward.  Thus  the 
backward  tilted  prograde  wave  of  phase  four  is  almost  certainly  M2,  but 
it  is  not  clear  why  this  phase  should  depart  from  the  common  forward 
tilt  for  the  site.  As  no  theory  exists  to  explain  the  tilted  orbits, 
it  is  premature  to  draw  conclusions.  It  is  eminently  possible  that  the 
direction  of  tilt  is  governed  by  fortuitous  elements  (as  is  the  direction 
of  rotation  of  a  bath  tub  vortex).  Once  established  there  will  be  a 
tendency  for  a  given  habit  to  hold,  even  though  a  new  transient  is  driving 
the  motion,  provided  a  smooth  transition  is  possible  from  one  habit  to  the 
other.  It  is  noteworthy  that  the  tilt  of  the  orbit  does  not  remain  constant 
for  a  given  phase  as  it  proceeds  outwards  from  the  source.  Kuo  indicated 
that  the  tilt  changes  linearly  with  distance,  and  there  is  some  indication 
of  a  similar  effect  in  the  present  results  for  the  prograde  phase.  The 
retrograde  phases  are  less  clearly  developed  and  conclusions  would  be 
premature,  but  no  linearity  is  obvious.  There  is  room  here  for  a  definitive 
investigation. 

We  may  proceed  to  examine  the  dispersion  curves  for  the  two  dominant 
phases  according  to  the  method  described  in  Appendix  E.  It  is  necessary  to 
make  a  subjective  selection  of  the  equivalent  regions  in  the  seismograms  from 
two  different  distances.  Seismograms  97 J  and  9&J  were  selected,  and  the 
corresponding  regions  defined  from  the  particle  trajectories.  These  regions 
were  then  digitized,  and  the  Fourier  Integrals  calculated  separately  for  each 
phase  at  each  station.  The  dispersion  curves  were  then  deduced  as  a  ’sheaf’ 
of  curves  for  each  phase,  due  to  the  indeterminate  multiplier  in  the  equa¬ 
tions  (see  Sato,  1955b). 
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The  resulting  sheaves  for  the  first  attempt  are  illustrated 
in  Fig.  24.  The  curves  for  the  prograde  mode  are  well  separated,  and 
it  is  comparatively  easy  to  select  the  correct  curve.  This  is  parti¬ 
cularly  important  since  the  airblast  arrival  is  usually  associated  with 
the  dominant  prograde  mode.  Selection  of  the  curve  may  be  made  on  the 
basis  of  the  wavelet  velocity  (as  defined  in  para  3.30),  thus  defining 
the  dispersion  curve  with  precision.  A  full  analysis  of  the  records 
from  this  1000  lb  charge  is  being  prepared  by  Hastrup  (1963)  in  which 
each  phase  is  analysed  separately  and  in  detail.  The  basic  conclusion  is 
unchanged  in  this  more  detailed  analysis,  but  there  is  some  evidence  that 
the  spectra  become  more  peaked  as  the  distance  increases  and  this  effect 
is  more  marked  for  the  pro grade  phase  than  for  the  retrograde  phase. 
Hastrup  also  discusses  the  change  in  the  peak  frequency  with  charge 
size  for  the  10  lb,  296  lb  and  1000  lb  charges.  Considerably  more  work 
is  required  in  this  area  and  use  will  be  made  of  later  shots  in  the 
U.  of  A.  — SES  MOHO  programme  to  define  the  dispersion  curves  and  the 
spectra  of  the  phases  in  detail.  The  preliminary  analysis  is  suf¬ 
ficient  to  show  that  the  conditions  for  Press-Ewing  coupling  exist, 
and  that  the  group  velocity  ~  phase  velocity  bracket  obtained  experi¬ 
mentally  with  the  surface-burst  charges  is  compatible  with  the  dispersion 
curve  deduced  from  the  buried  charges. 

*  Records  89D-93 D,  94J-9&J  and  the  buried  10  lb  charges  functioned 
at  SES  for  Dr  Kisslinger,  Saint  Louis  University. 
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EXPERIMENTAL  DISPERSION  CURVES  FOR  PROGRADE 
AND  RETROGRADE  PHASES  OF  A  IOOO  LB.  BURIED  CHARGE 
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3.60  The  Waves  Folio-wing  the  Airblast 

Following  the  arrival  of  the  airblast,  we  may  expect  to  find 
(a)  an  indication  of  the  action  of  the  airblast  on  the  instrument, 
together  with  either  (b)  a  continuation  of  the  virtually  unaffected 
surface  wave  or  (c)  an  aircoupled  wave  followed  by  uncoupled  surface 
waves. 

The  theory  of  aircoupling,  and  a  discussion  of  the  appearance 
to  be  expected  when  an  air-coupled  wave  is  superimposed  upon  a  dispersed 
wave,  is  given  elsewhere  in  this  dissertation.  Attention  is  also  drawn 
to  the  experimental  work  of  Anderson  and  Henriksen  (1962) .  In  this 
section  we  need  only  look  at  the  experimental  results  obtained  on  the 
Suf field  sites.  It  may  be  said  immediately  that  in  all  the  close-in 
records  there  is  some  indication  of  the  airblast  coupling  to  either 
a  prograde  or  retrograde  wave.  In  all  cases  except  the  most  distant 
records,  there  is  a  distinct  change  in  the  record  at  the  arrival  of 
the  airblast.  Preceding  the  airblast,  the  record  is  either  below  the 
recording  capability  of  the  instrument,  or,  alternatively,  is  a  smooth, 
sinusoidal  wave  train  showing  slight  dispersion.  After  the  arrival 
of  the  airblast  there  is,  in  all  cases  where  a  smooth  train  has  preceded 
the  airblast,  a  disturbed  train  following  the  airblast.  In  cases  where 
there  is  no  preceding  wave,  there  is  a  wave  recorded  after  the  arrival 
of  the  airblast.  In  many  cases  (for  example  seismograms  65W  and  66W) 
the  airblast  arrival  is  followed  by  many  cycles  of  virtually  constant 
amplitude,  constant  period  waves,  which  in  the  vertical  component  looks 
like  classical  air- coupling .  We  have  also  shown  (in  para.  3.50)  that 
the  experimental  phase  dispersion  curve  for  a  prograde  mode  on  the 
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Watching  Hill  test  site  crosses  the  sonic  (and  blast)  velocities,  so 
that  Press-E-wing -type  coupling  may  exist.  The  amplitude  of  the  air- 
coupled  wave  cannot  be  predicted,  and  in  individual  cases  it  is  pos¬ 
sible  to  have  either  a  small  amplitude  air-coupled  train  superimposed 
upon  a  pre-existing  surface  wave,  not  necessarily  the  same  mode,  or 
the  pre-existing  surface  wave  mode  driven  into  large  amplitude  reson¬ 
ance  . 

We  may  use,  as  our  first  illustrations,  selected  records 
obtained  with  the  200,000  lb  charge.  We  start  at  12,000  ft  along  the 
A  line,  seismogram  69W,  in  which  the  apparent  effect  of  the  airblast 
arrival  is  negligible  in  the  record.  The  arrival  coincides  with  a 
high  frequency  disturbance  which  is  rapidly  damped  -  a  typical  effect 
at  all  distances,  which  is  attributed  to  the  direct  action  of  the  air- 
blast  on  the  case.  After  this  high  frequency  effect  has  died  out, 
there  is  a  relatively  large  amplitude  surface  wave  which  has  departed 
from  the  pre-airblast  sinusoidal  wave  form.  The  particle  trajectory  is 
shown  in  Fig.  25.  As  indicated  earlier,  the  pre-airblast  orbit  is  a 
very  clean  ellipse,  with  the  major  axis  horizontal  and  a  ratio  between 
the  axes  of  about  two.  The  actual  arrival  of  the  airblast  does  not 
modify  the  orbit  instantaneously,  but  shortly  afterwards  the  orbit  is 
distorted  from  its  elliptical  form,  with  a  tendency  to  indicate  a  tilt¬ 
ing  of  the  major  axis  toward  the  vertical.  This  orbit  then  breaks 
sharply  into  an  approximately  orthogonal  tilt,  and  then  settles  down 
again  into  a  smooth,  almost  circular  prograde  orbit.  One  concludes 
that  there  is  no  direct  evidence  of  Press-Swing  coupling,  but  evidence 
of  two  successive  prograde  habit  phases,  modified  or  induced  by  the 
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Moving  closer  to  Ground  Zero.,  we  look  at  seismogram  6$W.  In 
this  seismogram  there  is  distinct  evidence  of  a  change  of  phase  in  the 
radial  component  at  the  arrival  of  the  airblast,  but  in  the  vertical  com¬ 
ponent  there  is  no  abrupt  change  in  the  wave  profile,  though  the  region 
after  the  arrival  of  the  airblast  is  of  larger  amplitude  and  less  smooth 
than  that  preceding  the  airblast.  Turning  to  the  particle  trajectories., 
we  see  that  the  orbit  preceding  the  airblast  (Fig.  1$  and  start  of  Fig.  26), 
is  a  smooth  prograde  orbit,  initially  elliptical  with  the  major  axis 
approximately  four  times  the  minor,  and  tilted  forward  of  the  vertical,  but 
becoming  virtually  circular  by  the  arrival  of  the  airblast.  Immediately 
following  the  airblast  there  is  horizontal  movement  towards  ground  zero, 
followed  by  a  large  amplitude  orbit  which  is  almost  circular,  but  has  a 
tendency  towards  the  form  of  a  forward- tilted  ellipse.  After  slightly  more 
than  two  cycles,  this  orbit  breaks  sharply  towards  the  vertical  and  then 
sets  into  an  almost  linear  oscillation  tilted  backwards  from  the  vertical, 
almost  orthogonal  to  the  tilt  of  the  pre-airblast  orbit.  In  the  confused 
transition  region  there  is  some  tendency  to  retrograde  motion,  but  the 
dominant  motion  is  prograde.  Thus,  in  this  seismogram  we  again  have 
evidence  of  two  prograde  modes,  generally  circular  in  form  but  with  a 
tendency  to  orthogonality  in  their  elliptical  regions.  There  is  little 
evidence  of  a  superimposed  Press-Ewing  train.  The  sharp  motion  in 
towards  ground  zero  cannot  be  explained  in  terms  of  the  action  of  the 
airblast  upon  the  instrument  itself,  and  is  without  doubt  due  to  the 
influence  of  the  airblast  on  the  ground. 

Moving  still  closer  to  GZ,  we  look  at  seismogram  65W.  This, 
with  66W,  constitutes  a  pair  of  seismograms  which  apparently  show  very 
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marked  air-coupling.  In  65W  we  see  that,  ahead  of  the  airblast, 
there  is  a  strong  surface  wave,  which  shows  a  change  of  phase  in 
the  vertical  approximately  0.3  sec.  ahead  of  the  airblast,  and  a 
very  marked  change  of  phase  at  the  arrival  of  the  airblast.  The 
particle  trajectory  is  shown  in  Fig.  27.  The  main  pre-airblast 
phase  shows  a  retrograde  orbit,  breaking  into  a  prograde  orbit 
0.3  sec.  ahead  of  the  airblast.  The  arrival  of  the  airblast  is 
followed  by  a  large  amplitude,  forward  tilted  prograde  orbit. 

There  are  approximately  three  and  a  half  cycles  of  this  orbit,  with 
evidence  of  a  superimposed  effect  at  the  top  of  the  orbit.  The  or¬ 
bit  finally  narrows  down  to  a  linear  oscillation,  and  then  breaks 
sharply  in  the  vertical  and  becomes  a  new  orbit  orthogonal  to  the 
previous  orbit,  but  still  retaining  the  prograde  habit.  At  a 
later  stage  this  orbit  becomes  circular,  then,  after  a  confused 
region,  settles  down  into  a  retrograde  orbit  tilted  forward  of  the 
vertical.  Thus,  with  the  exception  of  the  final  phase  (which  is 
presumably  below  the  recording  capability  at  the  distant  stations) 
and  the  initial  phase,  both  of  which  are  retrograde,  this  seismogram 
indicates  that  the  pair  of  orthogonal  phases  were  already  a  dominant 
effect  at  this  distance,  where  there  is  strong  coupling  to  the  air- 
blast  in  the  immediate  post-airblast  region. 

At  the  closest  distance,  2000  ft,  seismogram  64W,  we  have 
a  record  from  a  low  gain  instrument.  The  particle  trajectory  is 
shown  in  Fig.  2$.  Once  again,  the  strong  motion  is  associated  with 
the  arrival  of  the  airblast.  The  first  motion  after  the  arrival  of 
the  airblast  consists  of  a  single  elliptical,  backward  tilted  orbit 


V 


Is- 

CVJ 

6 


100  TON  PARTICLE  TRAJECTORIES  —  WATCHING  HILL  SITE 


I 


00 

00 

e> 

Li_ 


100  TON  PARTICLE  TRAJECTORIES  —  WATCHING  HILL  SITE 


-  70  - 


with  a  retrograde  habit.  This  widens  out  into  a  large  amplitude, 
almost  circular  orbit  before  changing  abruptly  into  a  forward-tilted, 
narrow  elliptical  orbit,  converting  into  a  prograde  phase.  Presumably, 
by  a  comparison  of  this  record  and  the  record  at  4000  ft  (65W),  the 
retrograde  mode  is  that  equivalent  to  the  first  prograde  mode  at 
4000  ft.  The  retrograde  mode  at  4000  ft  is  indicated  in  the  2000  ft 
record  by  the  low  amplitude  signal  visible  in  the  radial  trace  between 
the  P  wave  and  the  airblast  arrival. 

We  see  that  these  200,000  lb  seismograms  show  quite  clearly 
a  succession  of  four  phases,  an  initial  retrograde  phase  followed  by 
the  two  dominant  phases  which  are  usually  prograde  and  orthogonal  to 
each  other.  The  prograde  habit,  however,  is  not  invariable,  and  these 
phases  may  adopt  the  retrograde  habit.  Finally,  the  strong  motion  ends 
with  a  relatively  low  amplitude  retrograde  phase. 

The  arrival  of  the  airblast  is  at  all  times  associated  with 
the  (normally  prograde)  second  phase.  The  actual  motion  in  this  phase 
tends  to  start  slightly  ahead  of  the  airblast,  but  at  the  short  epi- 
central  distance  the  association  is  close.  Classical  air-coupling  in 
the  Press-Ewing  sense  does  not  appear  to  be  dominant.  Nevertheless, 
it  would  appear  that  the  air  blast  is  driving  a  surface  wave  -  or  waves  - 
with  large  amplitude.  This  is  made  abundantly  clear  in  the  20  ton  seismo¬ 
grams  59W,  60W,  6lW  and  62W.  Before  we  look  at  the  particle  trajectory 
for  these  waves,  it  is  rewarding  to  compare  the  vertical  phases  of  these 
seismograms.  These  are  re-plotted  in  Fig.  29  j  with  the  airblast  arrival 
aligned.  Clearly,  the  arrival  of  the  airblast  is  driving  the  post-air- 
blast  wave.  A  dispersive  wave  is  gaining  ground  upon  the  airblast,  as 
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FIG.  29 

VERTICAL  COMPONENTS  FROM  20  TON  SEISMOGRAMS 
SHOWING  EFFECT  OF  AIR  BLAST 
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a  smoothly  dispersed  train,  but,  despite  this  ’shedding’,  due  to 
dispersion,  there  is  no  indication  whatsoever  that  the  large  ampli¬ 
tude  post-airblast  train  at  the  shortest  distance  (59W)  is  gaining 
on  the  airblast  even  at  the  greatest  distance  (62W).  A  simple  dis¬ 
persive  train  would  certainly  have  done  so.  The  particle  trajectory 
for  the  3000  ft  seismogram  is  shown  in  Fig.  30.  The  impulsive  action 
on  the  instrument  case  is  clearly  evident  in  the  right-hand  side  of 
the  orbit.  The  orbit  again  shows  a  retrograde-prograde- retrograde- 
prograde  succession,  but  the  phases  are  neither  so  clearly  separated 
nor  so  well  formed  as  in  the  200,000  lb  records  -  possibly  due  to  the 
shorter  epicentral  distances  for  a  given  airblast  pressure. 

Precisely  the  same  association  of  the  airblast  and  a  strong 
post-airblast  surface  wave  is  demonstrated  by  the  seismograms  from 
10,000  lb  charges  on  the  Watching  Hill  site.  Typical  examples  are 
50W  and  51W,  at  2000  and  3000  ft  respectively  on  the  same  azimuth. 

The  association  is  repeated  in  the  2000  -  3000  ft  pair  from  another 
charge  and  different  azimuth  and  ground  zero  in  44W  and  45W,  and  many 
examples  at  2000  ft  in  different  azimuths  (e.g.  47W  et  seq)  indicate 
that  the  association  is  quite  invariable.  There  can  be  no  doubt  of 
the  reality  of  the  phenomenon  on  the  Watching  Hill  test  site,  even 
though,  due  to  the  pre-existing  surface  wave  and  the  multiplicity  of 
phases,  the  classical  Press-Ewing  form  is  not  clearly  evident. 

The  theoretical  dispersion  curves  given  in  Appendix  C 
indicate  that  Press-Ewing  coupling  may  be  expected  for  the  layering 
system  of  the  test  site,  for  waves  in  the  period  bracket  excited  by 
the  multi-ton  charges.  Small  charges,  which  excite  high  frequency 
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waves  (4-10  cps)  will  produce  seismograms  with  low  group  velocities 
and  little  evidence  of  air-coupling.  Larger  charges,  which  excite 
waves  of  longer  period,  will  exhibit  wave  trains  in  which  the  air- 
blast  is  associated  with  wave  groups  of  the  same  velocity  due  to 
the  flat  ’plateau1  region  in  the  group  velocity  curves.  Group 
periods  of  about  0.5  sec.  will  be  associated  with  phase  waves  (and 
hence  Press-Ewing  waves)  of  about  0.3  sec.  period.  Detailed  inter¬ 
pretation  of  the  interaction  requires  that  all  the  individual  phases 
from  a  variety  of  surface  and  buried  charges  be  analysed  by  the  methods 
of  Appendix  E.  For  this  purpose  it  is  necessary  to  digitize  the  com¬ 
plete  records,  plot  the  particle  trajectories,  and  then  select  regions 
of  each  phase  which  are  free  of  interference  by  the  adjacent  phases. 

The  digital  data  must  then  be  corrected  to  allow  a  ’smooth1  start  and 
end  to  the  phase,  before  proceeding  with  the  dispersion  analysis. 

Since  there  are  several  steps  which  require  subjective  judgment,  it 
is  essential  that  the  analysis  be  repeated  on  several  seisnograms  from 
different  epicentral  distances,  and  from  a  variety  of  charges.  The 
mechanics  of  the  operation  have  been  reduced  to  the  technical  level, 
but,  due  to  the  very  considerable  labour  involved,  it  will  be  some 
time  before  sufficient  analyses  of  individual  phases  are  available 
for  a  scientific  appraisal  of  the  results.  The  report  by  Hastrup 
(1963b)  may  be  taken  as  the  first  step  in  this  direction  for  a  single 
buried  charge.  A  similar  analysis  has  been  undertaken  on  the  immediate 
post-airblast  phase  of  seismograms  60W,  6lVJ  and  62W.  Details  are  given 
in  Appendix  F,  as  a  typical  example  of  the  process  of  reduction. 

The  fact  that  the  post-airblast  record  is  strongly  influenced 
by  the  phase  relationship  of  the  airblast  and  the  pre-airblast  seismic 
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wave  is  illustrated  by  three  seismograms  at  different  azimuths,  but 
the  same  range,  from  the  200,000  lb  charge.  In  seismograms  67V/  and 
71W  the  airblast  is  associated  with  the  crest  of  the  pre-airblast 
wave,  and  the  immediate  post-airblast  record  in  both  cases  is 
strikingly  similar,  and  of  large  amplitude.  Seismogram  74W,  in 
contrast,  is  from  the  same  distance  but  a  radically  different  azi¬ 
muth.  In  this  case  the  airblast  is  associated  with  a  point  well 
below  the  crest  of  the  pre-airblast  wave.  The  post-airblast  region 
is  entirely  different  in  form  from  the  common  pattern  of  67W  and 
71W,  and  is  of  comparatively  low  amplitude.  The  correspondence 
between  these  effects  and  the  effects  shown  in  the  T synthetic 
seismograms1  of  Fig.  22  leaves  little  doubt  that  the  explanation 
given  is  adequate.  An  air-coupled  wave  has  been  superimposed  upon 
the  pre-airblast  strong  motion.  The  classical  form  of  constant 
period  train  in  the  post-airblast  region  is  hardly  to  be  expected* 

An  additional  illustration  of  the  effect  may  be  seen  in  a  comparison 
of  the  40,000  lb  seismograms  from  Watching  Hill,  which  show  marked 
evidence  of  the  influence  of  the  airblast,  and  the  records  from  the 
40,000  lb  charge  on  Drowning  Ford,  a  test  site  in  which  the  evidence 
of  airblast- coupling  is  less  marked.  In  both  these  pairs  of  records 
the  interval  between  the  P  wave  arrival  and  the  airblast  arrival  is 
constant,  indicating  that  the  phase  change  is  not  due  to  variation 
in  the  wind,  but  to  slight  differences  in  the  dispersion  curves. 

It  was  noted  by  Press  (1963)  that  reversal  of  polarity 
occurs  between  the  main  shock  surface  wave  and  the  collapse  shock 
surface  waves.  It  is  possible  that  we  should  look  to  crater  re-bound 
to  explain  some  of  the  features  of  the  multiple  phases  in  both  the 
x  see  also  Smith,  S.W.,  J.G.R.  Vol.  68,  No.  5,  1963. 


vt  ....  1 


-  74  - 


buried  and  surface  charges. 

3.61  It  is  possible  to  follow  the  close  association  of  the  airblast 
and  the  ground  shock  into  the  very  high  pressure  regions,  but,  at  the 
present  time,  it  is  not  possible  to  do  so  with  instrumentation  com¬ 
parable  to  that  used  at  more  moderate  pressure  levels.  At  Suffield 
comparisons  have  been  made  between  the  arrival  time  of  the  air  blast 
and  the  ground  shock  along  lines  reaching  to  within  a  few  feet  of  the 
charge  itself.  The  arrival  time  of  the  airblast  was  measured  by 
foil-contactor-type  switches  and  by  piezo-electric  gauges  mounted  just 
above  the  surface.  The  ground  shock  arrival  was  detected  by  ground 
shock  switches  in  which  a  seismic  mass  actuated  the  closure.  These 
ground  shock  switches  respond  to  radial  acceleration  by  design,  but 
may  have  some  cross-axial  sensitivity.  It  is  not  at  present  certain 
what  the  time  resolution  or  the  threshold  acceleration  is  for  actu¬ 
ation.  The  switches  have  almost  ’dead  beat’  damping,  and  a  natural 
period  of  between  1  and  3  cps  .  Table  7  gives  the  travel  time  data 
obtained  for  the  40,000  lb  1963  trial  (F.E.  556).  The  data  are  also 
plotted  in  the  foxm  of  a  travel-time  curve  in  Fig.  31 « 

It  is  clear  that  there  is  very  close  association  between 
the  airblast  and  the  ground  shock  right  in  to  the  crater  area,  with 
the  two  effects  decelerating  according  to  very  similar  curves.  In 
the  very  high  pressure  region  the  airblast  leads  the  ground  shock,  as 
indicated  by  the  ground  shock  switch.  It  is  not  certain  that  the  switch 
detects  the  first  onset  of  the  ground  shock,  but  certainly  the  first 

*  A  report  on  the  design  of  these  gauges  (and  additional  results)  is 
in  preparation  but  no  publication  date  can  be  given  (S.T.N.  No.  63). 
The  development  is  being  undertaken  by  the  Instrumentation  Group  of 
the  Physics  and  Meteorology  Section,  SES. 
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ground  shock  with  sufficient  amplitude  to  operate  the  switch  follows 
very  closely  upon  the  airblast.  At  intermediate  pressure  levels,  the 
two  effects  become  virtually  synchronous.  A  stage  is  finally  reached 
where  the  airblast  is  lagging  a  ground  wave  of  sufficient  amplitude 
to  trigger  the  switch,  and  thereafter  the  travel  time  curve  for  the 
ground  shock  is  effectively  linear,  while  the  airblast  continues  to 
decay  in  velocity.  This  may  be  compared  with  the  events  recorded 
by  the  Sprengnethers  at  close  ranges.  At  the  shortest  ranges  the 
Sprengnethers  indicate  that  the  airblast  is  synchronous  with  the 
large  amplitude  event,  but  as  one  moves  away  from  GZ  a  moderately 
large  amplitude  dispersed  train  gains  ground  on  the  airblast,  even 
though  the  airblast  is  still  associated  with  the  largest  amplitude. 

In  nuclear  effects  studies  the  early  region  is  sometimes  divided  into 
the  ’  superseismic’  region  where  the  airblast  precedes  the  seismic 
event,  and  the  subsequent  sub-seismic  region.  The  ground  shock 
’precursor  wave’  may  explain  the  ’precursor  air  shock’  occasionally 
seen  in  advance  of  the  main  air  shock. 

3.62  During  several  of  the  large  trials  at  Suffield,  attempts  were 
made  to  record  the  electromagnetic  radiation  associated  with  the 
detonation  (Dobbie  and  Hamilton,  1963) .  No  definite  success  was 
achieved  in  this  phase  of  the  work,  but  consequent  upon  the  attempt, 
it  was  considered  worth  investigating  the  possibility  that  the  electri¬ 
cal  potential  between  various  parts  of  the  test  site  could  change  in 
association  with  the  detonation.  Evidence  from  the  Nevada  Test  Site 
indicated  that  potential  changes  synchronous  with  a  surface  wave  could 
occur.  No  evidence  is  known  to  the  writer  of  similar  changes  in  associ¬ 
ation  with  H.E.  detonations,  and  the  effect  was  therefore  looked  for  in 
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the  1963  multi- ton  detonations . 

Potential  probes  were  installed  at  a  depth  of  20  ft,  at  a 
spacing  of  200  ft,  in  pairs  radial  to  GZ  and  tangential  to  GZ,  at 
4000  ft  from  GZ.  The  intersection  of  the  two  probe  lines  coincided 
with  the  position  of  a  Sprengnether  seismograph. 

The  data  available  up  to  the  present  time  are  not  conclusive. 

No  short  period  variation  was  detected  which  could  be  associated  with 
the  surface  wave  periods.  However,  a  definite  change  in  the  potential 
difference  between  the  probes  was  detected  synchronous  with  the  arrival 
of  the  airblast.  This  was  a  long  period  effect  (order  of  seconds), 
and  there  is  no  conclusion  possible  at  this  stage.  At  the  present  time 
the  effect  appears  to  be  real,  and  cannot  be  traced  down  to  the  recording 
equipment  itself.  Fig.  32  illustrates  the  type  of  record  obtained  with¬ 
out  further  comment*. 

A  second  attempt  to  record  this  effect,  using  a  10,000  lb 
charge,  in  August,  1963*  did  not  confirm  the  existence  of  potential 
change.  Thus  the  effect  may  be  spurious,  or  may  be  associated  with 
large  charges  or  specific  soil  conditions.  Definitive  experiments  will 
be  carried  out  on  the  1963  Bell  Telephone  40,000  lb  trial  and  the 
1,000,000  lb  trial  in  1964. 

•  ••••••oo*ooooo0ooooo«ooo**o*«o«*«ooo*oooe9ooooo*ooooooooooae«*o«o*« 

NOTE.  The  BTL  40,000  trial  was  carried  out  15  August,  1963.  A  record 
similar  to  that  from  the  previous  40,000  lb  trial  was  obtained,  thus 
justifying  further  effort  to  investigate  the  phenomenon. 


The  writer  records  his  appreciation  of  the  effort  being  made  by 
the  Physics  and  Meteorology  Section’s  Instrumentation  Group  in 
this  work. 
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3.70  The  Transverse  Motion 

One  of  the  dominant  features  of  the  seismograms  from  the 
Suffield  test  sites  is  the  relatively  large  amplitude  transverse 
motion.  Such  transverse  motion  is  almost  automatically  dubbed 
"Love  Wave  Motion" .  Surprisingly  few  data  appear  to  be  available 
on  transverse  motion  for  short  periods,  or  for  explosive  sources. 

Virtually  the  whole  of  the  EJP  discussion,  and  the  references  cited 
by  these  authors,  deal  with  continental  and  oceanic  Love  waves  with 
periods  greater  than  about  eight  seconds.  An  exception  is  the  Lg 
phase,  a  short  period  transverse  motion  coupled  to  a  vertical  motion, 
but,  in  this  case  also,  the  discussion  is  reserved  for  the  continental 
phase  with  periods  in  the  range  1-6  seconds. 

Love’ s  theory  treats  the  transverse  motion  as  being  strictly 
independent  of  motion  in  the  other  two  components.  Numerical  compu¬ 
tations  have  been  undertaken  of  Love-type  motion  by  Harkrider  and 
Anderson  (1962)  and  others  for  specific  cases  relevant  to  long  period 
waves.  Horizontal  anisotropy  in  the  medium  will  lead  to  coupling 
between  the  body  wave  modes,  removing  the  division  into  P  and  S  type 
waves.  Presumably  a  similar  coupling  will  occur  in  the  case  of  surface 
waves  but  no  satisfactory  theoretical  analysis  appears  to  be  available. 

In  this  section,  attention  is  concentrated  on  the  experimental 
observations  by  the  writer.  These  data  are  more  extensive  than  any  hitherto 
available  for  the  range  of  periods  below  one  second,  for  a  range  of  charge 
sizes.  Strong  motion  studies  of  the  transverse  wave  for  a  limited  charge 
size  have  been  undertaken  by  the  Saint  Louis  explosion  seismology  group 
(Kisslinger  1959,  1961,  1962  and  Kisslinger  et  al  1961).  Discussion  of 
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some  of  the  Suffield  results  have  been  given  by  Jones  (1962b,  c) . 

We  consider  first  the  evidence  for  repeatability  of  the  wave 
form  in  a  given  azimuth  from  repeated  charges  on  the  same  GZ.  This  is 
a  matter  of  some  importance  as  several  of  the  theories  dealing  with  the 
origin  of  transverse  motion  (see,  for  example,  Press  and  Archambeau, 
1962)  suggested  that  we  must  look  to  the  crater  region  itself.  These 
theories  cannot  account  for  a  consistent  radiation  pattern  on  a  given 
ground,  and  must  fail  if  there  is  evidence  of  such  consistency.  The 
F.E.  536  series  of  10,000  lb  shots  provides  a  direct  check  on  this 
point.  Seismograms  47W,  52!W  and  5$W  were  all  obtained  at  a  distance 
of  2000  ft  from  10,000  lb  charges  detonated  on  the  same  GZ,  the 
craters  being  excavated,  filled,  and  partially  re-consolidated  between 
shots.  The  transverse  surface  wave  portion  of  these  records  is  re¬ 
plotted  in  Fig.  33(a)  on  a  common  time  scale.  Virtually  exact  dupli¬ 
cation  of  the  wave  form  is  evident,  both  in  amplitude  and  frequency. 

Two  records  were  also  obtained  on  another  common  bearing,  but  unfor¬ 
tunately  one  of  these  records  is  the  only  known  instance  of  a  recorder 
being  run  with  batteries  too  weak  to  bring  the  time  mark  motor  up  to 
speed.  However,  it  is  possible  to  match  two  recognizable  peaks  and  so 
correct  the  time  scale  in  this  case.  The  corrected  and  uncorrected 
record,  together  with  the  record  from  the  other  trial,  are  shown  in 
Fig.  33(b).  It  is  quite  evident  that  there  is  virtually  complete 
matching.  Even  one  case  of  matching  is  sufficient  to  prove  the  point, 
without  labouring  the  issue.  Kisslinger  et  al  (1961)  have  also  detec¬ 
ted  a  similar  repeatability  on  another  test  site,  so  that  the  ’radiation 
pattern’  appears  to  have  a  distinct  geographical  orientation,  presumably 
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based  upon  sub-surface  features.  A  similar  repeatability  is  evident 
in  the  records  from  the  Nevada  Test  Site  (Seismograms  $0  -  83N),  where 
there  is  no  air-blast  effect. 

Kisslinger  et  al,  using  very  small  charges  which  resulted  in 
short  and  comparatively  simple  seismograms,  were  able  to  demonstrate  a 
lobed  radiation  pattern  similar  to  that  predicted  on  the  Knoppoff  and 
Gilbert  Model  (Knoppoff  and  Gilbert,  I960) .  Nodal  lines,  of  the  type 
demonstrated  by  Kisslinger  et  al,  must  exist  for  any  transverse  motion 
which  is  not  a  pure  rotation.  Nevertheless,  the  present  writer  is  not 
convinced  that  these  nodal  lines  are  permanent  in  direction.  The  few 
azimuthal  variation  data  available  from  the  present  work  are  re-plotted 
in  Fig.  34  in  the  form  of  a  time- displacement  plot  on  seven  bearings.  It 
will  be  seen  that,  while  there  is  certainly  a  tendency  to  larger  ampli¬ 
tudes  on  the  25°  line,  there  is  no  real,  evidence  of  a  nodal  line.  In 
the  experience  of  the  writer,  the  overall  pattern  of  amplitude  in  the 
transverse  trace  is  quite  independent  of  the  bearing  of  the  seismic 
line.  The  shape  of  the  profile  repeats  along  a  given  line,  and  varies 
from  azimuth  to  azimuth,  but  the  duration  of  the  train,  and  the  mean 
amplitude,  appears  to  be  essentially  independent  of  the  particular  azi¬ 
muth  selected.  Dogmatic  statements  based  upon  the  few  data  available 
on  azimuthal  variation  on  the  Suffield  site  would  be  quite  out  of  order. 
It  is  proposed  that  a  very  detailed  study  be  undertaken  on  the  Watching 
Hill  test  site  with  moderate  size  charges  (about  100-300  lb)  to  investi¬ 
gate  the  radiation  pattern  at  intervals  of  not  more  than  five  degrees. 
With  the  limited  instrumentation  available,  this  will  require  a  minimum 
of  100  shots.  Due  to  the  cost  of  such  trials,  it  is  essential  that  the 
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programme  be  linked  with  other  programmes  to  make  full  use  of  the  charges, 
and,  until  such  a  coordinated  programme  can  be  developed,  no  time 
schedule  can  be  given.  An  alternative  method  would  be  to  use  a  simpler 
fom  of  geophone  together  with  existing  tape  facilities.  This  would 
reduce  the  required  number  of  shots  to  about  ten,  but  would  require  the 
rental  or  purchase  of  new  equipment.  There  are,  however,  distinct 
advantages  in  using  mechanical  seianographs  for  this  type  of  study, 
since  this  avoids  the  necessity  of  analogue  studies  to  determine  the 
true  response  of  the  system. 

Attention  has  been  drawn  (Jones,  1962c)  to  the  apparent  coup¬ 
ling  between  the  transverse  waves  and  the  radial-vertical  waves  on  the 
Suffield  test  sites.  With  buried  charges  there  is  an  association  of 
Love  and  Rayleigh  type  events,  but  this  association  may  not  amount 
to  coupling  so  much  as  a  similarity  in  the  dispersion  curves.  The 
coupling  becomes  evident  in  the  post-air-blast  region  of  the  seismo¬ 
grams  from  surface-burst  charges.  Attention  is  drawn  to  seismograms 
65W  and  71W  as  examples  of  this  effect,  which  does  not  appear  to  be 
invariable.  It  should  also  be  noted  that  even  where  direct  coupling 
is  not  evident,  as  in  the  records  from  the  10,000  lb  charges  on  the 
Watching  Hill  test  site,  the  Radial-Vertical  wave  onset  and  the  Love 
wave  onset  are  invariably  associated,  independent  of  azimuth.  It 
seems  unlikely  that  such  an  association  can  be  explained  by  a  mechan¬ 
ism  of  formation  which  depends  upon  either  a  pre-existing  or  a  develop¬ 
ing  crack  structure  in  the  ground.  It  would  appear  that  some  coherent 
flexure  of  the  surface  layers  is  required. 

The  initial  analyses  by  Hastrup  (1963b)  of  the  two  main 
phases  from  the  buried  charges  appear  to  indicate  that  the  transverse 
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phase  associated  with  the  dominant  prograde  mode  in  the  radial-vertic al 
components  may  be  closely  coupled  to  each  other,  since  both  the  spectra 
and  the  dispersion  curves  are  almost  precisely  similar  to  each  other. 

It  would  appear  that  this  does  not  apply  to  the  lower  velocity  phases, 
but  considerable  routine  analysis  is  required  before  definitive  con¬ 
clusions  can  be  drawn. 
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3.30  Remote  Effects 

The  remote  seismic  effects  produced  by  these  explosions  are, 
generally  speaking,  beyond  the  scope  of  this  dissertation.  One  aspect 
of  remote  action,  however,  is  of  interest  due  to  the  muisance  value. 
Local  areas,  comparatively  remote  from  the  charge,  are  occasionally 
T  shaken*  by  the  explosions,  with  or  without  actual  damage  to  struc¬ 
tures.  This  phenomenon  is  essentially  due  to  atmospheric  focussing 
of  the  blast  wave  back  to  the  surface,  but  it  is,  in  the  opinion  of 
the  writer,  likely  that  the  observed  effects  are  seismic  in  immediate 
origin.  The  essentially  acoustic  wave  is  focussed  on  to  a  limited 
area  and  produces  an  air- induced  wave,  possibly  described  more 
closely  by  the  Bateman  theory  than  the  Press-Ewing  theory. 

The  coupled  wave  observed  by  the  vrriter  at  55^000  ft  from 
the  first  40,000  lb  charge  (para.  2.43)  was  almost  certainly  of 
Press-Ewing  type.  However,  the  first  example  of  interest  to  the 
present  discussion  involved  a  10,000  lb  charge  on  the  Drowning  Ford 
site.  Immediately  after  the  explosion  a  report  was  received  from 
a  Mr  Walker,  who  farms  land  East  of  the  river,  some  fifteen  miles 
from  the  test  site.  He  claimed  that  his  house  and  other  buildings 
were  badly  shaken,  with  some  cracking  of  the  wall  panels.  As  a 
second  10,000  lb  trail  was  due  a  week  later,  the  writer  and  Mr  B.J. 
Perry,  Superintendent  of  Research  at  SES,  arranged  to  be  present  at 
TWalkerfs  FarmT  .  An  examination  was  made  of  the  structure  of  the 
house,  and  the  alleged  f cracking1  investigated.  In  the  opinion  of 
the  writer  these  cracks  were  not  fresh,  but  may  have  been  made  more 
obvious  by  dust  being  shaken  out  of  then.  Mr  Walker  claimed  that 


. 
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previous  trials  had  not  bothered  the  family,  though  they  were  always 
aware  when  trials  were  being  conducted.  The  SES  representatives, 
expecting  to  find  the  case  wildly  exaggerated,  decided  to  stand  out¬ 
side  the  house  and  observe  the  arrival  from  there.  The  radio  count¬ 
down  was  monitored,  and  a  stopwatch  started  at  zero .  The  cloud  was 
observed  within  seconds,  and  at  the  predicted  arrival  time  (assuming 
sonic  travel  for  this  estimate)  the  observers  were  shaken  out  of 
their  complacency  by  a  loud  ’thunderclap1  and  pronounced  shaking  of 
the  ground.  There  were  distinct  sound  effects  from  within  the  build¬ 
ing,  but  no  additional  damage  was  detected.  Mr  Walker  assured  the 
observers  that  the  effects  on  this  occasion  were  relatively  mild. 

Since  this  time  great  care  has  been  taken  to  ensure  that  no  large 
charges  are  functioned  under  conditions  when  focussing  is  predicted, 
and  in  general  no  complaints  have  been  received,  with  two  exceptions. 

The  first  exception  involved  the  Gust  farm,  some  fourteen  miles  from 
the  test  site,  where  a  reaction  similar  to  the  earlier  one  at  Walker’s 
farm  was  obtained.  The  second  exception  occurred  when  a  10,000  lb 
charge  was  fired  at  a  point  somewhat  nearer  the  village  of  Ralston. 

In  this  case,  minor  damage  (a  broken  plate  glass  window)  was  reported 
from  Medicine  Hat,  about  30  miles  away.  It  is  interesting  that  in  this 
case  people  in  the  street  did  not  hear  the  explosion,  while  those  within 
the  building  assumed  at  first  that  a  child  had  broken  the  window. 

Whether  or  not  this  damage  is  attributable  to  the  blast  is  a  moot  point. 

It  is  very  difficult  to  obtain  reliable  information  on  these 
effects  since  trials  are  carried  out  under  conditions  which  minimize 
focussing.  Attempts  have  been  made  to  compare  normal  zones  with  zones 
of  enhanced  pressure,  but  even  with  the  aid  of  a  helicopter,  it  is 


-  84  - 


difficult  to  arrange  for  instrumentation  to  be  at  the  correct  locations. 
From  the  limited  data  available,  it  appears  that  pressures  of  up  to 
0.01  psi  are  detectable  as  a  rumble  like  distant  thunder,  but  there  is 
no  detectable  ground  roll.  When  the  pressure  level  approaches  0.03 
the  audible  effect  is  like  rolling  thunder  and  there  is  some  sensation 
of  the  ground  shaking.  Presumably,  effects  such  as  those  at  Walker’s 
Farm  correspond  to  slightly  higher  pressure  levels,  possibly  as  high 
as  0.1  psi.  This  pressure  level  in  itself  is  too  low  to  cause  damage, 
but  may  cause  sufficient  ’ air- induced’  vibration  to  cause  very  minor 
structural  damage.  It  would  be  extremely  helpful  if  a  serious  attack 
could  be  made  upon  this  problem,  preferably  by  studying  smaller  charges 
detonated  under  conditions  favouring  focussing  upon  uninhabited  areas. 

At  the  present  time  it  is  very  doubtful  that  sufficient  data  are  avail¬ 
able  to  verify  the  computer  programmes  used  for  predicting  anomalous 
blast  effects.5*' 


The  investigation  of  these  effects  has  been  undertaken  by  Mr.  J.J. 
Vesso  (in  the  field)  under  the  direction  of  predictions  made  by 
Messrs  Johnson,  Gilbert  and  Clink.  At  the  present  time  data,  are 
being  collected  and  no  final  report  is  within  sight .  Attention 
is  drawn  to  Johnson  et  al  (1962)  which  gives  some  information  on 
the  methods. 
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4.00  CONSOLIDATION  OF  THE  OBSERVATIONS 

4.10  The  experimental  data  given  in  this  dissertation  will  require 
several  man-years  of  routine  technical  reduction  before  the  fine 
details  have  been  elucidated.  It  is  obviously  necessary,  on 
scientific  grounds,  to  carry  the  experimental  approach  well  into 
the  plastic  zone  in  order  to  determine  the  mechanism  by  which  the 
various  modes  are  excited.  Similarly,  some  limited  work  is  required 
at  longer  ranges  where  the  modes  are  well  separated.  The  eventual 
ramifications  of  the  programme  will  depend  in  part  upon  these  scien¬ 
tific  considerations,  and  in  part  upon  economic  and  political  decisions 
regarding  seismological  research  within  the  Defence  Research  Board. 

It  is  therefore  not  irrelevant  to  quote  Wilson  and  Caless  (1963) * 

"Determination  of  the  precise  way  in  which  signals  are 
generated  by  a  seismic  source  is  one  of  the  least  explored  research 
areas  of  seismology" . 

The  proceedings  of  the  Conference  of  Experts  in  Geneva 
indicate  quite  clearly  the  need  for  both  experimental  and  theoretical 
investigations  in  this  vital  area.  At  the  beginning  of  this  research, 
little,  if  any,  knowledge  existed  of  the  form  to  be  expected  in  the 
seismograms  from  the  strong  motion  region  close  to  high  yield  weapons 
or  other  explosive  devices  in  the  region  of  the  Western  Plains,  an 
area  which  is  not  small  or  unimportant  in  the  field  of  continental 
defence.  While  much  of  the  fine  detail  remains  obscure,  the  data 
in  the  present  dissertation  are  such  that  a  clear  picture  has  now 
emerged  of  the  general  characteristics  of  the  ground-air  interaction 
and  the  sequence  of  events  to  be  expected  in  the  strong  motion  region. 
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The  importance  of  this  as  a  source  of  input  data  in  structural  effects 
work  cannot  be  over-estimated. 

4.20  We  may  consolidate  what  we  now  know  regarding  the  strong  motion 
effects  of  large  explosions  in  the  plains  area  into  a  series  of  state¬ 
ment  s,  which  include  the  following: 

(1)  The  P  and  S  waves,  at  ranges  where  they  have  separated  from  the 
surface  waves,  may  be  neglected  completely  as  input  parameters  in  structural 
damage  studies.  The  amplitudes  are  far  too  small  to  be  significant  compared 
with  either  the  direct  air  blast  or  the  surface  waves. 

(2)  The  air  blast  and  the  seismic  waves  are  closely  associated  at  all 
distances  within  the  strong  motion  region,  which  we  may  define  for  this 
purpose  as  the  circular  area  centred  on  Ground  Zero  and  extending  to  an 
airblast  overpressure  of  about  one-tenth  psi.  This  definition  does  not 
contravene  the  more  usual  seismic  definition  of  the  strong  motion  region 

as  being  that  in  which  instruments  with  gains  not  greater  than  one  hundred 
produce  usable  records. 

(3)  Within  the  plastic  zone,  the  air  shock  and  the  ground  shock 
velocities  decay  according  to  similar  curves,  which  indicates  that  the 
ground  shock  is  generated  essentially  by  the  passage  of  the  airblast 
over  the  surface.  The  transition  between  this  situation  and  the  situ¬ 
ation  existing  in  the  case  of  contained  charges  requires  clarification, 

(4)  At  a  certain  distance,  which  corresponds  with  about  1.25  psi 
peak  airblast  pressure  in  the  case  of  a  40,000  lb  charge,  the  ground 
shock  and  air  shock  curves  cross  each  other,  so  that  the  ground  shock 
front  leads  the  airblast.  In  this  region  the  ground  shock  is  of  approxi¬ 
mately  constant  velocity,  while  the  air  shock  continues  to  decay  in 


. 
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velocity.  The  ground  shock  becomes  a  ’precursor  wave1  to  the  events 
associated  with  the  airblast.  The  corresponding  pressure  level  from 
the  200,000  lb  charge  was  monitored  by  a  Sprengnether  seismograph  and 
shows  a  similar  pattern  of  a  short  precursor  wave  immediately  ahead 
of  the  airblast. 

(5)  The  maximum  amplitude  in  all  seismograms  occurs  in  the  region 
immediately  following  the  arrival  of  the  airblast. 

(6)  Surface  waves  generated  at  the  Suf field  test  sites  (and  else¬ 
where  according  to  the  limited  data  of  Kisslinger,  Kuo  and  Howell)  ex¬ 
hibit  both  retrograde  (Rayleigh  type)  and  prograde  (Sezawa  type)  orbits. 
The  common  pattern  at  Suffield  consists  of  a  large  amplitude  prograde- 
retrograde  pair,  preceded  and  followed  by  low-amplitude  retrograde 
phases.  In  many  cases,  however,  both  the  large  amplitude  phases  exhibit 
prograde  orbits. 

(7)  Prograde  and  retrograde  orbits  are  normally  orthogonal  to 
each  other.  In  certain  cases,  however,  both  prograde- retrograde  and 
retrograde-prograde  conversions  may  occur  with  the  tilt  of  the  orbit 
remaining  constant. 

(S)  The  prograde  orbit  is  normally  tilted  forward  of  the  vertical 
and  the  retrograde  phase  tilted  backward. 

(9)  The  pre-airblast  phase  is  retrograde  at  points  close  to  the 
epicentre.  This  phase  converts  into  a  prograde  phase  at  greater 
distances,  through  a  transition  zone  where  the  orbit  is  virtually  linear 
and  tilted  from  the  vertical. 

(10)  Both  buried  and  surface-burst  charges  excite  similar  sequences 
of  phases  showing  prograde -retrograde  succession,  more  than  one  such 
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pair  being  evident. 

(11)  Where  a  phase  continues  for  several  cycles  beyond  a  transition 
zone,  the  tilted  elliptical  orbit  tends  to  become  more  circular  and 
vertical,  indicating  that  the  pure  form  is  a  vertical  ellipse  with  a 
small  ratio  of  minor  to  major  axes. 

(12)  Both  retrograde  and  prograde  phases  exist  which  exhibit 
dispersion  curves  which  cross  the  sonic  (and  blast)  velocity,  leading 
to  the  possibility  (and  theoretical  necessity)  of  Press-Ewing  type 
coupling  between  the  airblast  and  the  surface  waves. 

(13)  At  points  remote  from  the  charge,  the  dominant  seismic  effect 
(and  the  only  one  detectable  by  low  gain  instruments)  is  a  Press-Ewing 
type  air- coupled  wave. 

(14)  At  points  comparatively  close  to  GZ  there  is  evidence  that  a 
Press-Ewing  wave  is  superimposed  upon  a  surface  wave  excited  by  a 
different  mechanism. 

(15)  The  precise  form  of  the  post  airblast  seismic  wave  form  in 
the  region  where  action  (14)  is  occurring  depends  sharply  upon  the 
relationship  between  the  airblast  arrival  and  the  phase  of  the  surface 
wave  at  the  moment  of  arrival.  Very  slight  changes  in  this  relation¬ 
ship  can  cause  radical  changes  in  the  wave  form  of  the  seismogram. 

(16)  In  the  strong  motion  region  the  dominant  seismic  effect  is 
apparently  closely  coupled  to  the  airblast.  It  is  probable  that  this 
coupling  is  a  direct  induction  of  the  surface  wave  by  a  TLambf  or 
TMilesT  type  mechanism,  but  at  the  present  stage  of  the  theory  this 
cannot  be  verified  for  a  multi-layered  medium  supporting  a  multiplicity 
of  phases. 


. 
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(17)  The  amplitude  of  the  post-airblast  surface  disturbance  does 
not  of  necessity  decay  with  distance,  but  may  actually  increase,  due 
partly  to  the  combining  of  independent  wave  groups  in  varying  phase 
relationship,  and  partly  due  to  the  rotation  of  the  tilted  orbits 
towards  the  vertical. 

(l$)  In  the  seismograms  from  contained  charges,  the  maximum  ampli¬ 
tude  in  the  vertical  trace  from  the  retrograde  phase  decreases  with 
distance  fairly  uniformly,  but  there  is  a  tendency  for  the  amplitude 
of  the  prograde  phase  to  increase.  The  length  of  the  major  axes, 
however,  appears  to  decrease  with  distance  in  both  phases,  and  it  is 
suggested  that  amplitude  scaling  should  always  be  performed  upon  the 
axes  of  the  orbit  rather  than  upon  a  single  component  of  the  motion. 

(19)  The  theoretical  dispersion  curves  for  models  which  approximate 
to  the  test  site  conditions  all  show  a  common  pattern  in  which  Press- 
Ewing  coupling  is  predicted  for  the  surface  waves  from  moderately  large 
charges,  but  is  not  expected  as  a  dominant  effect  for  small  charges. 

(20)  The  periods  actually  excited  by  the  range  of  charges  used 

up  to  the  present  fall  within  the  allowable  range  for  the  Rayleigh  mode 
computed  by  the  Anderson  programme. 

(21)  The  peak  frequencies  in  the  immediate  post-airblast  phase  of 
the  seismograms  correspond  very  closely  to  the  two  frequencies  associ¬ 
ated  with  Press-Ewing  coupling  in  the  dispersion  curve  for  a  theoreti¬ 
cal  model  similar  to  the  test  site  conditions. 

(22)  Only  the  layering  down  to  and  including  the  first  Upper 
Cretaceous  layer  is  significant  for  charges  of  moderate  size  (up  to 
the  order  of  105  lb  TNT).  It  is  probable  that  explosions  in  the 
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kiloton  (and  larger)  range  will  excite  long  period  waves  which  will 
depend  upon  the  deeper  layers. 

(23)  Very  slight  changes  in  the  upper  layers  of  the  system  can 
produce  very  large  changes  in  the  seismic  effects.  Thus  the  system  is 
sensitive  to  changes  in  the  weathered  layer. 

(24)  In  all  components,  the  dominant  frequency  is  sharply  dependent 
upon  the  charge  size.  The  larger  the  charge  the  longer  the  period  of 
the  wave  excited  as  a  measurable  strong  motion  effect. 

(25)  Small  charges,  of  the  order  of  10  lb  TOT  excite  low  velocity 
waves  which  show  a  dominant  period  of  about  0.2  sec.  Larger  charges  do 
not  excite  this  period  as  a  dominant  one,  but  tend  to  excite  waves  in 
the  period  range  0.3  sec.  to  1  sec.  The  precise  range  depends  upon  the 
charge  size. 

(26)  Due  to  the  relatively  long  period  waves  excited  by  the  larger 
charges,  the  commonly  available  instrumentation  is  not  well  suited  to 
the  study  of  the  surface  waves.  The  one  second  period  region  corres¬ 
ponds  to  a  transition  zone  in  the  instrumentation  which  is  commercially 
available  at  the  present  time.  This  range  is  of  importance  in  the  study 
of  the  surface  waves  from  tactical  weapons. 

(27)  At  very  long  ranges  the  dominant  seismic  effect  appears  to  be 
aircoupling,  and  it  is  probable  that  the  effectiveness  of  anomalous  blast 
propagation  in  causing  structural  damage  depends  upon  how  efficiently  the 
focussed  wave  excites  an  air-coupled  wave.  Attention  is  drawn  to  the 
recent  paper  by  Frantti  (1963)  which  shows  the  correspondence  between  the 
spectra  of  the  microseismic  background  and  the  air-coupled  wave  at  100  km 

from  the  Suf field  site. 

* 


"The  Nature  of  High  Frequency  Noise  Spectra",  Geophysics,  Vol.  2$, 
No.  k,  PP  547-562. 
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(2S)  The  transverse  motion  in  a  given  azimuth  is  repeated  very 
closely  for  repeated  shots  on  the  same  Ground  Zero,  indicating  that 
the  radiation  pattern  is  governed  by  sub-surface  features  of  greater 
scale  than  the  crater  zone. 

(29)  There  is  strong  evidence  of  coupling  between  the  Rayleigh 
and  Love  type  waves  in  the  immediate  post-air-blast  region  of  the 
seismograms  from  surface-burst  charges.  In  the  seismograms  from 
contained  charges  there  is  an  association  between  the  two  wave  types, 
but  the  coupling  is  not  tight.  There  is  some  evidence  from  the  spectra 
and  dispersion  curves  that  the  prograde  radial-vertical  mode  and  its 
associated  transverse  mode  are  coupled.  The  evidence  is  against  coup¬ 
ling  between  the  retrograde  mode  and  its  associated  transverse  mode. 

The  last  two  statements  are  very  tentative,  being  based  solely  upon  the 
analysis  by  Hastrup  (1963b)  of  the  1000  lb  buried  charge  records. 

(30)  Dispersion  calculations  in  regions  composed  of  more  than  one 
phase  are  of  doubtful  validity  even  when  obtainable  in  an  ’apparently 
correct’  form.  It  is  a  pre-requisite  of  successful  dispersion  curve 
calculations  that  they  be  performed  on  a  single  phase  which  propagates 
without  spectral  distortion.  This  means,  in  general,  that  the  analysis 
of  strong  motion  seismograms  must  be  associated  with  the  analysis  of 
records  of  the  same  events  at  greater  epicentral  distances.  It  is 
necessary  therefore  to  follow  the  events  at  successive  closely  spaced 
distances  out  to  such  a  distance  that  acceptable  separation  of  the 
phases  has  occurred.  This  entails,  unfortunately,  the  use  of  a  variety 
of  recorders  of  increasing  gain  as  the  epicentral  distance  increases. 

It  is  essential  that  the  seismometer  spacing  be  kept  sufficiently  small 
to  allow  wave  velocity  estimates  to  be  possible,  due  to  the  similarity 
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in  the  wave  forms  at  successive  stations. 

4.30  Experimental  development  of  the  research  into  the  strong  motion 
seismology  of  large  explosions  may  be  expected  to  take  the  course 
outlined  below. 

During  the  twelve  months  following  August,  19&3 >  explosions 
at  Suffield  will  include  about  ten  charges  of  1000  lb  TNT  buried  at  a 
depth  of  50  ft  on  the  Watching  Hill  Site.  These  charges  are  designed 
specifically  as  energy  sources  for  the  SES-U.  of  A.  MOHO  programme, 
but  advantage  will  be  taken  of  them  to  improve  upon  the  results  obtained 
with  the  earlier  1000  lb  charges.  Specifically,  close  arrays  of  seismo¬ 
meters  will  be  installed  at  ranges  of  between  4000  and  6000  ft,  the 
seismometer  spacing  being  varied  between  50  ft  and  500  ft.  In  this  way 
it  is  expected  that  close  definition  of  the  spectra  and  the  dispersion 
curves  of  each  phase  can  be  obtained. 

In  July,  1964,  Suffield  will  be  the  site  of  a  detonation  of 
1,000,000  lb  of  TNT  on  the  surface  of  the  Watching  Hill  range.  During 
this  trial  an  attempt  will  be  made  to  obtain  a  corresponding  set  of 
records  to  that  obtained  from  the  1000  lb  charges.  That  is  to  say, 
no  attempt  will  be  made  to  cover  the  complete  strong  motion  region,  and 
attention  will  be  concentrated  instead  upon  the  limits  of  the  region 
where  the  phases  will  be  moderately  well  separated.  A  minimum  of  four 
seismometers  will  be  used  at  a  spacing  of  200  ft.  One  seismometer  will 
be  located  closer  to  GZ  (at  about  S000  ft)  where  the  effect  of  the 
airblast  will  be  marked. 

These  records,  and  the  existing  records  included  in  this 
dissertation,  will  be  the  subject  of  continuous  routine  analysis 
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according  to  the  methods  of  Appendices  E  and  F,  a  process  which  may  be 
expected  to  take  about  two  years. 

Theoretical  development  should  proceed  by  an  extension  of 
the  Miles  formulation  of  the  blast  wave  problem  to  the  case  of  a 
multilayered  medium.  Due  to  the  complexity  of  the  blast  wave  propa¬ 
gation  and  the  vicinity  of  the  zone  of  plastic  deformation,  it  is  probable 
that  only  qualitative  agreement  may  be  expected  between  the  theoretical 
predictions  and  the  seismograms  included  herein.  It  is  therefore  sug¬ 
gested  that  a  concurrent  attempt  be  made  to  carry  out  a  model  study  on 
the  surface  waves  from  airburst  charges,  which  may  be  expected  to  fit 
the  theory  more  exactly.  Little  model  work  on  surface  waves  has  been 
reported  in  the  literature  up  to  the  present  time.  It  appears,  however, 
that  the  technology  is  now  sufficiently  advanced  to  justify  such  work. 

The  model  should  be  of  comparatively  large  size  (order  ten  feet  square 
and  correspondingly  deep),  and  modern  plastic  materials  should  be  used 
which  offer  the  possibility  of  approximately  correct  dynamical  scaling, 
including  the  scaling  of  shear  velocity  contrasts  between  layers.  The 
use  of  small  airburst  charges  at  about  the  scale  equivalent  of  a  1  lb 
charge  at  10  ft  will  prevent  plastic  deformations  and  allow  the  elastic 
response  to  be  studied  by  repeated  explosions.  In  order  to  obtain 
correct  dynamical  scaling,  it  may  be  necessary  to  depart  from  air  at 
atmospheric  pressure  as  the  medium  for  the  detonation.  This  will  depend 
upon  the  plastic  material  chosen  for  the  upper  layer  of  the  model. 
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5.00 


THE  SEISMOGRAMS 


The  seismograms  in  the  following  Section  have 
been  prepared  by  S.A.  Cyganik  under  the 
direction  of  the  writer.  Preparation  has 
included  separation  of  the  traces  to  avoid 
overlap,  removal  of  the  0.02  second  time  marks 
due  to  the  difficulty  in  reproduction  at  small 
scale,  and  the  addition  of  explanatory  captions. 
The  seismograms  may  be  used  for  all  interpretive 
purposes  except  those  involving  very  low  ampli¬ 
tude  or  high  frequency  events.  The  careful 
reproduction  by  the  SES  Photographic  Section 
and  by  N.  Bonin  of  the  SES  Library  staff  is 
gratefully  acknowledged. 
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TYPICAL  5  TON  SEISMOGRAMS  —  WATCHING  HILL  SITE 
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OPERATION  STAGECOACH 
SEISMIC  RECORDS  — 3000  FT.  FROM  G.Z. 


OPERATION  STAGECOACH 


SEISMIC  RECORDS -4000  FT.  FROM  G.Z. 


NOTE  : 


OPERATION  STAGECOACH 
SEISMIC  RECORDS  —  5000  FT.  FROM  G.Z. 
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SEISMOGRAMS  FROM  record  vertical  scale  (inches) 

SPHERICAL  IOOO  LB.  CHARGE  BURIED  23.5  FT. 

WATCHING  HILL  SITE 
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SEISMOGRAMS  FROM 

CYLINDRICAL  IOOO  LB.  CHARGE  BURIED  50  FT.  —  SOIL  TAMPED 

WATCHING  HILL  SITE 
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SEISMOGRAMS  FROM 

SPHERICAL  1000  LB.  CHARGE  ON  SURFACE 
WATCHING  HILL  SITE 
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TABLE  1 


PREDICTED  PRESSURE-DISTANCE  DATA  FOR  200, 

,000  LB  TNT 

SURFACE  BURST  AT 

THE  WATCHING  HILL 

SITE 

Distance 

Pressure 

Distance 

Pressure 

(ft) 

.  (psi) 

(ft) 

(psi) 

200 

9g . 549 

g500 

.254 

300 

37.302 

9000 

.239 

400 

19.7g4 

9500 

.226 

500 

12.524 

10000 

.214 

600 

g.gl9 

10500 

.203 

700 

6.661 

11000 

.193 

goo 

5.2g3 

11500 

.lg4 

900 

4.342 

12000 

.176 

1000 

3.666 

13000 

.162 

1250 

2.60g 

14000 

.150 

1500 

2.007 

15000 

.140 

1750 

1.624 

16000 

.131 

2000 

1.361 

17000 

.123 

2500 

1.024 

igooo 

.116 

3000 

,g20 

19000 

.109 

3500 

.6g2 

20000 

.104 

4000 

.5g4 

22500 

.092 

4500 

.511 

25000 

.0g2 

5000 

.454 

27500 

.075 

5500 

.40g 

30000 

,06g 

6000 

.370 

40000 

.051 

6500 

.339 

50000 

.041 

7000 

.313 

75000 

.027 

7500 

.290 

100000 

.020 

gooo 

.271 

TABLE  2 

HOPKINSON’S  LAW  SCALING  FACTORS  FOR  VARIOUS 

CHARGES 

Charge 

Scale  to  Scale  to 

10, 

000  lb  scaled 

100  ft  for  Charge 

(lb) 

1  lb  10,000  lb 

to  Charge 

scaled  to  10,000  lb 

g 

2.00  .092g 

10. g 

logo 

30 

3.107  .1442 

6.93 

693 

60 

3.936  .lglg 

5.50 

550 

500 

7.93g  .36g4 

2.71 

271 

10000 

21.54  1.0 

1.0 

100 

40000 

34.20  1.5g7 

.6299 

63 

200000 

5g.4g  2.714 

.3g 

3g 

TABLE  3 

CLASSIFICATION  OF  SOIL  SAMPLES  FROM  BORING  AT  GZ  FOR  1961  100-TQN  TEST  SHOT 


Sample  No. 

Depth 

(ft) 

Group 

Symbol 

Classification  and  Remarks 

Water 

Content 

1 

1 

ML 

tan  clayey  silt,  soft 

7.6 

2 

5 

CL 

tan  silty  clay,  soft 

11.2 

3 

7 

CL 

tan  silty  clay,  med. 

13.3 

4 

10 

CL 

tan  silty  clay,  med. 

14.4 

5 

12 

CL 

tan  silty  clay,  med. 

14.5 

6 

15 

CH 

gray  clay,  stiff 

26.0 

7 

17 

CL 

gray  clay,  med. 

23.9 

8 

20 

CL 

gray  silty  clay,  med. 

21.8 

9 

22 

CL 

gray  clay,  soft 

34.4 

10 

25 

CH 

gray  clay,  med. 

35.8 

11 

27 

CH 

gray  clay,  med. 

33.4 

12 

30 

CH 

gray  clay,  soft 

36.6 

13 

32 

CL 

gray  silty  clay,  soft 

34.1 

14 

35 

CH 

gray  clay,  med. 

38.8 

15 

37 

CH 

gray  clay,  soft 

37.3 

16 

40 

CH 

gray  clay,  soft 

25.1 

17 

42 

CH 

gray  clay,  w/coarse  sand 
(jar  broken 

18 

45 

GC 

gray  gravel,  sand,  clay, 
mixed  and  soft 

15.2 

19 

48.5 

CH 

gray  clay,  w/few  small  gravel, 
med. 

27 .0 

Atterburg  Limits 


Sample 

LL 

PL 

PI 

2 

36.2 

18.0 

18.2 

4 

48.0 

20.1 

27.9 

7 

47.6 

20.7 

26.9 

10 

51.4 

22.1 

29.3 
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TABLE  $ 


AIRBLAST  TRAVEL-TIME  FROM  SEISMOGRAMS 


Seismogram _  Measured  Corrected  to  Scaled  Scaled 

Code  Wt  Dist.  Time  Zero  Wind  to  Distance 

(lb)  (ft)  (sec)  and  32°F  10,000  lb  (ft) 


3D 

g.5 

50 

.025 

.0239 

.25g5 

540 

4D 

11 

75 

.055 

.0526 

.5676 

gio 

3D 

it 

100 

.ogo 

.0771 

.g329 

logo 

6D 

ti 

125 

.105 

.1007 

1.0g76 

1350 

7D 

it 

150 

.130 

.1251 

1.3515 

1620 

10D 

n 

300 

.255 

.2444 

2.640 

3240 

11D 

30 

100 

.Og 

- 

- 

693 

12D 

11 

150 

.13 

.1237 

.g575 

1039 

13D 

it 

200 

.145 

.13g2 

.957g 

13g6 

14D 

11 

250 

.19 

.lgl7 

1.2590 

1732 

15D 

n 

300 

.24 

.2290 

1.5g69 

2079 

16D 

it 

400 

.355 

.3207 

2.2222 

2772 

17D 

it 

500 

.405 

.3gg4 

2.6917 

3465 

lgD 

n 

600 

.4g 

.4601 

3.lgg7 

41 5g 

19D 

it 

700 

.57 

.5459 

3.7g36 

4g51 

20D 

60 

250 

.195 

.lg93 

1.0412 

1375 

21D 

it 

400 

.325 

.3140 

1.7270 

2200 

2gD 

11 

600 

.51 

.4919 

2.7057 

3300 

29D 

n 

goo 

.6g 

.6557 

3.6066 

4400 

34D  10,000 

2000 

1.5g 

- 

- 

2000 

35D 

it 

2500 

2.16 

2.1249 

2.1249 

2500 

36D 

11 

3000 

2.5g 

2.5755 

2.5755 

3000 

37  D 

« 

3500 

3.065 

3.0110 

3.0110 

3500 

3gD 

ti 

4000 

3.4g 

3.4473 

3.4473 

4000 

39D 

it 

4500 

4.00 

3.9519 

3.9519 

4500 

40D 

it 

5450 

4.7g5 

4.762g 

4.762g 

5450 

AID 

11 

6000 

5.52 

5.4415 

5.4415 

6000 

42D 

it 

6500 

6.00 

5.9146 

5.9146 

6500 

43D 

11 

7000 

6.3g 

6.324 

6.324 

7000 

46W 

n 

5000 

4.354 

4.2646 

4.2646 

5000 

52W 

it 

2100 

1.755 

1.7g39 

1.7g39 

2100 

53W 

11 

2000 

1.55 

1.5g05 

1.5805 

2000 

56W 

ti 

2000 

1.66 

1.67g7 

1.67g7 

2000 

57W 

n 

1974 

1.655 

1.6476 

1.6476 

1974 

59W  40,000 

3000 

2.43 

2.303g 

1.4350 

lg90 

60W 

11 

4000 

3.31 

3.1329 

1.9515 

2520 

62W 

it 

7000 

5.95 

5.6390 

3.5125 

4409 

. 

- 

TABLE  6 


DATA  RELATING  TO  P  WAVE  AM)  EVENT  A 

(Vertical  Component) 


.smog  ram 

No. 

Distance 

(ft) 

P  Wave 
Arrival 
(sec) 

Period 

(sec) 

Event  "j 

12D 

150 

.Og 

13D 

200 

.0g5 

1  AD 

250 

.09 

15D 

300 

.12 

.Og 

16D 

400 

.13 

.Og 

17D 

500 

.14 

.Og 

lgD 

600 

.16 

.Og 

19D 

700 

.205 

.Og 

20D 

250 

.11 

.Og? 

21D 

400 

.125 

•0g5 

22D 

4L9 

.0g5 

2  3D 

450 

.Og 

24D 

4gO 

.Og 

25D 

510 

.09 

26D 

541 

.0g5 

27D 

571 

.095 

2gD 

600 

.17 

.0g5 

29D 

goo 

.21 

.0g5 

30D 

500 

.10 

31D 

1095 

•0g5 

32D 

965 

.Og 

3  3D 

1030 

.Og 

34D 

2000 

.33 

.095 

35D 

2500 

.42 

.10 

1.26 

36D 

3000 

.50 

.095 

1.56 

37  D 

3500 

.54 

.10 

1.5g 

3SD 

4000 

.57 

.105 

l.g 

39D 

4500 

.67 

.10 

l.g9 

40D 

5450 

.76 

.11 

2.2 

AID 

6000 

.g4 

.105 

3.14 

42D 

6500 

.g9 

.105 

3.24 

43D 

7000 

.92 

? 

3.4 

44W 

1000 

.75 

45W 

3000 

.10 

46W 

5000 

.70 

.10 

47W 

2000 

.09 

4gW 

no  wave 

recorded  due 

to  temporary 

49W 

2000 

.Og 

50W 

51W 

3000 

.09 

52W 

2100 

.35 

53W 

2000 

.33 

.0g5 
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TABLE  6  (Cont*d) 


DATA  RELATING  TO  P  WAVE  AND  EVENT  A 

(Vertical  Component) 


snogram 

No. 

Distance 

(ft) 

P  Wave 
Arrival 
(sec) 

Period 

(sec) 

54W 

2000 

.04? 

55W 

2000 

.03 

56W 

2000 

.33 

57W 

1974 

.34 

.065? 

53W 

2012 

.09 

59W 

3000 

.A3 

.03 

60W 

4000 

.53 

6lW 

5000 

.04? 

62W 

7000 

.91 

.12 

63W 

1000 

.15 

Event  "A" 


TABLE  7 


AIRBLAST  AND  GROUND  SHOCK  ARRIVAL  TIMES  AT  SHORT  RANGES 


F.E.  556  (40,000  lb) 
1962 


Air  Blast 


Ground  Shock 


Switch 

No. 

Dist . 
Feet 

Diaphragm 

Type 

Time 

m.secs. 

Piezo-Electric 

Type 

Time 

m.secs. 

Ground 

Switch 

Dist. 

(ft) 

Time 

m.secs. 

1 

7.52 

0.37 

2 

10. 25 

0.47 

1 

10.06 

X 

3 

11.95 

0.55 

4 

15.05 

0.77 

2 

15.06 

1.11 

5 

15.75 

X 

6 

17.97 

0.99 

7 

20.30 

1.19 

3 

20.05 

1.65 

5 

21.24 

1.25 

9 

25.15 

1.62 

4 

25.05 

X 

10 

27.23 

1.54 

12 

30.97 

2.25 

13 

33.96 

2.A5 

14 

35.11 

2.60? 

15 

40.94 

3.47 

5 

40.01 

X 

16 

43.74 

3.59 

6 

45.04 

4.69 

17 

46.11 

4.25 

15 

45.74 

4.69 

19 

52.55 

5.25 

21 

57.22 

6.07 

22 

60.03 

6.54 

23 

64.35 

7.33 

24 

65.90 

5.17 

7 

70.00 

9.22 

25 

73.02 

5.96 

5 

75.11 

X 

26 

75.21 

9.99 

27 

52.01 

10.79 

29 

90.94 

12.69 

30 

95.93 

13.74 

31 

99.94 

X 

9 

100.00 

X 

32 

105.25 

15.93 

10 

105.20 

17.50 

34 

115.03 

19.27 

35 

123.02 

20.75 

37 

137.90 

25.75 

35 

145.55 

25.55 

39 

150.52 

30.75 

11 

150.00 

35.06 

40 

154.94 

32.55 

12 

155.06 

37.14 

\ 

. 


- 

TABLE  7  (Gont1  d) 


Air  Blast 


Ground  Shock 


Diaphragm  Piezo-Electric 


Switch  Dist. 
Noo  Feet 

Type  Type 

Time  Time 

m.secs.  m.secs. 

Ground 

Switch 

Dist. 

(ft) 

Time 
m.secs . 

41 

42 

43 

44 

45 

46 

47 
4$ 

49 

50 

157.13 
160.27 
173.04 

135.14 
192.02 

199.91 

209.91 
220.11 
229.34 

240.15 

33.43 

34.34 

40.68 

46.61 

50.69 

54.42 

60.00 

65.90 

71.70 

77.73 

13 

250.00 

86.48 

51 

264.04 

92.45 

14 

255.07 

89.90 

52 

275.23 

99.51 

53 

282.77 

104.36 

54 

290.00 

109.06 

55 

303.24 

117 .36 

56 

339.96 

173.91 

15 

400.00 

192.62 

57 

414.76 

201.26 

16 

405.00 

194.96 

53 

459.30 

243.56 

59 

523.05 

230.01 

17 

600.00 

343.91 

60 

675.17 

401.21 

18 

605.06 

X 

61 

719.94 

437.61 

62 

399.71 

535.31 

19 

900.00 

594.76 

63 

974.35 

662.64 

20 

904.00 

X 

64 

1049.96 

703.46 

65 

1169.31 

812.60 

66 

1279.59 

906.21 

67 

1339.70 

1000.21 

21 

1400.00 

1011.56 

68 

1499.31 

1092.00 

22 

1405.52 

X 

69 

1629.31 

1230.00 

70 

1759.44 

1317.26 

71 

1369.57 

1417.56 

72 

2049.41 

1568.56 

73 

2219.57 

1716.56 

74 

2399.42 

1871.56 

75 

2749.23 

2267.56 

73 

4299.02 

3605.56 
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APPENDIX  A 


FORMAL  MATHEMATICAL  THEORY  OF  AIR  COUPLING 
A1  Original  Concept  of  ’  Air-coupled  VJaveT 

It  would  appear  that  the  study  of  air-coupled  waves  was 
introduced  by  Bateman  (1938)  who  described  an  oscillation  of  the 
surface  of  a  semi-infinite  medium  due  to  the  effect  of  a  sound  wave 
in  the  air  above  the  medium.  Bateman’s  approach  was  essentially  as 
given  below ,  though  the  notation  and  one  or  two  of  the  steps  in  the 
derivation  are  somewhat  different.  Bateman’s  original  formulation 
is  not  too  easy  to  follow,  due  to  his  introduction  of  many  abbreviating 
terns.  In  addition,  Bateman’s  development  included  the  effect  of  a 
surface  wind  and  went  on  to  discuss  the  focussing  effects  of  wind-shear. 
Thus  he  compounds  two  problems  which  are  normally  considered  as  quite 
distinct  effects  in  his  discussion  of  the  ’air-coupled  wave’.  Bateman’s 
solution  is  two-dimensional,  but  may  be  converted  to  the  case  of  axial 
symmetry.  The  axially  symmetric  case,  however,  is  an  added  complication 
without  throwing  much  additional  light  on  the  problem.  We  treat  a 
semi-infinite  half  space  in  association  with  a  semi-infinite  atmosphere, 
the  plane  of  separation  being  the  plane  z  =  0.  Let  C  be  the  velocity 
of  sound  in  the  air  (treated  as  a  constant  velocity)  and  let  and 
$  be  the  density  and  a  velocity  potential  in  the  air.  The  other  medium 
is  characterized  by  densities,  displacement  potentials  and  elastic 
constants  in  the  usual  notation.  We  wish  to  investigate  the  effect  on 
the  plane  boundary  of  the  lower  half  space  of  a  small  disturbance  in 


the  air. 


' 


- 


A  (ii) 


We  may  write  a  small  irrotational  disturbance  in  the  air  as 
satisfying : 

^  +  _  i 

dx2  dz2  c5”  St2 

(Note  -  This,  of  course,  removes  us  immediately  from  the  case  of  the 
supersonic  blast  wave) . 

We  now  look  for  solutions  of  the  type  <f>Q  =  f(z)e^<^x  ~ 

i2k2f(z)eik(x  "  ^  +  d2f(z)eik(x  "  ^ =  1  f (z)v2c2k2eik(x  "  ^ 

dz2  c2 

k2  ^4  -  f(z)  ♦  ^2f(g)  -  0 

Cc2  j  dz2 

If  we  put  f  =  AcmZ  we  get 


m  =  ±  k /I 


-  ±.  -  ±  ik  /4  - 1 


f(z)  =  A  elk^v2/02  -  +  3  ^(v^c2  -  l)'/2: 


V2 

Now  put  1— 
c2 


-  1  =  M2  and  we  get 


A  e^z  *  B  e-^z?  eik(x  ‘  vt) 
o  o  j 


A.l 


Similarly  we  may  write  for  the  potentials  in  the  ground; 


f>,  -  ^e^z  ♦  B,e-ikRz}  e11^*  ‘  vt) 
y,  =  [az^sz  *  B2eikSz  ] 

p  ? 

where  R2  =  —  -  1  and  S2  =  —  -  1,  and  we  note 
a2  |32 

that  these  potentials  are  displacement  potentials. 

Now  we  are  interested  in  a  solution  which  corresponds  to  a 


surface  w ave,  so  that  the  displacement  potentials  must  diminish 


- 


. 


, 


■ 


A  (iii) 


indefinitely  with  distance  from  the  interface.  This  means,  in  the 
usual  way,  that  the  one  constant  must  be  zero  and  the  exponential  tern 
in  the  other  part  must  have  a  real  and  negative  exponent. 

This  gives  us  immediately: 


j  ik(R„  +  x  -  vt) 

P\  =  A,e  z 


A. 2 


■tf,  -  k2e^Sz  +  x  -  Yt> 


A. 3 


where  R  and  S  are  imaginary,  that  is  v<  a,[B. 

At  the  plane  z  =  0  we  have  the  boundary  condition  that  the  velocities 
nomal  to  the  plane  must  be  the  same  in  both  media.  As  we  are  using 
a  mixed  system  of  potentials,  this  leads  to  the  simple  relationship 


dh  =  L.  l 

<9z  <3t  C_<9z  <3x  J 

which  gives 

ikM(Ao  -  B0)eik(x  -  =  -i2k2RvA,  ‘  vt)  .  i2k2vA2eik(x  ' 

or  iM(A0  -  B0)  =  kv(RA,  +  A2)  A. 4 

We  have  the  further  condition  that  there  is  no  horizontal  stress 
in  the  interface,  since  the  upper  medium  is  assumed  to  be  an  ideal 
fluid.  (Note  -  This  implies  that  there  is  no  ’drag  force’ .  In  the 
case  of  a  travelling  blast  wave  there  is  in  any  practical  case  a  very 
marked  drag  force  on  the  ground,  but  the  writer  is  not  aware  of  any 
mathematical  treatment  in  which  this  condition  is  included).  On  this 
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assumption  of  no  drag  force  we  may  write: 


2£f L  +  =  0 


<3x8  z  8X2  8z2 


which  leads  to: 


2i2k2RA,  +  i2k2A2  -  i2k2S2A2 


or,  finally. 


A2(S2  -  1)  -  2RA,  =  0 


A. 5 


We  now  assume,  as  is  relevant  in  our  case,  that  the  motion 


in  the  solid  medium  is  induced  directly  by  the  forces  deriving  from  the 
motion  of  the  air.  With  our  mixed  system  of  potentials  it  is  easier  to 
apply  this  condition  in  the  form  of  Newton* s  Law  relating  the  inverse 
effect,  that  is,  the  acceleration  in  the  air  is  assumed  to  derive  from 
the  normal  stress  in  the  ground,  thus: 


This  gives  us 


ik(x  -  vt) 


e 


i.k(x  -  vt) 


ik(x  -  vt) 


which  reduces  easily  to: 


f>Qiv(A0  +  B0)  =  XkA,(l  +  R2)  +  2pk(R2A,  +  SAg) 


A. 6 


Equations  4>  5  and  6  give  in  effect  relationships  between  the  integration 
constants,  of  which  there  are  four.  Thus  we  require  an  additional  re¬ 
lationship  before  proceeding.  The  relationship  introduced  by  Bateman 
appears  to  be  the  weakest  part  of  his  derivation,  since  he  assumes  that, 
without  loss  of  generality,  we  may  assume  there  is  some  height  H  in  the 
atmosphere  at  which  the  normal  component  of  velocity  vanishes.  It  is  by 
no  means  clear  that  this  is  acceptable,  but  one  may  assume  that  Bateman 
considered  the  disturbance  from  a  finite  source  likely  to  decay  to  small 
values  at  distances  sufficiently  far  from  the  boundary.  Alternatively, 
in  the  case  of  wind  shear  and  inversion,  which  are  included  in  Bateman’s 
analysis,  it  may  be  assumed  that  this  height  H  is  that  height  at  which 
the  ’sound  ray’  has  curved  so  that  the  tangent  to  its  path  is  parallel 
with  the  surface  and  the  ray  is  just  ’bending  down’  to  a  focus. 

Be  that  as  it  may,  if  such  a  condition  exists  we  may  write* 

A  =  Ke1™  ,  B  =  Ke"1™ 
o  0 

This  is  equivalent  to  writing: 

AQ  +  BQ  =  2K  cos  kMH 

A0  -  B0  =  2iK  sin  kMH 

Substituting  for  these  terms  we  get: 

-2KM  sin  kMH  =  kv(RA,  +  A2)  A.7 

2/°QivK  cos  KMH  =  AkA, (l  +  R2)  +  2|ik(R2A1  +  pk2)  A.g 


' 

.  :  . 
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Using  equation  5  to  substitute  for  A1  and  dividing  7  by  $  we  get: 


-M  ,  ,  ,„r  (S2  +  1)R 

— — y  tan  kMH  =  - - - : - 

\(1  +  R2)(S2  -  1)  +  2m.  [R2(32  -  1)  +  2RS} 


-mb2  tan  kMH  = 
PQiv4 


R 


A(1  +  R2)(S2  -  1)  +  2|i(R2S2  -  R2  +  2RS) 


R 


\  4  -  z\  4 +  M  -Ta  -  4  -  4 + 1  -  4  * 1  ♦ 

(-  B2  J  a2  L  a2R2  a2  R2  a. 


_  2S  2-  +  2^  v 
a2 


a232  a2  p2 


R 


v 


-  2vf  ?  +  2^  ^  _lL  -  2vf.  _  if  +  2  +  2Rs) 


2  2  2  \  /  2 o2  2  o2 

a  (3  a  J  a  p  a  p 


R 


y2  p'[^ 


2v‘ 

a‘ 


r]  *  2P,e*  { 


2KS  -  ^  +  2 


-MB2  p, 
IV  PQ 


tan  kMH 


RB* 


v4  -  4V2 B2  +  434RS  +  434 


-Mft 

iv4p 


tan  kMH  - 


R 


v4  -  4V2  B2  +  434RS  ♦  43^ 


or  finally 


M  P1  (v4  -  4V2  B2  +  434rs  +  434)  tan  kMH  -  -Riv4pQ 


A. 9 


Since  R,  S  are  imaginary ,  we  may  change  the  above  to  make  the  terns 
in  the  equation  real;  by  writing  iR  =  R1  ,  iS  =  S1 
This  gives: 

Mf^  (v4  -  4v2B2  -  4B4R1S1  +  4B4)  tan  kMH  =  R’v4^ 

Nov/  note  that  the  term 


Vr4  -  4Vj^2B2  -  AB^’S1  +  4B4  =  0  is  the  Rayleigh  Wave  condition. 


- 


'  • 


A  (vii) 


If  we  have  the  condition  that  in  our  equation  v  <  c,  and  both 
quantities  are  less  than  the  value  v^  which  satisfies  the  Rayleigh 
condition,  then  both  sides  of  the  equation  are  positive.  In  particular 
cases  therefore,  we  may  find  a  value  of  v  corresponding  to  a  given  H, 
or  alternatively,  if  we  know  the  value  of  v  we  can  determine  H  uniquely. 

It  will  be  seen  that  Bateman’s  solution,  when  written  in  the 
above  form,  is  essentially  a  derivation  for  a  Stonely  Wave  (see  Stoneley, 
1924)  in  the  case  of  a  semi-infinite  fluid  ocean  over  a  semi-infinite 
elastic  half  space.  Further,  Bateman’s  ’coupled  wave’  is  a  case  of  free 
vibrations  of  the  system,  not  the  case  of  Impulsive  forcing.  Bateman 
points  out  that  solutions  will  only  be  possible  for  certain  ranges  of 
the  elastic  constants,  and  for  rather  specialized  atmospheres.  On 
the  other  hand,  it  is  pointed  out  by  EJP  (page  112)  that  these 
’Generalised  Rayleigh  Waves’,  ’Stoneley  Waves’  or  Bateman’ s  ’Coupled 
Waves’  are  always  possible  at  the  interface  of  a  fluid  and  solid  medium, 
and  their  velocity  is  smaller  than  that  of  regular  surface  waves.  This 
agrees  with  the  Bateman  derivation. 

It  should  be  noted  that  there  is  a  fundamental  difference 
between  Bateman’s  coupled  waves  and  the  Press-Ewing  waves  discussed 
below.  The  Bateman  wave,  or  Stonely  wave,  may  exist  in  a  semi-infinite 
medium,  that  is  in  a  medium  which  is  not  normally  considered  to  be 
dispersive  to  Rayleigh  waves.  On  the  other  hand,  in  the  normal  derivation 
of  the  Stonely  wave  the  two  media  are  strictly  ’ semi- infinite’  so  that 
only  one  velocity  is  possible.  In  the  Bateman  case  a  variation  in  the 
velocity  can  occur,  depending  upon  the  variation  in  the  effective 
height  H.  Presumably  for  certain  types  of  inversion  and  wind  shear 


. 


' 
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conditions  it  would  be  possible  to  produce  a  form  of  ’dispersion’ 
with  distance  from  the  point  of  origin. 

A2  The  Press-Ewing  Air-Coupled  Wave 

The  type  of  coupled  wave  produced  by  a  constant  velocity 
pressure  step  travelling  over  the  boundary  of  a  dispersive  half  space 
is  usually  called  the  ’Press-Ewing’  air-coupled  wave,  since  the  detailed 
analysis  for  this  case  was  given  by  Press  and  Ewing  (1951) .  However, 
several  other  papers  bearing  on  similar  problems  appeared  in  a  very 
brief  span  of  time,  and  some  of  the  properties  of  the  wave  were  dis¬ 
cussed  in  papers  prior  to  the  main  source.  Reference  should  be  made, 
for  example,  to  Press,  Crary,  Oliver  and  Katz  (1951) >  which  deals  with 
a  similar  form  of  coupling  in  the  case  of  a  floating  ice  sheet,  and  to 
Haskell  (1951)  who  deals  with  the  case  where  the  ground  surface  is 
constrained  to  move  sinusoidally.  The  essential  part  of  these  dis¬ 
cussions  is  available  in  detail  in  EJP.  The  essential  difference 
between  this  treatment  and  that  of  Bateman  is  the  requirement  of  a 
dispersive  medium.  A  general  discussion  would  be  prohibitively  lengthy 
for  this  appendix,  and  attention  is  therefore  restricted  to  the  simplest 
system,  and  further  restricted  only  to  those  aspects  which  are  essen¬ 
tial  to  a  superficial  understanding  of  the  process.  As  mentioned 
previously,  we  require  as  a  minimum  a  homogeneous  atmosphere  overlying 
a  half  space  which  can  maintain  dispersive  waves.  For  simplicity,  we 
take  the  case  where  the  dispersion  is  introduced  by  an  elastic  layer 
welded  to  the  surface  of  an  otherwise  uniform  half  space.  Our  intention 
is  to  show  that  a  ’coupling  term’  enters  into  the  frequency  equation,  so 
that  at  some  frequency  a  form  of  ’resonant’  increase  in  amplitude  will 
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occur.  We  shall  use  the  notation  of  EJP,  which  is  now  more  or  less 
standard.  The  simplest  case  would  be  to  consider  a  plane  wave 
propagating  in  the  air,  but  it  is  more  realistic  to  include  the  effect 
of  a  point  source  in  the  preliminary  derivation.  Let  this  point  source 
be  at  height  h  above  the  plane  z  zero  and  let  p.f  a.,  (3j  (j  =  0,1,2) 
be  the  densities  and  wave  velocities  in  the  three  media  as  and  when 
applicable.  Let  the  elastic  layer  welded  to  the  half  space  be  of 
thickness  H.  In  order  to  write  the  distribution  of  potentials  in 
a  convenient  form,  it  is  convenient  to  divide  the  atmosphere  into 
two  zones  separated  by  the  plane  z  =  -h  (Note:  the  z  axis  is  taken 
positive  in  the  solid  medium,  negative  in  the  fluid).  Then  the  upper 
zone  will  have  a  potential  depending  only  upon  the  source  (or  alter¬ 
natively  we  may  say  dependent  only  upon  the  ’interface’  z  =  -h), 
while  other  zones  will  have  potentials  which  are  the  sums  of  two 
terms  relating  to  the  two  boundary  planes,  except  for  the  lowest 
zone,  the  half  space,  which  will  again  have  single  term  potentials. 
Further,  in  the  fluid  there  will  be  a  single  potential  only  corres¬ 
ponding  to  compressional  waves,  while  in  the  solid  there  will,  in 
general,  be  both  ’compressional’  and  ’shear’  potentials.  As  we  are 
assuming  perfectly  elastic,  homogeneous  media,  there  will  not  be, 
in  general,  any  coupling  between  the  two  modes  -  we  may  treat  the 
shear  waves  and  the  compressional  waves  as  independently  travelling 
waves,  though,  of  course,  the  displacement  at  any  point  will  be  a 


summation  of  the  two  modes. 


- 
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If  now  we  denote  by  and  ^  the  compressional  and  shear 
displacement  potentials  we  may  write  six  equations  corresponding  to 
the  four  ’distributions’  of  potentials  of  compressive  type  and  two  shear 
type,  as  below* 


<f>y0  =  Aevoz  JQ(kv) 

=  ^Bevoz  +  Ce“Vozj 

JQ(kv) 

■  [O'  ♦  O'* 

^0(kv) 

iz  -  Fe~V2z  Jo(kv^ 

Yi  =  [Meviz  +  NE"vizjjo(kv) 

Ya  =  Pe~v^z  JQ(kv) 

In  order  to  evaluate  the  nine  coefficients  we  require  nine  boundary 
conditions.  These  we  obtain  from  the  continuity  of  stress  and  dis¬ 
placement  at  the  welded  interface  z  =  H  (four  equations),  the  con¬ 
tinuity  of  vertical  displacement  and  stress  at  the  plane  z  =  0 
(two  equations),  the  vanishing  of  horizontal  stress  at  z  =  O(one 
equation)  and  an  expression  of  the  fact  that  at  the  plane  z  =  -h 
there  is  a  discontinuity  in  the  vertical  displacement  (one  equation) 
and  continuity  in  pressure  (one  equation).  The  condition  that  dis¬ 
placements  vanish  at  infinite  distance  from  the  plane  z  zero  is 
included  in  the  equations.  We  get: 
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pov^(B  +  C)  +  (2p,k2  -  P1w2)(D  +  E)  +  2p, v|k2(M  -  N)  -  0 
2vi(D  -  E)  +  (v^  +  k2)(M  +  n)  =  0 

(2^i1k2  -  P1w2)(DeviH  +  Ee_vi H)  _  (2p2k2  -  P2w2)Fe"V2H  +  2^  v|k2(MeviH  -  Ne"viH) 

+  2M.2v»k2  pe'v2H  =  0  (7  -  15) 


2(J.,  v,  (DeVlH  -  Ee-VlH)  +  2|j.2v2Fe  VaH  +  ii/vj2  +  k2)(MeV,H  +  Ne  V,H) 

"  Mv|2  +  k2)Pe"v2H  =  0 
vq(B  -  C)  -  v, (D  -  E)  -  k2(M  +  N)  =  0 


Vl(DeV]H  -  Ee"VlH)  +  v?Fe“V2H  +  k2(MeviH  +  Ne~viH  -  k2Pe~V^H  =  0 


De 


viH  +  Ee~viH  -  Fe"v2H  +  v{  (MeVl  H  -  Ne  V’h)  +  viPe-v^H  =  0 


-v!H. 


AeVoh  -  BeVoh  -  Ce-voh  =  0 

v  AeVoH  +  v  (Ce-voh)  =  2Z 
o  u 

In  the  last  equation,  Z  represents  the  strength  of  the  source, 
but  we  need  not  define  it  further  at  the  present  time.  The  last  two 
equations  enable  us  to  eliminate  B  and  C  in  terms  of  A  and  Z.  In 
theory  we  can,  following  the  method  of  Lamb,  solve  explicitly  for  the 
potentials  (and  hence  the  displacements)  for  any  particular  case, 
including  that  of  the  pulse  type  point  source,  provided  we  generalize 
the  Z  term  in  terms  of  Fourier-Bessel  integrals.  Practically  this  is 
difficult,  and  for  our  present  purpose  we  restrict  ourselves  to  the 
frequency  equation.  For  this  we  need  only  consider  the  vanishing  of 
the  determinant  on  the  LHS  of  the  foregoing  set  of  equations. 
Effectively,  we  eliminate  the  wave  amplitude  function.  We  could 
proceed  with  the  determinant  as  it  stands,  but  it  has  been  reduced 
(Jardetzky  and  Press,  1952)  to  the  form: 
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2|i,  k2  -  ^w2 
PoW2  cosh  v,H 


2v,  2k2  3?  -  w2 


0 


0 


sinH  v,H  Bf  cosh  vjH 


0  Piw2  -Pi'w2  0 

0  0  0  -  P,  w2 


Q  ~2^^2~  iW8 
Pq-w2  cpsh  VjH 


-w2 


0 


2^v{k2 


p  w2  sinh  v|H 


0 


0 


0 


0 


0  P2w2+2(ni“h>) k2  2(n2-^1)v^kt 

-^v?2  2(h2-ja,)v  P2^2+2(p1-p2)k2 


0 


0 


2v  pfcosh  v'H  p  v^sinh  vJH 
o  *  o 


0  Vjtanh  VjH  v1  coth  VjH 


0 


1  v] tanh  v’H  coth  vJH 


-1 


-k‘ 


After  purely  algebraic  (but  considerable)  manipulation,  the  determinant 
may  be  expanded  into  a  single  frequency  equation.  Details  are  given  in 
EJP,  and  to  avoid  tedium  in  presentation,  we  merely  state  that  the  final 
form  may  be  written s 

Lq  +  L,  sinhC^kH)  sinh(n3kH)  +  L2  sinhCn^H)  cosh(n3kH)  +  L3  coshCn^H) 

sinh(n3kH)  +  L4  coshCn^H)  cosh(n3kH)  =  0 
This  deceptively  simple  equation  is  complicated  by  the  fact  that  each 
of  the  L  and  n  terms  is  far  ffom  simple.  We  confine  ourselves  to 
discussing  the  L  terms,  which  include  the  effect  of  the  air.  The 
effect  does  not  enter  into  L0  but  1^  in  expanded  form  ist 
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o 


where  v  =  c/p, 

We  note  that  the  last  tern  in  this  formidable  expression  involves 
p0  and  a0,  and  thus  represents  the  effect  of  the  air.  The  air 
coupling  is  involved  through 


fo  Po 


A  similar  expansion  of  the  other  L  terms  results  in  each  case 
in  a  term  of  the  above  form.  We  see  that,  provided  the  phase 
velocity  "c"  in  the  ground  is  different  from  the  velocity  of 
T sound1  in  the  air,  then  the  coupling  term  will  be  very  small, 
due  to  the  small  value  of  the  ratio  /^/P,  .  If,  on  the  other  hand, 
the  velocity  of  the  pulse  in  the  air  approaches  the  phase  velocity 
in  the  ground,  for  some  frequency,  then  the  air-coupling  term  may 


2(H  -  1)}  2 

Hi  J 
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become  very  large,  and  possibly  dominant.  In  such  a  case  we  would 
expect  to  see  this  particular  frequency  excited  strongly  so  that  in  a 
displacement  record  at  the  surface,  the  small  amplitude  dispersed  train 
of  surface  waves  would  be  interrupted  by  a  large  amplitude  constant 
frequency  train.  This  train  is  the  T Press  Ewing’  coupled  wave.  This 
theory  is  verified  in  the  case  of  points  reasonably  remote  from  snall 
air-burst  explosive  charges.  Obviously,  in  the  case  of  the  explosions 
considered  in  this  dissertation,  the  main  departure  from  the  theory  is 
that  we  do  not  have  a  constant  value  for  the  velocity  a0,  and  so  the 
first  two  equations  in  1-6  are  not  strictly  valid.  We  require  to 
replace  them  by  a  pair  of  equations  similar  in  form  to  the  pulse  used 
by  Miles.  Analagous  to  the  above,  however,  we  should  expect  that  in  a 
given  region  (short  compared  to  the  rate  of  decay  of  velocity)  we 
should  obtain  a  coupling  effect  to  the  relevant  frequency.  We  should 
expect,  however,  that  the  particular  frequency  selected  may  (a)  change 
with  distance  from  a  given  charge  and  (b)  change  with  charge  size,  since 
the  rate  of  decay  decreases  with  increase  in  charge  size.  The  experi¬ 
mental  data  given  in  this  dissertation  verify  that  this  does,  indeed,  occur. 
In  addition,  at  points  very  remote  from  the  charge  the  effective  ’dispersion- 
in  the  coupled  wave  train  should  have  reduced  the  period  of  the  train  to 
that  predicted  by  the  above  theory  for  the  constant  velocity  pulse.  Thus 
in  the  case  of  normal  dispersion  we  would  expect  to  have  a  relatively 
high  frequency  train  at  points  remote  from  the  charge,  and  a  relatively 
long  period  train  close  to  the  charge. 
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SEISMIC  SURVEY  OF  THE  TWO  MAIN  SITES 

The  information  given  in  this  Appendix  is  based  upon  the  re¬ 
port  submitted  by  the  Century  Geophysical  Corporation  of  Canada  upon  a 
contract  completed  in  1961*  Mr  J.  L.  Robinson,  a  vice-president  of  the 
Corporation,  took  an  active  part  in  the  planning  and  interpretation,  and 
the  field  operation  was  directed  by  Mr  D.  R.  McDonald.  The  Corporation 
showed  great  interest  and  collaborated  closely  with  the  present  writer, 
undertaking  both  field  and  interpretive  work  beyond  that  called  for  in 
the  contract.  The  detailed  report  is  summarized  in  this  appendix  to 
give  only  immediately  relevant  information. 

B1  Operational  Details 

A  series  of  seismic  measurements  was  conducted  by  the 
Corporation,  party  No.  10,  on  two  sites  at  Suf field  Experimental 
Station,  on  the  week-end  of  7th  June,  1961.  The  objectives  were 
(a)  to  obtain  bedrock  profiles  beneath  the  sites,  (b)  to  identify 
ground-waves  generated  by  buried  charges  and  (c)  to  measure  the 
velocities  in  the  various  ground  waves. 

The  ground  locations,  elevations  and  bearings  of  the  two 
seianic  lines  were  provided  by  SES  surveyors.  Shot  points  and  geophone 
locations  and  elevations  were  then  obtained  by  alidade  plane  table  and 
chain  survey.  A  Mayhew  1000  Special  Drilling  Rig,  owned  by  the  Cor¬ 
poration,  was  used  for  drilling  the  shot  holes.  Drilling  conditions 
were  fair.  Very  little  gravel  was  encountered  on  the  Watching  Hill 
site.  However,  near  surface  sand  and  gravel  required  50-ft  of  casing 
to  be  set  in  VH3  and  60-ft  of  casing  in  VH4  on  the  Drowning  Ford  site. 

Recording  instruments  used  included  Century  Model  510B  ampli¬ 
fiers,  Century  Model  301BZ  6.5  cps  geophones  and  a  Century  24  trace 
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mirror- type  oscillograph  and  camera.  The  above  instruments  were  also 
connected  to  a  Fortune  Model  SR3  A.M.  tape  recorder.  The  filter  setting 
used  for  shooting  was  NF-100,  so  that  the  effective  low  filter  would  be 
the  6.5  cps  geophones.  The  AVC  setting  was  on  the  fast  position  for  all 
shooting.  Pre-suppression  was  used  where  needed.  Input  gains  varied  to 
keep  trace  interference  to  a  minimum. 

The  shooting  was  carried  in  groups  as  follows: 

Group  IA  Refraction  Profile  S.P’s  1-5  Watching  Hill 

IB  Refraction  Profile  S.P.’s  38-37  5-A  Watching  Hill 
IC  Air  shooting  6000  ft  spread  from  S.P.  38  Watching 
Hill 

IIA  5-ft  geophone  spacing  shooting  from  S.P.l  Watching 
Hill 


IIB  5-ft  geophone  spacing  shooting  from  S.P. 4  Watching 
Hill 

IIIA  Surface  3-component  phone  survey  25-ft  spacing, 
S.P.’s  11-27  Watching  Hill 

IIIB  50-ft  spacing  surface  phones  from  S.P.l  Watching 
Hill 

IIIC  Sub-surface  pressure  phone  survey  in  VH  1  and  VH  2, 
S.P.’s  11-27  Watching  Hill 
IVA  Refraction  Profile  S.P.’s  6-10  Drowning  Ford 
IVB  Short  spreading  air  shooting  from  S.P. 6  Drowning 
Ford 

VA  Surface  3-component  phone  survey  S.P.’s  28-36 
Drowning  Ford 

VB  Sub- surface  pressure  phone  survey  VH  3  and  VH  4 
S.P.’s  28-36  Drowning  Ford 
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Two  geophones  were  connected  in  parallel  for  each  station  for 
all  refraction  profiling.  Ranks  of  amplifiers  were  parallel  whenever 
possible,  to  provide  high  and  low  sensitivity  traces.  Deep  holes  con¬ 
taining  pressure  phones  were  kept  filled  with  water  during  that  phase 
of  the  programme.  Good  ground-to-geophone  coupling  was  achieved  in  the 
3-component  phase  by  first  softening  the  ground  with  water.  Explosive 
used  was  CIL  Hi- Velocity  Gelatin  60$. 

B2  Interpretation  Report,  Watching  Hill  Site 

Fig.  B1  (this  appendix)  shows  the  location  of  the  velocity  and 
refraction  survey  on  Watching  Hill  and  Drowning  Ford.  The  refraction 
survey  was  shot  for  bedrock  control.  The  velocity  survey  was  shot  to 
measure  the  velocity  and  direction  of  the  seismic  waves,  both  surface 
and  sub-surface.  Fig.  B2  shows  the  design  of  the  refraction  and  velocity 
profiles.  Reading  from  the  top  we  have  first  the  refraction  lines  with 
phones  spread  250  ft  apart  to  obtain  bedrock  control.  Secondly,  a  group 
of  stations  spaced  5  ft  apart  at  S.P.’s  1  and  4  to  obtain  dispersion  data. 
Third,  a  design  to  analyse  wave  types.  Three- component  phones  were  set 
out  near  S.P.  11  to  respond  to  shots  at  S.P.’s  11-27.  Finally,  the  velo¬ 
city  profile  was  obtained  by  lowering  the  two  cables  containing  pressure 
phones  in  VH1  and  VH2,  to  respond  to  energy  from  S.P.’s  11-27.  This 
determined  the  velocity  and  direction  of  waves  propagating  below  the 
surface. 

Fig.  B3  shows  the  refraction  profile,  together  with  the 
lithology  actually  logged  by  the  driller.  Water  was  found  at  55  ft  in  a 
thin  gravel  seam,  although  the  soft  grey  clay  above  and  below  it  appeared 
relatively  dry.  This  suggests  a  false  water  table,  since  the  stream  level 
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of  the  South  Saskatchewan  is  at  least  250  ft  below  the  Watching  Hill 
surface*. 

Fig.  B4  shows  the  refraction  time-distance  plot  obtained  from 
the  horizontal  measurements.  Here  the  velocity  distribution  is  roughly 
5500  fps,  7500  fps  and  8500  fps.  The  base  of  the  5500  fps  material  is 
considered  to  be  the  bedrock.  Fig.  B5  gives  the  data  obtained  from  the 
vertical  velocity  survey. 

In  any  given  velocity  bed  there  seems  to  be  little  evidence 
of  anisotropy  when  wave  front  directions  change. 

(Additional  data  and  interpretation  are  available  in  the  full 
report.  This  appendix  merely  abstracts  the  essential  items.  No 
detailed  interpretation  is  yet  available  for  the  Drowning  Ford  site, 
but  no  unexpected  features  have  been  noted  in  the  layering  system) . 


See  comments  in  Section  2.20.  This  water-bearing  gravel  has  fre¬ 
quently  been  broached  by  high  explosive  trials.  In  these  cases 
the  craters  filled  briefly  with  water,  the  water  subsiding  to 
leave  dry  craters.  The  bed  seems  to  tap  local  surface  water,  the 
available  volume  being  small. 
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APPENDIX  C 


THEORETICAL  RAYLEIGH  WAVE  DISPERSION  CURVES  FOR  VARIOUS  MODELS 

This  appendix  gives  the  details  of  sane  of  the  dispersion 
curves  computed  for  the  writer  on  the  Toronto  IBM  7090  by  Dr  Fraser, 
the  DRB  representative  at  the  Computation  Centre.  The  programme  used 
was  provided  to  the  writer  by  Dr  D.  L.  Anderson,  Seianology  Labs., 
Pasadena.  It  was  noted  that  the  programme  will  not  handle  input  data 
in  which  the  starting  or  ending  points  of  the  computation  are  fixed 
outside  the  asymptotic  limits,  and,  if  such  starting  or  ending  points 
are  included  considerable  machine  malfunction  occurs,  leading  to  reams 
of  nonsensical  output.  No  such  difficulty  occurs  provided  the  asymptotic 
velocities  are  not  approached  closer  than  one  phase  velocity  decrement 
(100  fps  in  the  present  case) .  The  group  velocity  calculations  are 
slightly  suspect  since  in  all  cases  the  values  computed  for  2000  fps 
and  1000  fps  are  apparently  in  gross  error,  departing  radically  from 
the  generally  smooth  curve.  The  reason  for  this  malfunction  has  not 
been  traced. 

Fig.  Cl  gives  a  pictorial  plot  of  eight  physical  models,  in 
which  the  velocity  and  density  parameters  approximate  to  those  of  the 
test  sites.  The  layering  system  of  the  first  model  (listed  as  Model  6) 
is  based  directly  upon  the  survey  data,  and  is  the  preferred  model  of 
the  Saint  Louis  group.  Model  $  modifies  this  model  at  the  high  velocity 
end  only,  and  model  9  and  model  10  similarly  retain  the  upper  layering 
constant  but  modify  the  high-velocity  system.  Model  11  introduces  two 
thin  low-velocity  layers  near  the  surface. 

The  outputs  for  these  models  are  given  as  photo-copies  of  the 
computer  output,  and  models  6  and  11  are  also  plotted  in  Fig.  C2a.  The 
points  marked  TKT  give  the  approximate  limits  of  the  computation  for 


. 
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this  model  by  Kisslinger,  using  the  Lamont  programme.  It  will  be  seen 
that  Model  11  departs  radically  from  Model  6  in  the  high  frequency 
region.,  a  maximum  shift  in  period  of  0.1  sec  at  1000  fps.  In  the  low 
frequency  region  Model  11  shifts  to  somewhat  higher  velocities  than 
those  of  Model  6.  It  should  be  noted  that  Models  8,  9  and  10,  which 
retain  the  upper  layering  of  Model  6,  conform  quite  closely  to  Model  6 
up  to  the  point  where  Models  6  and  11  cross,  and  then  conform  to 
Model  11.  Thus,  the  general  shift  towards  lower  frequencies  in  the 
range  from  2  cps  to  10  cps  between  Model  6  and  Model  11,  is  purely  the 
result  of  the  introduction  of  two  thin  low  velocity  layers  at  the  sur¬ 
face.  This  zone  is  equivalent  to  the  silty-clay  zone  above  the  water 
bearing  gravel,  and  is  liable  to  large  changes  in  water  content,  depend¬ 
ing  upon  the  precipitation  history  of  the  area.  Thus  the  probability  of 
shifts  in  the  dispersion  curves  of  the  order  shown  in  the  figure  is  quite 
large.  In  addition,  these  surface  layers  are  precisely  those  most  likely 
to  be  modified  by  the  minor  relief  in  the  plains  areas. 

The  remaining  three  models  (Models  12,  15  and  16)  indicate  the 
variation  which  may  be  expected  when  we  change  from  a  ’best  selection’ 
to  an  arbitrary  layering  system.  Model  12  uses  approximately  the  same 
parameters  as  Model  $,  but  fixes  the  layer  interval  at  50  ft.  Models 
15  and  16  are  identical  with  Model  12,  except  that  the  layering  is  cut 
short  at  the  third  and  second  layer  respectively,  the  half  space  become 
the  third  or  fourth  level  of  Model  12.  The  computer  outputs  are  in¬ 
cluded,  and  the  curves  plotted  in  Fig.C2b.  It  will  be  seen  that  the 
curves  in  C2b  are  similar  in  pattern  to  the  curves  in  C2a,  and  generally 
speaking  the  phase  velocity  curves  fall  between  those  of  Models  6  and  11. 
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Thus,  the  selection  of  a  Tbest  model’  on  the  basis  of  limited  ranges 
in  the  experimental  data  is  not  satisfactory.  If  only  a  narrow  range 
of  frequencies  is  available,  particularly  if  this  range  is  at  the 
high  frequency  end,  almost  any  model  which  is  physically  not  too  far 
removed  from  the  correct  one  will  satisfy  the  data.  Selection  of  the 
correct  model  requires  a  wide  range  of  frequencies. 

It  will  be  noted  that  in  all  the  models  the  possibility  of 
Press-Ewing-type  air-coupling  exists.  A  ’group  period’  of  between 
0.5  and  0.6  sec.,  pre-existing  in  the  ground,  will  be  interrupted  by, 
or  have  superimposed  upon  it,  a  constant  frequency  train  of  about 
0.3  sec.  period. 

Attention  is  drawn  to  the  Group  Velocity  plateau  between 
0.3  sec.  and  0.6  sec.  in  these  models. 
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(b)  THEORETICAL  DISPERSION  CURVES 
MODELS  12.  15  AND  16 
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ANISOTROPIC  RAYLEIGH  WAVE  DISPERSION 
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APPENDIX  D 


EDITED  EXTRACTS  FROM  SEiII-ANNUAL  REPORTS  ON  CONTRACT  AF  19 (604) -7402 

D1  The  contents  of  this  appendix  are  entirely  abstracted  from  the  above 
semi-annual  reports  submitted  by  Prof.  Carl  Kisslinger,  Saint  Louis 
University,  on  the  contract  ’Seismic  Waves  Generated  by  Chemical 
Explosions’  No.  AF  19(604) -7402  Project  Code  No.  $100  -  Task  2  (Vela- 
Uniforai).  The  present  writer  takes  no  credit  and  accepts  no  responsi¬ 
bility  for  the  contents,  which  are  due  to  Kisslinger,  Mateker,  and 
McEvilly  —  except  for  the  fact  that  the  records  discussed  were 
obtained  at  the  SES  test  site  under  the  writer’ s  supervision#  It  is 
felt  that  there  is  some  merit  in  including  these  verbatim  transcriptions 
herein  in  order  to  make  the  dissertation  as  comprehensive  as  possible. 

D2  Semi-Annual  Report  No.  2 

The  report  starts  with  a  detailed  description  of  the  partici¬ 
pation  at  SES.  Then  follows  a  general  description  of  the  results, 
which  continues* 

"The  amplitude  of  the  first  arrival  at  2000  ft  on  the  vertical 
is  about  300  microns,  with  a  pulse  duration  of  0.037  second.  This  was 
determined  to  be  a  head  wave  from  the  top  of  the  hard  clay  with  a  velo¬ 
city  of  5600  fps.  The  relatively  small  radial  component  agrees  with 
the  interpretation  that  this  is  primarily  refracted  energy  coming 
from  below.  The  maximum  peak-to-peak  motion  is  about  4  mm,  with  a 
period  of  0.5  second,  corresponding  to  an  acceleration  of  about 
0.03  g.  The  maximum  motion,  as  at  all  distances,  occurs  after  the 
arrival  of  the  airblast . 


' 
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Transverse  motion,  of  particular  interest  to  the  (Saint 
Louis)  investigators  is  small,  with  a  maximum  amplitude  in  the  body 
waves  of  $0  microns,  or  about  0.27  of  the  maximum  P  wave.  The 
maximum  peak  to  peak  SH  amplitude  is  754  microns,  20$  of  the  maximum 
vertical  motion.  At  4000  ft  the  compressional  waves  arriving  by 
various  paths  are  well  separated.  The  first  arrival,  which  can  be 
correlated  over  the  rest  of  the  profile  to  15,000  ft,  gives  a 
velocity  of  8600  fps.  The  event,  which  was  the  first  arrival  at 
2000  ft,  is  a  prominent  event  following  the  onset  by  0.13  sec. 

.  the  first  motion  at  10,000  ft  is  noteworthy.  The 

period  is  0.18  sec.,  and  the  peak  to  peak  amplitude  of  27  microns 
is  greater  than  the  extension  of  the  attenuation  curve  for  the  first 
7000  ft  would  predict.  At  12,000  ft  and  15,000  ft  the  period  becomes 
0.3  sec.  The  apparent  surface  velocity  is  very  high,  about  50,000  fps. 
If  we  accept  a  value  of  20,000  fps  in  the  Mississippian,  a  dip  of 
about  15  degrees  is  required  to  reconcile  the  values.  Since  the 
profile  could  not  be  reversed  there  is  no  confirmation  of  this 
interpretation. 

(At  4000  ft  and  6000  ft  the  transverse  motion) .  has 

separated  into  three  distinct  events,  all  of  which  can  be  correlated , 
though  with  difficulty  because  of  the  small  amplitude, to  12,000  ft. 
The  ratio  of  the  maximum  SH  to  maximum  Vertical  P  is  .26  at  4000  ft 
and  .23  at  12,000  ft.  The  velocities  of  3&40  and  5060  fps  correspond 
with  5600  and  8600  fps  for  the  P  waves  in  the  same  refractor.  This 
clear  cut  development  of  refracted  SH  waves  from  the  explosive  source 
is  of  considerable  interest. 


. 
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If  the  $500  fps  arrival  is  a  shear  wave,  the  corresponding 
P  wave  would  have  a  velocity  of  about  14*000  fps.  No  such  arrival  could 
be  identified  on  the  verticals.  Since  the  velocity  and  time  intercept 
of  this  event  are  in  close  agreement  with  the  $500  fps  event  on  the 
verticals,  and  its  amplitude  is  small  compared  to  the  other  SH  arrivals, 
it  is  likely  that  this  is  actually  a  component  of  a  P  wave  returning 
from  a  dipping  refractor,  so  that  it  has  a  small  component  not  in  the 
vertical  plane  through  the  profile . 

The  surface  wave  motion  preceding  the  airblast  is  very 
similar  on  lines  A  and  B.  The  motion  on  Line  C  is  slightly  larger. 

The  apparent  time  delay  in  the  arrivals  at  this  station  may  be 
explained  by  the  fact  that  the  air  blast  travel  times  have  not  been 

* 

corrected  for  the  slight  wind  blowing  in  the  direction  of  the  C  line  . 
The  motion  over  the  first  two  seconds  is  in  phase  on  all  three,  but 
the  event  arriving  at  5  sec.  is  reversed  in  phase  on  the  C  line  record 
as  compared  to  the  other  two.  This  indicates  that  the  two  long 
period  events  have  different  origins,  one  of  which  has  a  node  between 

the  B  and  C  lines . .  The  presence  of  refracted  SH  waves,  about 

one-fourth  as  large  as  the  corresponding  P  waves  is  of  interest  as  it 
adds  evidence  that  SH  motion  must  be  expected  as  a  normal  result  of 
an  explosion  in  or  on  the  earth. 

Small  Scale  Tests  at  SES 

As  part  of  the  effort  to  complete  work  on  explosion 
generated  SH  waves,  it  was  desired  to  extend  the  observations  into 

*  See  discussion  in  para  3-60.  The  above  is  incorrect. 
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new  media . the  experiments  were  designed  to  be  of  the  same  type 

as  those  carried  out  previously  for  clay  sandstone  and  limestone. 
However,  the  availability  of  ten  seismographs  made  possible  more 
complete  observations  than  in  any  of  the  earlier  work.  The  layout 
was  completed  by  the  Saint  Louis  team. . .but .. .the  work  was  carried 
out  by  the  SES  staff . 

The  analysis  of  the  results  (at  Saint  Louis)  has  been 
directed  to  the  buried  charges.  The  surface  and  air-bursts  yielded 
records  that  are  markedly  different,  and  open  a  new  area  of  enquiry. 

The  transverse  motion  was  derived  from  the  seianograms . . .and  the 
presence  of  clear-cut  nodes  with  polarity  reversal  on  either  side 
is  evident.  The  presence  of  high  frequency  small  amplitude  motion 
near  the  nodes  is  also  interesting .. .The  results  confirm  the  previous 
conclusion  of  unsyrametrical  radiation  of  SH  motion. 

Another  interesting  result  from  the  same  set  of  records 
concerns  the  P  wave  radiation.  The  pattern  shows  a  bulge  corresponding 
to  an  SH  node. ... .according  to  Knoppoff  and  Gilbert’s  model  of  crack 
propagation  there  should  be  a  lobe  of  P  energy  in  the  direction  of 
propagation. . . .The  amplitude  of  the  long  period  motion  for  the  shallow 
shot  was  found  to  be  less  than  from  the  deeper  shot... which  contradicts 
the  Florissant  results. .. .and  is  obviously  connected  with  the  generating 
mechanism. .. .the  surface  wave  portion  of  the  longitudinal  and  vertical 
traces  for  the  10  lb  shot  at  1$.5  ft,  recorded  from  60-&00  ft  have 
been  examined.  Rayleigh  type  motion  is  well  developed  with  an  early 
prograde  event  showing  little  dispersion,  and  a  prominent  packet  of 
Rayleigh  waves.  The  period  of  the  train  is  fairly  constant  at  about 
0.20  sec... and  the  maximum  amplitude  moves  through  the  packet  with 


' 


a  group  velocity  of  about  310  fps.  The  observed  phase  velocities  range 
from  600  to  BOO  fps ...... 

D3  Semi-Annual  Report  No.  3 

...  Results  obtained  at  SES,  and  by  other  investigations, 
indicate  that  the  amount  of  energy  in  the  seismic  waves  changes 
radically  as  a  shot  depth  is  varied  from  just  below  to  just  above... 

The  spectra  of  the  entire  radial  motion  at  about  15  metres 
from  shots  in  sandy  clay  (SES)....  were  calculated.  In  the  sandy  clay 
at  SES  the  spectra  peaked  at  A  cps,  the  predominant  Rayleigh  wave 
frequency,  with  secondary  peaks  at  6-7  cps  and  9  cps  corresponding 
to  the  prograde  Rayleigh  waves  and  possibly  body  waves  respectively. 

...  Fourier  spectra  of  the  horizontal  transverse  motion 
(Love  wave)  have  been  calculated  for  the  profiles  at  SES.  Charges 
of  10  lb  at  a  depth  of  lB  ft  and  7»5  ft,  and  a  500  lb  shot  at  the 
surface  have  been  compared.  For  the  shot  at  lB  ft  the  energy  in  the 
wave  train  is  contained  in  the  2-11  cps  band  with  prominent  peaks 
at  3-5  and  5  cps. 

For  the  shot  at  7.5  ft,  the  general  energy  distribution  is 
similar  to  the  deeper  source  with  the  exception  that  the  peaks  at  3.5 
and  5  cps  no  longer  stand  out  above  the  overall  energy  background. 
Record  amplitudes  for  this  shot  are  visibly  smaller  than  for  lB  ft 
depth. 

Spectra  from  the  500  lb  shot  are  very  similar  in  shape  to 
those  of  the  lB  ft  shot  with  the  significant  features  occurring  at 
frequencies  1  cps  lower  for  the  surface  shot. 
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The  transverse  components  of  the  three  sets  of  azimuthal 
records  (2  lb  at  10.5  ft,  2  lb  at  5  ft  and  4  lb  at  5  ft  recorded 
at  100  ft)  were  digitized  and  the  corresponding  radiation  patterns 
for  the  first  large  amplitude  motion  were  plotted.  The  position  of 
the  nodes  was  found  to  be  about  the  same  for  the  three  sets  of  data. 

Spectral  analysis  of  surface  wave  portions  of  the  records... 
have  been  completed.  At  SES  dispersed  Rayleigh  waves  and  Love  waves 
exist,  as  well  as  a  Rayleigh- type  prograde  wave.  Ultimately,  these 
data  will  be  used  to  contribute  the  necessary  parameters  for  Fourier 
re-synthesis  of  the  source  motion. 

Attenuation  coefficients  based  upon  the  relationship 
Ar  =  AqT-1  //2e_ar 

were  determined  for  individual  Fourier  components.  At  SES  a  varies 
from  0.004/meter  to  0.006/meter  in  the  range  3.5  to  5  cps . 

One  of  the  most  significant  features  of  the  SES  data  that 
as  yet  is  unexplained  completely  is  the  absence  of  the  normal  4  cps 
Rayleigh  wave  associated  with  the  low  velocity  surface  layer  from  the 
100  ton  shot.  This  event,  usually  contributing  the  biggest  motion 
on  the  vertical  and  radial  records  for  a  buried  charge,  could  not 
be  found  in  the  100  ton  records.  The  characteristic  of  the  Rayleigh 
wave  are  known  from  an  analysis  of  two  buried  10  lb  charges  and  a 
500  lb  surface  charge.  The  event  presents  an  almost  identical 
appearance  for  these  three  sets  of  records,  so  it  is  not  the  fact 
that  the  100  ton  charge  was  surface-burst  that  accounts  for  the 
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An  inspection  of  the  SES  records  from  earlier  trials 
indicates  that  the  Rayleigh  wave  appears  for  an  all  shots,  but  not 
large  ones,  of  5  tons  or  greater.  There  is  Rayleigh- type  motion 
following  the  airblast  arrival,  but  it  is  much  too  early  to  be  the 
same  event  as  that  recorded  from  the  smaller  shots. 

The  Rayleigh  wave  observed  in  the  snail  scale  shots  has  a 
very  narrow  spectrum,  a  few  cps  wide,  sharply  peaked  at  4  cps,  with 
a  group  velocity  at  this  frequency  of  about  460  fps.  One  possible 
explanation  of  the  absence  of  the  wave  is  that  the  larger  shots, 
including  the  100  ton,  do  not  generate  appreciable  energy  at  this 
frequency.  This  point  is  worthy  of  further  study.  The  Rayleigh 
waves  that  are  observed  for  the  larger  shots  are  probably  associated 
with  deeper  interfaces  with  higher  velocities.  In  the  case  of  the 
100  ton  shot,  the  entire  thickness  of  the  low  velocity  layer  was 
involved  in  the  source  region,  so  that  the  input  impulse  was  applied 
directly  to  the  top  of  the  hard  clay.  This  would  account  for  the 
generation  of  the  higher  speed  waves  not  seen  in  smaller  sources". 
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A  MODIFIED  SATO  PROGRAMME  FOR  DETERMINING  PHASE  VELOCITIES 

(Hastrup  Programme) 

At  some  stage  in  the  analysis  of  seismograms  it  becomes  necessary 
to  investigate  the  spectra  and  the  dispersion  curves  for  the  various  pulses. 
The  most  elegant  way  of  doing  this  is*  without  doubt,  some  modification  of 
the  method  of  Sato  (1955) >  but  the  precise  form  in  which  the  method  is 
applied  will  depend  upon  the  type  of  record  and  the  particular  computing 
facilities  available.  In  this  appendix,  methods  of  determining  the  spec¬ 
tra  and  the  dispersion  curves  which  are  peculiarly  suited  to  the  Suffield 
records  are  described.  The  formulation  is  due  to  0.  Hastrup,  of  the  Danish 
Defence  Research  Board,  and  was  prepared  during  his  tenure  of  a  DRB  Fellow¬ 
ship  at  SES  during  the  summer  of  1963.  Examples  of  the  computations  using 
these  programmes  are  given  earlier  in  this  dissertation,  and  will  also 
appear  in  greater  detail  for  specific  cases  elsewhere  (Hastrup,  1963a,  b) . 
Numerical  Calculation  of  Fourier  Integrals 


The  Fourier  Integral  is  defined  by 


E.l 


where  f(t)  is  the  transient  function, 

F(u>)  is  generally  a  complex  function  of  the  frequency  c o  and 


can  be  written  either  as 


E.2 

E.3 


or  as 


F  («)  =  |F(o3)|  e16^ 


where 


and 


*  The  factor  ^ .  is  i 
v/2n 


_  is  introduced  to  make  the  Fourier  integral  and  the 


inverse  transformation 


form 


. 
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E  (ii) 


0(a))  =  arc  tan 


ReCF(")H 


E.4 


If  f(t)  is  a  complicated  function,  or  is  only  known  graphically, 
it  is  necessary  to  evaluate  E.l  numerically.  The  present  method  is  based 
upon  linear  interpolation  of  the  function. 

Consider  the  signal  between  t  ^  =  (n  -  l)a.t  and  tn  =  n&t;  then 
f(t)  in  this  region  may  be  approximated  by 

f(t)  RJ  g(t)  -  fn_,  -  6n(-t  -  n  +  1)  E.5 


where  6n  =  f  -  f  _  . . 
n  n  n-1 


The  Fourier  Integral  for  this  small  trapezoid  is  then  given  by 

k"  ■  ■*  -icot 


AF"(W)  '  vk  C,  ^  *  6n(^ ' n  *  ^  e_i 


dt 


E.6 


This  may  be  re-written  as 

Re[>Fn(^  =  —=  ^  Esin  y(fn  cos(2n  -  l)f  sin(2n  -  If)  + 
v2n 

&n  sin  2(n  -  l)  ijr\ 


Im]jiFn(a))3|  =  -y  E  "sin  "  l)f  +  ~  cos(2n  -  l)?j  + 


cos  2(n  - 

2 


coat 


where  we  have  introduced  the  angle  Y  =  2  * 

To  get  f(o))  we  have  to  sum  (6)  over  all  the  intervals  to  get 

N  N 

F(o))  =  £  ReEAFn(w)E  +  i  E  In£^Fn(a))E 
n=1  n=1 


E.7 


. 
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The  input  tape  for  this  calculation  on  the  SES  Stantec 


Zebra  computer  is  given  below  in  Zebra  Simple  Code. 


+  <oo 

2IT 

(“o 

Initial  frequency) 

+  AW 

(  AtO 

Frequency  increment) 

2TI 

+  ^m 

2n 

(  ^ 

'  m 

Majcimum  frequency) 

Y 

Blank  tape 

+  Identification  Number 
+dt 

+  h  (bias  level) 


Blank  tape 
+  f. 


+  f 


1  C  First  record  ordinates 


+  f 


N 


Y1Y00 

+  Ld o 

2n 

+  AU 

2n 

+  f£gi 

2n 

Blank 

+  Iden.  No. 
+  at 
+  h 
Y 

Blank 

+  *o 

+ 


Route  Selection  by  Key 
Stop  if  Y1Y00 

Stop  after  reading  new  Iden.  No 
and  data. 

Punch  F(a>)  and  B(^) 

Punch  R^F(6a)3  and  ImQ’(^)II 


Programme  may  be  re-entered  at 
punch  by  dialling  Q14.  U4  then 
selects  output. 


+  Iden.  No. 
+  at 

+  h 
Y 

Blank 

+  fo 
+  fi 


+  % 

Y 

Blank 

YY11Y00  (STOP) 


Blank 

YY11YOO  (STOP) 


No  purpose  would  be  served  by  listing  the  actual  computer  programme.  The 
output  from  this  calculation  may  then  be  used  as  the  input  into  the  dis¬ 
persion  programme. 


APPENDIX  F 


SPECTRAL  ANALYSIS  OF  POST  AIR  BLAST  PHASE  OF  TWO  SEISMOGRAMS 

In  this  appendix  is  given  an  example  of  the  method  of  detailed 
reduction  which  must  be  applied  to  all  the  individual  phases  in  the 
seismograms  before  complete  interpretation  is  possible.  The  example  is 
chosen  as  illustrating  both  the  power  and  the  weakness  of  the  method. 

It  was  decided  to  attempt  an  analysis,  post  airblast  phase, 
of  three  seismograms  from  the  Watching  Hill  40,000  lb  charge,  F.E.  535* 
seismograms  60W,  6lW  and  62W  being  selected  as  apparently  consistent 
records  showing  air-coupling.  The  analysis  was  undertaken  without 
reference  to  previous  analyses  of  these  records  or  associated  records, 
in  order  to  illustrate  the  complete  approach.  The  analysis  therefore 
proceeded  in  the  following  stages. 

Stage  1.  Digitization 

The  records  were  digitized  onto  punched  tape  including  the 
following  data*,  identification  number,  charge  size,  epicentral  distance, 
bias  from  base  line  to  the  transverse  trace,  vertical  trace,  and  radial 
trace,  at  the  zero  position;  number  of  counts  at  0.02  second  intervals 
from  zero  to  the  start  of  digitization  (in  this  case,  the  arrival  of  the 
airblast) |  digital  data  in  sets  of  three,  giving  transverse,  radial  and 
vertical  in  that  order,  to  a  point  well  beyond  the  phases  of  interest; 
total  number  of  counts. 

Stage  2.  Verification 

The  punched  paper  tape  from  stage  1  was  used  as  the  input 
to  a  plotting  programme  on  the  Stantec  Zebra  computer.  Errors  in 
digitization  could  thus  be  detected  by  direct  comparison  between  the 
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original  trace  and  the  plotted  trace,  or  by  machine  failure  in  the 
event  of  gross  error.  The  computer  plotted  traces  for  the  vertical 
components  of  the  three  seismograms  are  given  as  Fig.  F.l. 

Stage  3.  Particle  Trajectories 

From  the  corrected  digital  data,  particle  trajectories 
were  plotted  for  the  immediate  post  airblast  region.  These  tra¬ 
jectories  are  given  as  Figs.  F2,  F3  and  F4.  It  was  immediately 
apparent  that  the  trajectories  were  not  completely  consistent  with 
each  other.  The  trajectory  for  seismogram  62W  is  grossly  different 
in  appearance  from  the  other  two,  and  can  not  be  correlated  with  them. 
On  the  other  hand,  the  orbits  from  60W  and  6lW  showed  a  basically 
similar  pattern,  although  there  are  differences  even  in  the  direction 
of  rotation  of  regions  of  the  orbits.  The  assumption  was  made  that 
comparison  was  possible  from  the  moment  of  arrival  of  the  airblast 
to  a  distinct  transition  in  the  orbits.  These  starting  and  ending 
points  were  selected,  and  are  indicated  in  the  figures. 

Stage  4.  Smoothing  of  Traces 

It  is  essential  that  the  analysis  be  performed  upon  phases 
which  start  and  end  at  zero  amplitude.  The  vertical  component  digital 
data  for  the  selected  region  were  therefore  abstracted  from  the  punched 
tape,  and  a  fresh  tape  made  in  which  the  first  and  last  few  values  were 
corrected  or  ’ smoothed1  so  that  the  first  and  last  reading  in  the 
smoothed  data  became  equal  to  the  bias  level  for  the  vertical  trace. 

Stage  $.  Spectral  Analysis 

The  smoothed  data  from  stage  4  were  then  used  as  input  data 


into  the  Hastrup  spectral  analysis  programme  of  Appendix  E.  The  output 
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data  for  the  two  seisnograras  are  tabulated  as  Fig.  F.5.  These  data 
were  not  used  directly,  but  used  as  an  input  to  stage  6. 

Stage  6.  Dispersion  and  Spectrum  Programme 

The  output  tapes  from  stage  6  were  used  as  the  input  data 
for  the  Hastrup  programme  of  a  Sato- type  analysis  (see  Hastrup,  1963b). 
This  programme  computes  firstly  the  spectral  amplitude  for  the  two 
seismograms  at  selected  frequencies  (in  this  case,  at  intervals  of 
0.25  cps  from  1  cps  to  6  cps).  Finally,  it  computes  the  family  of 
dispersion  curves,  provided  the  spectral  distribution  is  consistent 
between  the  two  seismograms.  The  output  for  the  first  part  of  this 
analysis  is  given  in  Fig.  F.6,  and  plotted  in  Fig.  F.y.  It  will  be 
seen  that  the  spectra  are  similar  in  pattern,  with  a  large  amplitude 
peak  in  the  vicinity  of  2  cps,  and  a  smaller  amplitude  peak  in  the 
vicinity  of  3.5  cps,  followed  by  several  minor  peaks  of  small  signi¬ 
ficance.  Visual  inspection  of  these  seismograms  immediately  suggests 
that  a  Sato-type  dispersion  analysis,  which  assumes  that  the  spectral 
distribution  remains  constant  between  stations,  will  not  give  a  con¬ 
sistent  picture  in  this  case.  The  actual  computation  verifies  this 
since  the  families  of  velocity  data  for  each  frequency  of  the  com¬ 
puted  values  are  negative  in  most  cases,  Wiich  is  physically  meaning¬ 
less. 

Stage  7.  Interpretation 

Interpretation  of  a  single  analysis  of  the  above  type  is, 
obviously,  a  dangerous  proceeding.  Nevertheless,  certain  deductions 
may  be  made.  In  the  first  place,  it  is  noted  that  the  mean  peak 
frequencies  for  the  two  spectra  are  2  cps  and  3*5  cps  for  the  two 
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dominant  peaks,  most  of  the  energy  being  in  the  2  cps  region.  Comparison 
with  the  theoretical  dispersion  curve  of  Model  6,  Appendix  C,  shows  that 
these  two  frequencies  correspond  exactly  with  the  sonic  velocities  in  the 
group  and  phase  dispersion  curves  respectively,  which  is  what  one  would 
expect  if  the  post  airblast  region  is  composed  of  a  Press-Ewing  air- 
coupled  wave  superimposed  upon  a  pre-existing  surface  wave.  Despite 
this  precise  agreement,  attention  is  drawn  to  the  fact  that  the  pre¬ 
existing  wave  is  a  prograde  wave,  not  a  retrograde  Rayleigh  wave. 

If  we  interpret  the  seisnograms  in  this  way,  that  is  as  an 
air-coupled  wave  superimposed  upon  a  pre-existing  wave,  the  appearance 
of  the  trajectories,  and  the  impossibility  of  applying  Sato’s  method  to 
compute  the  phase  dispersion  curve  may  be  explained  readily.  The  precise 
appearance  of  the  orbit  will  depend  upon  the  individual  orbits  of  the  two 
components,  and  changes  in  the  relative  phase  of  these  two  components 
may  modify  the  composite  orbit  to  such  an  extent  that  changes  of  the  order 
shown  in  Figs.  F2  and  F3  may  occur. 

An  analysis  of  the  above  type  is  a  comparatively  slow  process. 
Digitization  of  the  records  takes  approximately  six  hours  for  a  pair  of 
seismograms.  Verification  of  the  digitization  requires  approximately 
30  minutes  of  computer  time  followed  by  about  two  hours  of  teleprinter 
time.  Plotting  of  the  particle  trajectories,  after  correction  of  the 
data  amounts  to  roughly  an  hour  per  trajectory.  The  spectral  analysis 
programme  takes  approximately  three  hours  for  records  of  the  length  used 
in  the  present  case.  The  final  stage,  computing  the  spectral  amplitude 
and  the  phase  velocities  takes  about  30  minutes  of  computer  time.  The 
time  schedule  given  assumes,  of  course,  that  we  are  dealing  with  clearly 
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separated  phases  which  may  be  correlated.  Where  it  becomes  necessary 
to  use  subjective  judgment  in  selecting  the  data  to  be  analysed,  the 
time  involved  may  be  increased  by  a  large  factor. 


COMPUTER  VERIFICATION  OF  DIGITIZATION  -  VERTICAL  COMPONENT 


FIG.  F.  2 

POST  AIRBLAST  ORBIT  —  60  W 


FIG.  F.  3 

POST  AIRBLAST  ORBIT  —  61  W 
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FIG.  F.  4 

POST  AIRBLAST  ORBIT—  62 W 
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NUMBER  OF 

TRAVEL 

DISTANCE 

FREQUENCIES 

TIME 

+0.600000000  + 

2 

+0.210000000 

+ 

2 

+0.334000000 

+ 

1 

+0.400000000  + 

4 

ANGULAR  FREQUENCY 

FREQUENCY  CPS 

REAL  PART 

IMAGINARY  PART 

+0.628318530  + 

1 

+1.000000000 

+ 

0 

-0.111980497 

+ 

2 

+0.212992636  + 

2 

+0.785398163  + 

1 

+0.125000000 

+ 

1 

+0.364943949 

+ 

2 

+0.336831966  + 

2 
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1 

+0.150000000 
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2 
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2 

+0.175000000 
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1 

-0.642616815 

+ 

2 

-0.255177842  + 

2 
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2 

+0.200000000 
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1 

-0.171778816 

+ 

2 

+0.109274870  + 

3 

+0.141371669  + 

2 

+0.225000000 
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1 

+0.143919460 

+ 

3 

+0.392652289  + 

2 

+0.157079632  + 

2 

+0.250000000 

+ 

1 

+0.653531312 

+ 

2 

-0.1  3071  3925  + 

3 

+0.172787595  + 

2 

+0.274999999 

+ 

1 

-0.864065750 

+ 

2 

-0.515563474  + 

2 

+0.188495559  + 

2 

+0.299999999 

+ 

1 

-0.329348742 

+ 

1 

+0.619080966  + 

2 

+0.204203522  + 

2 

+0.324999999 

+ 

1 

+0.700658625 

+ 

2 

-0.292614352  + 

2 

+0.219911485  + 

2 

+0.349999999 

+ 

1 

-0.223402680 

+ 

2 

-0.726944982  + 

2 

+0.235619448  + 

2 

+0.374999999 

+ 

1 

-0.440285867 

+ 

2 

+0.600213594  + 

1 

+0.251327411  + 

2 

+0.399999999 

+ 

1 

+0.107742141 

+ 

2 

+0.794781424  + 

1 

+0.267035375  + 

2 

+0.424999999 

+ 

1 

-0.395627615 

+ 

1 

-0.238077516  + 

2 

+0.282743338  + 

2 

+0.449999999 

+ 

1 

-0.186687164 

+ 

2 

-0.158362856  + 

1 

+0.298451301  + 

2 

+0.474999999 

+ 

1 

+0.516252486 

+ 

0 

-0.399059385  - 

1 

+0.314159264  + 

2 

+0.499999999 

+ 

1 

-0.884915091 

+ 

1 

-0.112995775  + 

2 

+0.329867227  + 

2 

+0.524999998 

+ 

1 

-0.130053843 

+ 

2 

+0.419438448  + 

1 

+0.345575191  + 

2 

+0.549999998 

+ 

1 

+0.170163987 

+ 

1 

+0.304274347  + 

1 

+0.361283154  + 

2 

+0.574999998 

+ 

1 

-0.398101381 

+ 

1 

-0.613955907  + 

1 

+0.376991117  + 

2 

+0.599999998 

+ 

1 

-0.698330204 

+ 

1 

+0.254067066  + 

1 

5  RECORD 

NUMBER  OF 

TRAVEL 

DISTANCE 

FREQUENCIES 

TIME 

+0.610000000  + 

2 

+0.210000000 

+ 

2 

+0.419000000 

+ 

1 

+0.500000000  + 

4 

ANGULAR  FREQUENCY 

FREQUENCY  CPS 

REAL  PART 

IMAGINARY  PART 

+0.628318530  + 

1 

+1.000000000 

+ 

0 

-0.221907582 

+ 

2 

-0.157205634  + 

2 

+0.785398163  + 

1 

+0.125000000 
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1 
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2 
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+0.150000000 
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+0.175000000 
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+0.349999999 
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1 

+0.213632857 
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2 
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2 

+0.235619448  + 
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1 
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2 
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+0.399999999 
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1 
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+ 

2 
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FOURIER  INTEGRAL  COMPUTATION  60W  AND  6lW 


. 


ANGULAR  FREQUENCY 


•+0.628318530  4-  l 
40.785398163  4-  l 
40.942477794  4-  1 
40.109955743  +  2 
40.125663706  4-  2 
40.141371669  4-  2 
40.157079632  4-  2 
40.172787595  4-  2 
40.188495559  4-  2 
40.204203522  4-  2 
40.219911485  4-  2 
40.235619448  4-  2 
40.251327411  4-  2 
40.267035375  4-  2 
40.282743338  4-  2 
40.298451301  4-  2 
40.314159264  4-  2 
40.329867227  4-  2 
+0.345575191  +  2 
40.361283154  +  2 
40.376991117  +  2 


SPECTRAL  AMPLITUDE 


60W 


61W 


FREQUENCY 

CPS 


4-1.000000000  4-  0 
40.125000000  +  1 
40.150000000  +  1 
40.175000000  +  1 
40.200000000  4-  1 
40.225000000  4-  l 
40.250000000  +  1 
40.274999999  +  1 
40.299999999  +  1 
40.324999999  +  1 
40.349999999  +  1 
40.374999999  +  1 
40.399999999  +  1 
40.424999999  +  1 
40.449999999  +  1 
40.474999999  +  1 
40.499999999  +  1 
40.524999998  +  1 
+0.549999998  +  1 
+0.574999998  +  1 
+0.599999998  +  1 


40.240635605  +  2 
+0.496628492  +  2 
+0.593841660  +-  2 

+0.691427582  +  2 
+0.110616802  +  3 
+0.149179654  +  3 
+0.146140898  +-  3 
+0.100618851  +  3 
+0.619956409  +  2 
+0.759306044  +  2 
40.760498366  +-  2 
+0.444358198  +  2 
+0.133884816  +-  2 
40.241342321  +  2 
+0.187357640  +  2 
+0.517792539  +  0 
+0.143522794  +  2 
+0.136650240  +  2 
+0.348623953  +  1 
+0.731728478  +-  1 
+0.743111800  +-  1 


+0.271949602  +-  2 
40.607230239  +  2 
+0.867287455  +  2 
+0.964653818  +-  2 
+0.909089996  +■  2 
+0.787094652  +  2 
+0.654494927  +  2 
+0.484800115  +  2 
+0.300880112  +  2 
+0.294523397  +  2 
+0.391116722  +-  2 
+0.375803941  +  2 
+0.242038139  +  2 
+0.117009452  +-  2 
+0.148545889  +  2 
+0.155689700  +  2 
+0.904938095  +  1 
+0.742496653  +  0 
+0.412059892  +-  1 
+0.379384479  +  1 
+0.599764562  +  0 
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SPECTRAL  AMPLITUDES  SEISMOGRAMS  60W  AND  61W 
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SPECTRA  60  W  AND  61  W 


